Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 932,340 (A)8→7 intergenic (‑290/‑255) trxB ← / → lrp thioredoxin reductase, FAD/NAD(P)‑binding/leucine‑responsive global transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913932,3330A.96.4% 221.1 / 1.1 55intergenic (‑283/‑262)trxB/lrpthioredoxin reductase, FAD/NAD(P)‑binding/leucine‑responsive global transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (27/26);  minor base C (1/1);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGTCTTCT  >  NC_000913/932196‑932465
                                                                                                                                         |                                                                                                                                    
cTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGgaagaa                                                                                                                                     >  1:953020/1‑139 (MQ=255)
cTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGgaagaa                                                                                                                                     >  1:891323/1‑139 (MQ=255)
     tAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                >  2:272577/1‑139 (MQ=255)
     tAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                <  2:161353/139‑1 (MQ=255)
     tAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                <  2:810812/139‑1 (MQ=255)
        aTTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaacag                                                                                                                             >  1:637715/1‑136 (MQ=255)
         tttCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaacagt                                                                                                                            <  1:839230/139‑5 (MQ=255)
             ccGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTc                                                                                                                       >  1:514161/1‑139 (MQ=255)
             ccGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTc                                                                                                                       >  1:565500/1‑139 (MQ=255)
                 tttATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTtat                                                                                                                   <  1:8786/139‑1 (MQ=255)
                   tataAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTatat                                                                                                                 <  1:989581/139‑1 (MQ=255)
                           gATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATaa                                                                                                         >  1:76942/1‑139 (MQ=255)
                           gATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATaa                                                                                                         >  2:47139/1‑139 (MQ=255)
                            aTTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              <  2:690517/118‑3 (MQ=255)
                            aTTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              >  1:690517/1‑116 (MQ=255)
                             ttAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                <  1:858114/115‑1 (MQ=255)
                             ttAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                >  2:858114/1‑115 (MQ=255)
                              tAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACca                                                                                                      <  2:514161/139‑1 (MQ=255)
                                aaTGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACcatg                                                                                                    <  2:255122/139‑1 (MQ=255)
                                aaTGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACcatg                                                                                                    <  1:113364/139‑1 (MQ=255)
                                    atgaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACcatgcatg                                                                                                >  2:113254/1‑139 (MQ=255)
                                      gaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                <  1:624399/106‑1 (MQ=255)
                                      gaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                >  2:624399/1‑106 (MQ=255)
                                      gaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTa                                                                                              >  1:22596/1‑139 (MQ=255)
                                       aaTAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGAAAATCCACTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTATTATATGCACATAAACATGCATGTaa                                                                                             >  2:728805/1‑139 (MQ=255)
                                           aaCGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                         <  2:76942/139‑1 (MQ=255)
                                                 ccTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAACTACCATGtt                                                                                   >  1:982493/1‑139 (MQ=255)
                                                 ccTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGtt                                                                                   >  2:930477/1‑139 (MQ=255)
                                                      tAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCg                                                                              <  2:22596/139‑1 (MQ=255)
                                                         tGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACca                                                                                                      >  1:67321/1‑112 (MQ=255)
                                                         tGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACca                                                                                                      <  2:67321/112‑1 (MQ=255)
                                                            atatCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              >  2:550776/1‑84 (MQ=255)
                                                            atatCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              <  1:550776/86‑3 (MQ=255)
                                                               tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTatat                                                                                                                 >  1:327796/1‑95 (MQ=255)
                                                               tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTatat                                                                                                                 <  2:327796/95‑1 (MQ=255)
                                                               tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGa                                                                     >  2:630222/1‑139 (MQ=255)
                                                                  ggCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCg                                                                              >  1:1011809/1‑127 (MQ=255)
                                                                  ggCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCg                                                                              <  2:1011809/127‑1 (MQ=255)
                                                                  ggCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAAt                                                                  >  2:352813/1‑139 (MQ=255)
                                                                     aTGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATgcg                                                                                                              >  2:110861/1‑92 (MQ=255)
                                                                     aTGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATgcg                                                                                                              <  1:110861/92‑1 (MQ=255)
                                                                     aTGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAA‑CAAAAAAACAGTATTCTTATATGCGCATAAc                                                                                                        <  2:1053853/98‑1 (MQ=255)
                                                                     aTGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAA‑CAAAAAAACAGTATTCTTATATGCGCATAAc                                                                                                        >  1:1053853/1‑98 (MQ=255)
                                                                       gttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACg                                                             >  1:270457/1‑139 (MQ=255)
                                                                        ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                            <  1:47139/139‑1 (MQ=255)
                                                                        ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                            <  1:930477/139‑1 (MQ=255)
                                                                        ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                            <  2:637715/139‑1 (MQ=255)
                                                                        ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                            <  2:582941/139‑1 (MQ=255)
                                                                        ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                            <  2:953020/139‑1 (MQ=255)
                                                                           tACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATg                                                         >  1:645230/1‑139 (MQ=255)
                                                                                  aaTCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTatat                                                                                                                 >  1:759877/1‑76 (MQ=255)
                                                                                  aaTCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTatat                                                                                                                 <  2:759877/76‑1 (MQ=255)
                                                                                  aaTCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGa                                                  >  1:1046486/1‑139 (MQ=255)
                                                                                  aaTCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGa                                                  >  1:765407/1‑139 (MQ=255)
                                                                                     cGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAg                                               >  2:490167/1‑139 (MQ=255)
                                                                                       aTGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGAc                                             <  2:982493/139‑1 (MQ=255)
                                                                                                 ttGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTg                                   <  2:565500/139‑1 (MQ=255)
                                                                                                    aCAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              >  2:915244/1‑44 (MQ=255)
                                                                                                    aCAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              <  1:915244/46‑3 (MQ=255)
                                                                                                    aCAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATg                                >  1:626232/1‑139 (MQ=255)
                                                                                                            cTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGc                        <  2:891323/139‑1 (MQ=255)
                                                                                                            cTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGc                        <  1:490167/139‑1 (MQ=255)
                                                                                                              ggTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTa                                                                                 <  2:347731/80‑1 (MQ=255)
                                                                                                              ggTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTa                                                                                 >  1:347731/1‑80 (MQ=255)
                                                                                                              ggTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTg                      >  1:597744/1‑139 (MQ=255)
                                                                                                               gtgtTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              >  1:661825/1‑33 (MQ=255)
                                                                                                               gtgtTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGc                     >  1:260126/1‑139 (MQ=255)
                                                                                                               ctgtTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              <  2:661825/34‑3 (MQ=255)
                                                                                                                  ttttGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAAcc                                                                                                       >  2:803697/1‑54 (MQ=255)
                                                                                                                  ttttGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAAcc                                                                                                       <  1:803697/54‑1 (MQ=255)
                                                                                                                      gCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAgg              >  2:308315/1‑139 (MQ=255)
                                                                                                                           aaCATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACg         <  1:272577/139‑1 (MQ=255)
                                                                                                                             cATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGt       <  2:260126/139‑1 (MQ=255)
                                                                                                                                  gAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGtcttct  <  2:270457/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                    
CTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGTCTTCT  >  NC_000913/932196‑932465

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑