Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,965 A→G 7.4% intergenic (+439/+101) narI → / ← rttR nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,9650AG 7.4% 606.6 / 13.9 229intergenic (+439/+101)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (80/120);  new base G (13/3);  total (100/129)
Fisher's exact test for biased strand distribution p-value = 2.64e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 229 total aligned reads displayed.

TCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG  >  NC_000913/1286817‑1287112
                                                                                                                                                    |                                                                                                                                                    
tcATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGAc                                                                                                                                                     >  1:5558439/1‑150 (MQ=255)
   tCCTTCCCCGCGTGGGCAGAATCTTTGATTGCGGCTTAGTTTGAGTATTCCGGGGCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCt                                                                                                                                                  <  2:5685704/150‑1 (MQ=14)
     cTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCg                                                                                                                                                >  1:5712132/1‑150 (MQ=255)
      ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                                               >  2:5936591/1‑150 (MQ=255)
      gtCCCCGCGTGGGCAGACGATTTGATTGGGGAGTCGTTTGAGTATTCGTGTGCTTTAGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                                               <  2:1687177/149‑1 (MQ=14)
         cccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGA‑gac                                                                                                                                           >  2:4676557/1‑148 (MQ=255)
          ccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGAGCTGAGCGCTTCGCCCTTCGGGTCGATGCCTGCAGCGACGCTCTCTCGCTGTCGCTCGA‑GACg                                                                                                                                          >  2:5532466/1‑150 (MQ=14)
            gCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGaa                                                                                                                                        <  2:3366295/150‑1 (MQ=255)
            gCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACAAAGAGCTGCGCGCTACGCCCGTCGGGTCCTTGCCGTCGGCAACCCACGCTAGCAGGCGCGCGA‑GTCGaa                                                                                                                                        >  1:6090516/1‑150 (MQ=11)
                 ggCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGA‑GTCGAACCttg                                                                                                                                   >  2:6105514/1‑149 (MQ=18)
                 ggCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTa                                                                                                                                   >  2:3644977/1‑150 (MQ=255)
                 ggCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTa                                                                                                                                   <  2:3438948/150‑1 (MQ=255)
                 ggCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTa                                                                                                                                   <  2:3490442/150‑1 (MQ=255)
                 ggCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTa                                                                                                                                   >  2:3521806/1‑150 (MQ=255)
                 ggCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCGCACACGCGGACGCTCGA‑GCCCaaacata                                                                                                                                   >  1:2515625/1‑145 (MQ=14)
                  gCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGA‑gt                                                                                                                                            >  1:3803453/1‑139 (MQ=37)
                   cAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCAGCGGCAAAGCAATCACGCTGGCGCTCGA‑GACGAAccgtggt                                                                                                                                 >  1:5830419/1‑145 (MQ=11)
                          tttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCATCGCTGCGCGCATCGACCTTCGGGACGTTGCCTGCGGCAAAGCACTCACGCAGACGCTCCA‑GCCGAAACGTAGTCGAAGCt                                                                                                                          >  1:5011654/1‑150 (MQ=14)
                                 gCGGATTCGTTTGAGAATTCCGGGGCTTTAGACAGTGATGGTGGGGGGGGAAGGATTACTCCGCGCTGCGCGCTTCGCCCGTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCGtctccgc                                                                                                                   <  2:1825669/150‑4 (MQ=14)
                                  cGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGA‑GCCGAACCATAGTCGAAGCGTCTCATcc                                                                                                                  >  2:6294909/1‑150 (MQ=255)
                                    gATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGACCCTCGA‑GCCGAACCATAGTCGAAGCTTAACATccat                                                                                                                >  2:2614854/1‑148 (MQ=14)
                                      ttCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTcc                                                                                                              <  2:2277422/150‑1 (MQ=255)
                                        cGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTcccc                                                                                                            <  2:4864413/150‑1 (MQ=255)
                                          tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGA‑GTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGc                                                                                                          >  2:698990/1‑150 (MQ=2)
                                          tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGc                                                                                                          >  2:4655786/1‑150 (MQ=255)
                                            tGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTAGCTGGCGCTAGA‑GCCGAACCATAGTCGAAGCTTAACATACTTAACCGcca                                                                                                        >  1:1875585/1‑148 (MQ=14)
                                              agaATTCGTGTGCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                                      <  1:4333706/150‑1 (MQ=255)
                                               gtaTTCCTGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                     <  2:6648694/148‑1 (MQ=255)
                                               gatttCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                     <  1:4526398/147‑1 (MQ=255)
                                               gaATTCCGGGGCTTTTGCAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                     <  1:2899613/150‑1 (MQ=255)
                                               gaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                     >  2:4594107/1‑150 (MQ=255)
                                                aaTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGAGTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGc                                                                                                    <  1:330654/150‑1 (MQ=255)
                                                  ttCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAg                                                                                                  <  1:2597654/150‑1 (MQ=255)
                                                    gcGGGGCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaa                                                                                                <  2:380173/149‑1 (MQ=255)
                                                        ggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                                            <  2:2805045/150‑1 (MQ=38)
                                                        ggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GACGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                                            >  2:5594679/1‑150 (MQ=25)
                                                         tcTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                           <  2:3275472/149‑1 (MQ=37)
                                                         tcTTATGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                           <  1:1116077/149‑1 (MQ=255)
                                                         gCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACCCACACACCCAGACGCTCCA‑GCCGAACCATAGTCGAATCTTATCAAAAATACCCCCAAGGGCAGAaaatat                                                                                           >  1:1959188/1‑148 (MQ=2)
                                                         gCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                           <  2:4400802/150‑1 (MQ=38)
                                                            tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCATGCGGCAACGCACTCTCGCTGACCCTCGATGCCCAA‑CATAGTCGAATCAACTCAAACTTACCCCCATGGGCAGAATAtttttt                                                                                        >  1:992775/1‑147 (MQ=11)
                                                              tGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTg                                                                                      <  2:4227238/150‑1 (MQ=39)
                                                               gAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGc                                                                                     >  2:566940/1‑150 (MQ=39)
                                                               gAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCACGCTGACGCTCGA‑GCCGAACCGTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGc                                                                                     >  2:1424292/1‑150 (MQ=255)
                                                                aaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCg                                                                                    <  2:5558439/150‑1 (MQ=39)
                                                                 aaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCgg                                                                                   <  1:2305043/150‑1 (MQ=39)
                                                                 aaGTGATGGTGGTGGGGGAAGGATTACTAAGCGCTGCGCGCTTCGCCCTTCTGGTCGTTGCGTGAGGCAAAGCTCTCTCGGTGACGCTGGA‑GGCGAAACTTTGTCGTAGGTTCTCATGCTTCCGCGCATGGGCGGTATATTTGTTTGCgg                                                                                   >  2:2804778/1‑150 (MQ=11)
                                                                        ggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCt                                                                            <  1:5074410/150‑1 (MQ=40)
                                                                        ggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCt                                                                            <  1:5356802/150‑1 (MQ=40)
                                                                         gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCGCACTAGCTGACGCGCGA‑GCCGAACCAAAGTCGCAGCTTCTCAACCATACCCCCAGGGGCAGAATGATTTTTTTCGGATTCGCtt                                                                           >  1:4826594/1‑150 (MQ=14)
                                                                           ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                         <  1:1306421/150‑1 (MQ=255)
                                                                           ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                         <  1:1752045/150‑1 (MQ=255)
                                                                                  gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                                                  <  2:6386279/150‑1 (MQ=255)
                                                                                  gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                                                  <  1:1790116/150‑1 (MQ=255)
                                                                                  gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                                                  <  1:1182545/150‑1 (MQ=255)
                                                                                   aaGGATTACTCCGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                 <  1:1851249/150‑1 (MQ=255)
                                                                                    aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                                <  1:870286/150‑1 (MQ=255)
                                                                                      gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTGGGGAGTTCAGGGa                                                              >  2:5018172/1‑150 (MQ=255)
                                                                                       aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGAc                                                             <  1:2984631/150‑1 (MQ=255)
                                                                                           cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACGCTCTCTCGCTGGCGCTCGA‑GCCGAACCTTAGTCGAAGCGTATCATCCTTACCCGCCTGGGCCGAATATTTGTTTGCGGATTCGCTGGAGTGTGGAGGGACtgtt                                                         >  2:59498/1‑150 (MQ=11)
                                                                                           cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                         <  1:2725732/150‑1 (MQ=255)
                                                                                            tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTg                                                        <  1:5215413/150‑1 (MQ=255)
                                                                                              aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCGTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGag                                                      >  2:3148464/1‑149 (MQ=255)
                                                                                               gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                     <  1:6518188/150‑1 (MQ=255)
                                                                                                 gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑g                                                                                                                                             >  2:2458733/1‑60 (MQ=35)
                                                                                                  cTGCGCGCTTCGCACTGCGGGGCGTTGACTGCGGCAACGCTCTCTCGCTGTCGCTCGA‑GCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTg                                                  <  2:5498128/150‑1 (MQ=14)
                                                                                                   tGCGCGCTTGGGCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATTTTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGa                                                 <  1:6165448/150‑1 (MQ=255)
                                                                                                   tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGa                                                 <  1:2640867/150‑1 (MQ=255)
                                                                                                     cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtg                                               <  2:1765946/150‑1 (MQ=255)
                                                                                                     cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtg                                               <  2:832870/150‑1 (MQ=255)
                                                                                                     cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCAGGGGCAGAATATTTGATTGCGGGTTAGGTTGTGAGTTGAGGGACGTTTGAAAGGGAtg                                               >  2:2437510/1‑150 (MQ=14)
                                                                                                     cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCGTCTCATCCATCCCCGCAGGGGCAGAATATTTGTTTGCGGATTCGCGTGTGAGGTCCGGGGCTTTTGGAAGGGAtg                                               >  2:3110718/1‑150 (MQ=14)
                                                                                                       cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggt                                             <  2:5107963/150‑1 (MQ=255)
                                                                                                            cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                                        <  1:3295936/150‑1 (MQ=255)
                                                                                                            cGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCTCGCTGACGCTCGA‑GGCGAACCTTAGTCGCAGCTTCTCATCCTTCCCCGCAAGGGCAGAATAATTGATTGCGGATTCGGTTGGGAGTTCAGGGACTTTTGAAAGTGAtggaggtgg                                        >  2:772171/1‑150 (MQ=14)
                                                                                                            cGCCATTCGGGTCGTTGCCTGCGGCATCGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                                        <  1:6625850/150‑1 (MQ=255)
                                                                                                             gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                                       <  2:2922979/150‑1 (MQ=255)
                                                                                                              cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                                      <  1:1048176/150‑1 (MQ=255)
                                                                                                                 gagagagcGTTGCCGCAGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                   >  2:5895194/8‑150 (MQ=17)
                                                                                                                  tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCGGACGCTCGA‑GCCGAACCTTAGGCGAAGCTTCTCAGCCGTCCCCGCATGGGCAGAATCTTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAg                                  <  2:4448058/150‑1 (MQ=17)
                                                                                                                  tCGGGTCGGTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGA‑GCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAg                                  <  2:4244500/150‑1 (MQ=255)
                                                                                                                    agagcGTTGCCGCAGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCATCCCCGCAAGGGCAGAAAATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                >  1:5443952/5‑150 (MQ=17)
                                                                                                                     ggTCGTTGCCTGCGGCAAAGCTCTCTAGCTGACGCTAGA‑GCCGAACCTTAGTCGAAGCTTCTAAACCTTACCCGCATGGGCCGAAGATTTGTTTGCGGATTCGGGTGTGGGTTCAGGGACTGTTGTAAGTGGTGGTGGTGGGGGCAggct                               >  2:1391025/1‑148 (MQ=255)
                                                                                                                      gTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGGGAGTTCAGGGACTTTTGATAGTGATGGTGGTGGGGGAAGGAtt                              >  2:4981892/1‑150 (MQ=255)
                                                                                                                       tCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa                             <  1:4271403/150‑1 (MQ=255)
                                                                                                                       tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATGTGATTGCGGATTCGCTTGCGAGTTCAGGGACTTTTGAAAGTGATGGTGGGGGGGGACGGATTa                             >  2:4009415/1‑150 (MQ=255)
                                                                                                                        cggtGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAc                            <  1:4596925/147‑1 (MQ=255)
                                                                                                                          ttGCCTGCGGCAACGCTCTCACGCAGACGCTCGA‑GCCGAACCATAGTCGAAGCTTCTCATACATCCCCCCCTGGGCAGAAAATTTGATTGCGGAGTCGCTTGAGAGTTGAGGGACTTTTGAAAGTGGTGGTGGTGGGGGGAGGATTAcac                          >  1:6185894/1‑148 (MQ=17)
                                                                                                                          gtGCCTGCGGCATCGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                          <  1:494341/149‑1 (MQ=255)
                                                                                                                           tGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCa                         <  1:3521806/150‑1 (MQ=255)
                                                                                                                              cTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTTGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCacag                      >  2:2314941/1‑147 (MQ=255)
                                                                                                                                    cAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgc                <  1:5574181/150‑1 (MQ=255)
                                                                                                                                    cAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgc                <  1:5941935/150‑1 (MQ=255)
                                                                                                                                        gCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTTGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCCTCAGCTCTGCGCGCtt            >  2:5866791/1‑150 (MQ=255)
                                                                                                                                        gCTCTCTCGCTGACGCTCGA‑GCCGAACCATAGTCGAAGCATCACAACCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACACAGCGCTGCGAGCtt            >  1:4416109/1‑150 (MQ=17)
                                                                                                                                          tctctcGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg          <  2:5713949/150‑1 (MQ=255)
                                                                                                                                          tctctcGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCAAACTTACCCGCAAGCGCAGAAAATTGGATAGCGGATTCGCATGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg          >  1:2454754/1‑150 (MQ=17)
                                                                                                                                            tctcGCTGCCGCTCGA‑GCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGcc        <  2:6441858/150‑1 (MQ=255)
                                                                                                                                               cGCTGTCGCTCGA‑GCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCtt     <  2:6319850/150‑1 (MQ=255)
                                                                                                                                                gCTGGCGCTCGA‑GCCGAACCTTCGTCGAAGCTTCTCATCCGTCCCCGCTTGGGCAGAATATGTGAGTGCGGATTCGCTTGAGAGTTCCGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc    <  1:6170465/150‑1 (MQ=11)
                                                                                                                                                    |                                                                                                                                                    
TCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGA‑GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG  >  NC_000913/1286817‑1287112

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑