Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,281,959 A→T 5.8% Q699L (CAA→CTA)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,281,9590AT 5.8% 203.1 / 3.4 70Q699L (CAA→CTA) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (29/36);  new base T (2/2);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

GGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTACCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGC  >  NC_000913/1281813‑1282107
                                                                                                                                                  |                                                                                                                                                    
ggAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGaacaac                                                                                                                                                    <  1:1110906/149‑1 (MQ=255)
   aCGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGGCTATCGCTGGCGAAGCGGAAAAAGCCGAGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCAACCGTATTGCGGCAGAACAACCAg                                                                                                                                                 >  2:855010/1‑149 (MQ=255)
   aCGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAg                                                                                                                                                 >  2:511451/1‑149 (MQ=255)
   aCGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTGTTGCGGCAGAACACCCAg                                                                                                                                                 >  1:895603/1‑149 (MQ=255)
       aTGGGCTGGCTGCCGTCTGCACCCCAGTTAGGCACTAACCCGCTGACACTCGCTGTCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGaaaa                                                                                                                                             <  2:414649/149‑1 (MQ=255)
         gggCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGAAGAACAACCAGAACACg                                                                                                                                           <  2:517929/149‑1 (MQ=255)
         gggCTGGCTGCCGGCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGTCTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACg                                                                                                                                           <  2:828253/149‑1 (MQ=255)
          ggctggctGCCGTCTGCACCGCAGTTAGGCGCTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAACACgg                                                                                                                                          >  1:708807/1‑149 (MQ=255)
            ctggctgCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTa                                                                                                                                        >  1:25157/1‑149 (MQ=255)
             tggctgCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTaa                                                                                                                                       >  1:1343126/1‑149 (MQ=255)
                ccgCCGTCTGTACCGCAGTTCGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGACAGGGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTaaaaa                                                                                                                                    <  1:324825/147‑1 (MQ=255)
                  gCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAacc                                                                                                                                  <  1:747231/149‑1 (MQ=255)
                      tcggCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGAAAAATGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCAccc                                                                                                                              <  1:1265353/146‑1 (MQ=255)
                        tGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCcgc                                                                                                                            <  2:1267160/149‑1 (MQ=255)
                        tGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCACAAAGACGGATCCATCCGTTTTACGGCAGAACAACCAGAAAAGGGTAAAAAACACCcgg                                                                                                                            >  1:1189875/1‑148 (MQ=255)
                          cACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACTACCAGAAAACGGTAAAAACCACCCgcgc                                                                                                                          >  2:771607/1‑149 (MQ=255)
                          cACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACTACCAGAAAACGGTAAAAACCACCCgcgc                                                                                                                          >  2:771452/1‑149 (MQ=255)
                           aCCGCAGTTAGGCACTAACCCGCTGACTATCGGTCGCGAAGCGGGAAAAGCCGGGGTGGATCCGGTGGACTAAACGGTGAAATCCCTGAAAGAGGGGTCCATCCGTATTGCGGCAGAACAACCACAAAACGGCAAAGACCACCCGCGCa                                                                                                                         >  2:592283/1‑149 (MQ=255)
                           aCCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCCGAACACGGTAAAAACCACCCGCGCa                                                                                                                         >  2:592241/1‑149 (MQ=255)
                           aCCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGGTCCATCCGTTTTGCGGCAGAACAACCAGAACACGGTAAAAACCACCCGCGCa                                                                                                                         >  2:98280/1‑149 (MQ=255)
                            ccGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCaa                                                                                                                        >  2:692628/1‑149 (MQ=255)
                                       cACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCt                                                                                                             <  1:732399/149‑1 (MQ=255)
                                           acgccGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCg                                                                                                         <  2:911219/146‑1 (MQ=255)
                                           aaCCCGCTGACTAGCGCTGGCTGAGCGGAAATCGCCGGGATGAATCCGGTGGACTATTCGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCg                                                                                                         <  2:672426/149‑1 (MQ=255)
                                             cccGCTGACTTTCTCTGTCGAAGCGGAAATAGCCCGGATGAATCCGGTGGACTATACGGTGAAATCCATGAAAGAAGGAACCATCCGTTTTGCGGCAGAACAACAAGAAGCCGGTAAATGCCACCCGCGCAACCGGTTCATCTGGCGtt                                                                                                       <  2:899098/149‑1 (MQ=255)
                                              ccGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTc                                                                                                      >  1:1215183/1‑149 (MQ=255)
                                              ccGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTc                                                                                                      >  1:156459/1‑149 (MQ=255)
                                                  tGACTATCGCTGGCGAAGCGGAAAAAGCCGGAATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTCGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAAc                                                                                                  <  1:905114/149‑1 (MQ=255)
                                                     cTATCGCTGGCGAAGCGCAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAACCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTg                                                                                               <  2:459315/149‑1 (MQ=255)
                                                       aTCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGGGGGTTCCATCCGTTTTGCGGCAGCAAAACCAGGATACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCt                                                                                             >  2:342462/1‑149 (MQ=255)
                                                           cgggCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCATCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGt                                                                                         >  2:692439/3‑149 (MQ=255)
                                                           cTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAGAGGTAAAAACCACCCGCGCAACCTGTTCAGCGGGCGTTCTAACCTGccgggt                                                                                         >  2:268016/1‑144 (MQ=255)
                                                               cGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCATCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCCCGGttctt                                                                                     >  2:696833/1‑149 (MQ=255)
                                                                gAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCACCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGttcttc                                                                                    <  1:235248/149‑1 (MQ=255)
                                                                  aGCGGAAAAAGCCGGGATGAATCCGGGGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACACCCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCg                                                                                  <  2:320601/149‑1 (MQ=255)
                                                                     ggAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTa                                                                               >  1:423526/1‑149 (MQ=255)
                                                                      gAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTaa                                                                              >  1:526637/1‑149 (MQ=255)
                                                                      gAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCGGTTCATCTGGCGTGCTAACCTGCTGGGTTCTTCCGGGaa                                                                              >  1:279466/1‑149 (MQ=255)
                                                                       aaaaaGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTaaa                                                                             <  1:1075019/149‑1 (MQ=255)
                                                                       aaaaaGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGGTCTAACCTGCTCGGTTCTTCCGGGaaa                                                                             >  1:509110/1‑149 (MQ=255)
                                                                           aGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTc                                                                         <  1:119790/149‑1 (MQ=255)
                                                                            gCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCGGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCa                                                                        >  1:21380/1‑149 (MQ=255)
                                                                             ccGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCAt                                                                       <  2:714226/149‑1 (MQ=255)
                                                                              cGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATg                                                                      <  1:81136/149‑1 (MQ=255)
                                                                               gggATTAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAACCGGTAAAAACCACCCGCGAAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGa                                                                     <  1:465816/149‑1 (MQ=255)
                                                                                   tGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTACGGAAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTATAACCTGCTCGGTTATTACGGTAAAGGTCATggcttt                                                                 >  1:88389/1‑144 (MQ=255)
                                                                                    gAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGCTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTa                                                                <  1:652365/149‑1 (MQ=255)
                                                                                      aTCCGGTGGACTATACGGTGAAATCCCTGAAAGAGAGGTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTTTAACCTGCTCGGTTCTTCCGGTAAAGGACATGAGTTTATg                                                              >  1:240603/1‑149 (MQ=255)
                                                                                       tCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGc                                                             >  1:1275339/1‑149 (MQ=255)
                                                                                       accGGTGGACTATACGGTGAAATCCCTGAAACAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGc                                                             >  2:990555/2‑149 (MQ=255)
                                                                                        tcGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGCGTTTATGCt                                                            <  2:1172509/148‑1 (MQ=255)
                                                                                          ggTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCa                                                          <  1:1256163/149‑1 (MQ=255)
                                                                                            tGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAg                                                        >  1:947523/1‑149 (MQ=255)
                                                                                            tGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACAGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAg                                                        <  1:392544/149‑1 (MQ=255)
                                                                                                          aaaTCCCTGAAAGAGGGTTCCATCCGTTTGGCGGCAGAACAACCAGAAAACGGTAAAACCCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTACCTGCTGGGGAc                                          <  1:1137337/149‑1 (MQ=255)
                                                                                                            aTCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTACCTGCTGGGGACgg                                        <  2:1259637/149‑1 (MQ=255)
                                                                                                             tCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTACCTGCTGGGGACGGa                                       <  2:495986/149‑1 (MQ=255)
                                                                                                                 tGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACTACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTACCTGCTGGGGACGGAGCAc                                   <  1:771607/149‑1 (MQ=255)
                                                                                                                 tGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACTACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTACCTGCTGGGGACGGAGCAc                                   <  1:771452/149‑1 (MQ=255)
                                                                                                                     agggggTCCCCTCCGTTTTGCGGCAGAACAACCAGAAAACGGCAAAAACCGCCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTGGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTACCTGCTGGGGACGGAGCACGGTa                               <  2:911449/146‑1 (MQ=255)
                                                                                                                      gagGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTAGCTGCTGGGGACGGAGCACGGTAt                              <  1:886281/149‑1 (MQ=255)
                                                                                                                         ggTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTACCTGCTGGGGACGGAGCACGGTATCCa                           <  2:333129/149‑1 (MQ=255)
                                                                                                                          gTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTACCTGCTGGGGACGGAGCACGGTATCCAg                          <  2:1343126/149‑1 (MQ=255)
                                                                                                                           ttCCGTCCGTTTTGCGGCAGAACAACCAGAAAACGGTAACAACCACCCGCGGAAACTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTACCTGCTGGGGACGGAGCACGGTATCCAgg                         <  2:88389/149‑1 (MQ=255)
                                                                                                                           ttCCCTCCGTTTGGCGGCAGAACAACCAGGAAACGGTAAAAACCACCCGCGCAACCTGTTCAGCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTACCTGCTGGGGACGGAGCACGGTATCCAgg                         <  2:708807/149‑1 (MQ=255)
                                                                                                                           ttCCATCCGTTTTGCGGCAGAACCACCCGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTACCTGCTGGGGACGGAGCACGGTATCCAgg                         <  1:312344/149‑1 (MQ=255)
                                                                                                                           ttCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTTGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTGAAGTACCTGCTGGGGACGGAGCACGGTATCCAgg                         <  1:683105/149‑1 (MQ=255)
                                                                                                                               aTCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAGACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTAATGCTCCAGTACCTGCTGGGGACGGAGCACGGTATCCAGGGTaa                     >  1:40902/1‑149 (MQ=255)
                                                                                                                                tCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAGCCTGTTCATCCGGCGTTCTAACCTGCTCGGTTCTTCCTGTAAAGGTCATGAGTTTATGCTTAAGTACCTGCTGGGGACGGAGCACGATATCCAGGGTaaa                    <  2:826640/149‑1 (MQ=255)
                                                                                                                                       tGCGGCAGAACAACCAGAAAACGGTAAAACCCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTACCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTgg             <  1:665049/149‑1 (MQ=255)
                                                                                                                                         cGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTACCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTgggg           >  1:693164/1‑149 (MQ=255)
                                                                                                                                               aacaacCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCGAGTAGCTGCTGGGGACGGAGCACGGGATCCACGGGGAAGATCTGGGGccgcgg     >  2:475771/1‑144 (MQ=255)
                                                                                                                                                  aacCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTACCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGggc  >  2:114145/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTACCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGC  >  NC_000913/1281813‑1282107

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