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)
             ccGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTc                                                                                                                       >  1:514161/1‑139 (MQ=255)
             ccGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTc                                                                                                                       >  1:565500/1‑139 (MQ=255)
                 tttATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTtat                                                                                                                   <  1:8786/139‑1 (MQ=255)
                   tataAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTatat                                                                                                                 <  1:989581/139‑1 (MQ=255)
                           gATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATaa                                                                                                         >  1:76942/1‑139 (MQ=255)
                           gATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATaa                                                                                                         >  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)
                              tAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACca                                                                                                      <  2:514161/139‑1 (MQ=255)
                                aaTGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACcatg                                                                                                    <  2:255122/139‑1 (MQ=255)
                                aaTGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACcatg                                                                                                    <  1:113364/139‑1 (MQ=255)
                                    atgaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACcatgcatg                                                                                                >  2:113254/1‑139 (MQ=255)
                                      gaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                <  1:624399/106‑1 (MQ=255)
                                      gaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                >  2:624399/1‑106 (MQ=255)
                                      gaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTa                                                                                              >  1:22596/1‑139 (MQ=255)
                                       aaTAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGAAAATCCACTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTATTATATGCACATAAACATGCATGTaa                                                                                             >  2:728805/1‑139 (MQ=255)
                                           aaCGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                         <  2:76942/139‑1 (MQ=255)
                                                 ccTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAACTACCATGtt                                                                                   >  1:982493/1‑139 (MQ=255)
                                                 ccTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGtt                                                                                   >  2:930477/1‑139 (MQ=255)
                                                      tAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCg                                                                              <  2:22596/139‑1 (MQ=255)
                                                         tGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACca                                                                                                      >  1:67321/1‑112 (MQ=255)
                                                         tGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACca                                                                                                      <  2:67321/112‑1 (MQ=255)
                                                            atatCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              >  2:550776/1‑84 (MQ=255)
                                                            atatCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              <  1:550776/86‑3 (MQ=255)
                                                               tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTatat                                                                                                                 >  1:327796/1‑95 (MQ=255)
                                                               tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTatat                                                                                                                 <  2:327796/95‑1 (MQ=255)
                                                               tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGa                                                                     >  2:630222/1‑139 (MQ=255)
                                                                  ggCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCg                                                                              >  1:1011809/1‑127 (MQ=255)
                                                                  ggCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCg                                                                              <  2:1011809/127‑1 (MQ=255)
                                                                  ggCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAAt                                                                  >  2:352813/1‑139 (MQ=255)
                                                                     aTGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATgcg                                                                                                              >  2:110861/1‑92 (MQ=255)
                                                                     aTGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATgcg                                                                                                              <  1:110861/92‑1 (MQ=255)
                                                                     aTGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAACAAAAAAACAGTATTCTTATATGCGCATAAc                                                                                                        <  2:1053853/98‑1 (MQ=255)
                                                                     aTGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAACAAAAAAACAGTATTCTTATATGCGCATAAc                                                                                                        >  1:1053853/1‑98 (MQ=255)
                                                                       gttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACg                                                             >  1:270457/1‑139 (MQ=255)
                                                                        ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                            <  1:47139/139‑1 (MQ=255)
                                                                        ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                            <  1:930477/139‑1 (MQ=255)
                                                                        ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                            <  2:637715/139‑1 (MQ=255)
                                                                        ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                            <  2:582941/139‑1 (MQ=255)
                                                                        ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                            <  2:953020/139‑1 (MQ=255)
                                                                           tACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATg                                                         >  1:645230/1‑139 (MQ=255)
                                                                                  aaTCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTatat                                                                                                                 >  1:759877/1‑76 (MQ=255)
                                                                                  aaTCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTatat                                                                                                                 <  2:759877/76‑1 (MQ=255)
                                                                                  aaTCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGa                                                  >  1:1046486/1‑139 (MQ=255)
                                                                                  aaTCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGa                                                  >  1:765407/1‑139 (MQ=255)
                                                                                     cGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAg                                               >  2:490167/1‑139 (MQ=255)
                                                                                       aTGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGAc                                             <  2:982493/139‑1 (MQ=255)
                                                                                                 ttGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTg                                   <  2:565500/139‑1 (MQ=255)
                                                                                                    aCAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              >  2:915244/1‑44 (MQ=255)
                                                                                                    aCAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              <  1:915244/46‑3 (MQ=255)
                                                                                                    aCAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATg                                >  1:626232/1‑139 (MQ=255)
                                                                                                            cTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGc                        <  2:891323/139‑1 (MQ=255)
                                                                                                            cTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGc                        <  1:490167/139‑1 (MQ=255)
                                                                                                              ggTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTa                                                                                 <  2:347731/80‑1 (MQ=255)
                                                                                                              ggTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTa                                                                                 >  1:347731/1‑80 (MQ=255)
                                                                                                              ggTGTTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTg                      >  1:597744/1‑139 (MQ=255)
                                                                                                               gtgtTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              >  1:661825/1‑33 (MQ=255)
                                                                                                               gtgtTTTGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGc                     >  1:260126/1‑139 (MQ=255)
                                                                                                               ctgtTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              <  2:661825/34‑3 (MQ=255)
                                                                                                                  ttttGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAAcc                                                                                                       >  2:803697/1‑54 (MQ=255)
                                                                                                                  ttttGCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAAcc                                                                                                       <  1:803697/54‑1 (MQ=255)
                                                                                                                      gCGAAAACATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAgg              >  2:308315/1‑139 (MQ=255)
                                                                                                                           aaCATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACg         <  1:272577/139‑1 (MQ=255)
                                                                                                                             cATTCGAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGt       <  2:260126/139‑1 (MQ=255)
                                                                                                                                  gAGGAAGAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGtcttct  <  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: