Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,543,231 A→C 6.5% M277R (ATG→AGG) ‡ narU ← nitrate/nitrite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,543,2310AC 6.5% 541.0 / 24.1 154M277R (ATG→AGG) ‡narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base A (72/72);  new base C (2/8);  total (74/80)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

Only 100 of 159 total aligned reads displayed.

CATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGC  >  NC_000913/1543083‑1543379
                                                                                                                                                     |                                                                                                                                                     
cATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCa                                                                                                                                                       >  1:15320/1‑149 (MQ=255)
 aTAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTTGCCAGCa                                                                                                                                                       >  4:571790/1‑149 (MQ=255)
 aTAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCAt                                                                                                                                                      >  2:706042/1‑149 (MQ=255)
 aTAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCAt                                                                                                                                                      >  3:1318487/1‑149 (MQ=255)
  tAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATg                                                                                                                                                     <  3:1431647/149‑1 (MQ=255)
  tAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATg                                                                                                                                                     >  1:659338/1‑149 (MQ=255)
     aaaTAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCa                                                                                                                                                 >  1:618973/1‑149 (MQ=255)
      aaTAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGCTGGCaa                                                                                                                                                >  3:784930/1‑149 (MQ=255)
        tAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCaaaa                                                                                                                                              <  2:618973/149‑1 (MQ=255)
        tAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTGCCATGCTGGCaaaaa                                                                                                                                             <  2:134372/149‑2 (MQ=255)
         aaaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGCTGGCaaaaa                                                                                                                                             <  4:784930/149‑2 (MQ=255)
         aaaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGCTGGCaaaaa                                                                                                                                             <  1:238393/149‑2 (MQ=255)
            gTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTGCCATGCTGGCAAAAaccca                                                                                                                                         <  4:1098635/149‑3 (MQ=255)
             ttGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCGCCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGc                                                                                                                                         >  2:962609/1‑149 (MQ=255)
             ttGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGc                                                                                                                                         <  1:1097886/149‑1 (MQ=255)
             ttGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGc                                                                                                                                         >  1:1367679/1‑149 (MQ=255)
             ttGATCAACGTCACCCGCACGCCACCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGc                                                                                                                                         <  3:1314691/149‑1 (MQ=255)
                  cAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGaaa                                                                                                                                    <  2:1300603/149‑1 (MQ=255)
                   aaCGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGaaaa                                                                                                                                   >  1:486329/1‑149 (MQ=255)
                    aCGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAAccc                                                                                                                                           <  1:1363504/140‑1 (MQ=255)
                      gTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAAcc                                                                                                                                >  4:748098/1‑149 (MQ=255)
                      gTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAAcc                                                                                                                                >  2:923421/1‑149 (MQ=255)
                         aCCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGAt                                                                                                                             >  3:167628/1‑149 (MQ=255)
                            cGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGa                                                                                                                              <  1:214885/145‑1 (MQ=255)
                                cgccgcCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGaa                                                                                                                      <  4:610515/149‑1 (MQ=255)
                                 gccgccGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAAc                                                                                                                     <  3:193140/149‑1 (MQ=255)
                                       gAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAgg                                                                                                               <  1:1192196/149‑1 (MQ=255)
                                       gAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAgg                                                                                                               <  2:807928/149‑1 (MQ=255)
                                        aaCTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGt                                                                                                              <  3:151678/149‑1 (MQ=255)
                                         aCTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACAAACCGAAGGTg                                                                                                             <  3:962052/149‑1 (MQ=255)
                                           ttATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGc                                                                                                           >  2:1173096/1‑149 (MQ=255)
                                              tCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAg                                                                                                        >  1:853960/1‑149 (MQ=255)
                                                 gAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTa                                                                                                     >  4:739490/1‑149 (MQ=255)
                                                 gAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTa                                                                                                     >  3:371218/1‑149 (MQ=255)
                                                 gAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTa                                                                                                     >  4:1165860/1‑149 (MQ=255)
                                                    aTAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAg                                                                                                  <  4:1328553/149‑1 (MQ=255)
                                                     tAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGc                                                                                                 <  3:748098/149‑1 (MQ=255)
                                                      aGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCa                                                                                                >  1:1409117/1‑149 (MQ=255)
                                                       gCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAg                                                                                               >  4:809994/1‑149 (MQ=255)
                                                        caccacCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAgg                                                                                              >  2:1123757/1‑149 (MQ=255)
                                                         accaccAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGc                                                                                             >  1:1222846/1‑149 (MQ=255)
                                                          ccaccaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCt                                                                                            >  2:685965/1‑149 (MQ=255)
                                                            accaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTca                                                                                          <  2:486329/149‑1 (MQ=255)
                                                              caACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTcagc                                                                                        >  3:711920/1‑149 (MQ=255)
                                                               aaCCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTcagca                                                                                       >  4:1415303/1‑149 (MQ=255)
                                                                  cgaccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGcc                                                                                    <  4:1307425/149‑1 (MQ=255)
                                                                    accgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCag                                                                                  <  2:763372/149‑1 (MQ=255)
                                                                     ccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCaga                                                                                 <  2:721663/149‑1 (MQ=255)
                                                                       gcgcgATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAgaga                                                                               <  2:437057/149‑1 (MQ=255)
                                                                        cgcgATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGAt                                                                              <  3:410264/149‑1 (MQ=255)
                                                                                cACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGt                                                                      <  2:265173/149‑1 (MQ=255)
                                                                                 aCCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGtt                                                                     >  1:202359/1‑149 (MQ=255)
                                                                                  ccGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTg                                                                    >  2:656768/1‑149 (MQ=255)
                                                                                  ccGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTg                                                                    >  1:554781/1‑149 (MQ=255)
                                                                                  ccGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTg                                                                    <  3:172363/149‑1 (MQ=255)
                                                                                   cGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCagag                                                                                >  1:860495/1‑136 (MQ=255)
                                                                                   cGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGt                                                                   <  4:1093520/149‑1 (MQ=255)
                                                                                    gATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATaaa                                                                                                                          <  1:1157276/93‑1 (MQ=255)
                                                                                    gATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATaaa                                                                                                                          >  2:1157276/1‑93 (MQ=255)
                                                                                    gATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTa                                                                  >  2:810483/1‑149 (MQ=255)
                                                                                    gATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTa                                                                  >  4:339388/1‑149 (MQ=255)
                                                                                     aTAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                                 <  2:1367679/149‑1 (MQ=255)
                                                                                     aTAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                                 <  2:222984/149‑1 (MQ=255)
                                                                                     aTAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGAATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                                 <  4:1181305/149‑1 (MQ=255)
                                                                                           gggCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAgg                                                           <  4:159067/149‑1 (MQ=255)
                                                                                             gCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTg                                                                                                             <  2:205969/97‑1 (MQ=255)
                                                                                             gCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTg                                                                                                             >  1:205969/1‑97 (MQ=255)
                                                                                              ccAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAg                                                        >  1:1363102/1‑149 (MQ=255)
                                                                                              ccAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAg                                                        <  1:270113/149‑1 (MQ=255)
                                                                                               cAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGc                                                       <  3:284776/149‑1 (MQ=255)
                                                                                                aaaGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCt                                                      >  2:606284/1‑149 (MQ=255)
                                                                                                 aaGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAgg                                                           >  1:782507/1‑143 (MQ=255)
                                                                                                  aGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTgg                                                    >  1:1381766/1‑149 (MQ=255)
                                                                                                    aaCGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTc                                                  <  1:1271434/149‑1 (MQ=255)
                                                                                                     aCGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAGCGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCg                                                 >  1:1066545/1‑149 (MQ=255)
                                                                                                     aCGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCg                                                 >  4:1142828/1‑149 (MQ=255)
                                                                                                       gCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                                 <  4:1051734/131‑1 (MQ=255)
                                                                                                       gCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                                 >  3:1051734/1‑131 (MQ=255)
                                                                                                         cAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCaa                                             >  1:479114/1‑149 (MQ=255)
                                                                                                             cgcAGAATATTGACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAATCCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCaa                                             >  1:1207054/1‑145 (MQ=255)
                                                                                                             cgcAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGa                                         <  2:936337/149‑1 (MQ=255)
                                                                                                             cgcAAAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCaa                                             <  2:1207054/145‑1 (MQ=255)
                                                                                                              gcAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAg                                        <  2:1381766/149‑1 (MQ=255)
                                                                                                              gcAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAg                                        >  4:1047127/1‑149 (MQ=255)
                                                                                                                aGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGgc                                      <  4:367544/149‑1 (MQ=255)
                                                                                                                 gAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGgcg                                     >  2:66953/1‑149 (MQ=255)
                                                                                                                 gAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGgc                                      <  3:1141462/148‑1 (MQ=255)
                                                                                                                  aaTATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGgcgc                                    >  3:1055693/1‑149 (MQ=255)
                                                                                                                   atatTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGCGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGcgcg                                   >  3:30769/1‑149 (MQ=255)
                                                                                                                              cGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGc                          >  2:120103/1‑147 (MQ=255)
                                                                                                                                    cTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcat                  <  3:718067/149‑1 (MQ=255)
                                                                                                                                       ggTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGC‑‑GCGTTGTAAGACAGGGAGc                                                       >  3:433684/1‑107 (MQ=255)
                                                                                                                                        gggttttGCCATGCTGGCAAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAg                                                        <  3:151807/106‑1 (MQ=255)
                                                                                                                                        gTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattcat              >  1:147105/1‑149 (MQ=255)
                                                                                                                                          tttttGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATc             <  3:1142828/149‑1 (MQ=255)
                                                                                                                                          tttttGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATc             <  4:145491/149‑1 (MQ=255)
                                                                                                                                          tttttGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATc             <  2:124416/149‑1 (MQ=255)
                                                                                                                                           ttttGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATcc            <  4:785578/149‑1 (MQ=255)
                                                                                                                                            tttGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATccc           <  1:517924/149‑1 (MQ=255)
                                                                                                                                                     aTGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATgc  <  4:451878/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGC  >  NC_000913/1543083‑1543379

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