Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,543,231 A→C 5.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 5.5% 636.0 / 23.3 183M277R (ATG→AGG) ‡narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base A (92/81);  new base C (2/8);  total (94/89)
Fisher's exact test for biased strand distribution p-value = 5.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

Only 100 of 188 total aligned reads displayed.

GAACCACTTCCCAGACCCGCAGTCAGAAACAGCCCCATAAATACGGCGTAAAAGGCGATGAAATTACCGGAGCCTGTGCCCGGTAAGGTAAGGAACAGCAGGGCACTGAAAATCGCCATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTA  >  NC_000913/1542967‑1543495
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                         
gAACCACTTCCCAGACCCGCAGTCAGAAACAGCCCCATAAATACGGCGTAAAAGGCGATGAAATTACCGGAGCCTGTGCCCGGTAAGGTAAGGAACAGCAGGGCACTGAAAATCGCCATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATACGGGAACTGGGTTTTTGCCAGCATGGCAAAAcc                                                                                                                                                                                                                                                                >  2:271148/1‑275 (MQ=255)
                                                                                                agcagGGCACTGAAAATCGCCATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGAGAACTGGGTTTTTGCCAGCATGGCAAAACCCGAAGAAAAACCGATAAACGAACCGAAGGGGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGAAGGAGTTGTAAGACAGGGAGCTGGTAGGCact                                                                                                                                                                >  1:33124/1‑273 (MQ=255)
                                                                                                 gcagGGCACTGAAAATCGCCATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGGTGAAAGACCGGGGCCTGGTCGGCCAtt                                                                                                                                                               >  4:42399/1‑275 (MQ=255)
                                                                                                 gcagGGCACTGAAAATCGCCATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCCCAGCCGCCAGAGATGCCAGCGGTGTAAGACAGGGAGCGGGTCGGGAAtt                                                                                                                                                               >  4:170842/1‑275 (MQ=255)
                                                                                                 gcagGGCACTGAAAATCGCCATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCCCAGCAGCCAGAGCTGCAGGCGTTGGAAGAACGGGAGCTGGGCGCCAAt                                                                                                                                                                >  4:56043/1‑274 (MQ=255)
                                                                                                                 cGCCATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTgg                                                                                                                                               <  2:42399/275‑1 (MQ=255)
                                                                                                                     aTAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTTGCCAGCa                                                                                                                                                                                                                                                                           >  8:571790/1‑149 (MQ=255)
                                                                                                                     aTAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCAt                                                                                                                                                                                                                                                                          >  7:159067/1‑149 (MQ=255)
                                                                                                                     aTAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCAt                                                                                                                                                                                                                                                                          >  6:706042/1‑149 (MQ=255)
                                                                                                                     aTAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCAt                                                                                                                                                                                                                                                                          >  7:1318487/1‑149 (MQ=255)
                                                                                                                      tAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATg                                                                                                                                                                                                                                                                         <  7:1431647/149‑1 (MQ=255)
                                                                                                                      tAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATg                                                                                                                                                                                                                                                                         >  5:659338/1‑149 (MQ=255)
                                                                                                                      tAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATg                                                                                                                                                                                                                                                                         <  6:1441762/149‑1 (MQ=255)
                                                                                                                         aaaTAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCa                                                                                                                                                                                                                                                                     >  5:618973/1‑149 (MQ=255)
                                                                                                                            tAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCaaaa                                                                                                                                                                                                                                                                  >  6:1192196/1‑149 (MQ=255)
                                                                                                                            tAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTGCCATGCTGGCaaaaa                                                                                                                                                                                                                                                                 <  6:134372/149‑2 (MQ=255)
                                                                                                                             aaaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGCTGGCaaaaa                                                                                                                                                                                                                                                                 <  5:238393/149‑2 (MQ=255)
                                                                                                                             aaaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAACAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAAc                                                                                                                                                                                                                                                                 <  5:1268252/149‑1 (MQ=255)
                                                                                                                                gTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTGCCATGCTGGCAAAAaccca                                                                                                                                                                                                                                                             <  8:1098635/149‑3 (MQ=255)
                                                                                                                                 ttGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCGCCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGc                                                                                                                                                                                                                                                             >  6:962609/1‑149 (MQ=255)
                                                                                                                                 ttGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGc                                                                                                                                                                                                                                                             >  5:1367679/1‑149 (MQ=255)
                                                                                                                                 ttGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGc                                                                                                                                                                                                                                                             <  5:1097886/149‑1 (MQ=255)
                                                                                                                                 ttGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGc                                                                                                                                                                                                                                                             >  7:307952/1‑149 (MQ=255)
                                                                                                                                 ttGATCAACGTCACCCGCACGCCACCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGc                                                                                                                                                                                                                                                             <  7:1314691/149‑1 (MQ=255)
                                                                                                                                  tGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATccc                                                                                                                               >  2:315276/1‑274 (MQ=255)
                                                                                                                                     tCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTGCCATGCTGGCAAAAaccca                                                                                                                                                                                                                                                             <  5:1331260/144‑3 (MQ=255)
                                                                                                                                       aaCGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGaaaa                                                                                                                                                                                                                                                       >  5:486329/1‑149 (MQ=255)
                                                                                                                                        aCGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAAccc                                                                                                                                                                                                                                                               <  5:1363504/140‑1 (MQ=255)
                                                                                                                                          gTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGCCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCTGAAAAAcc                                                                                                                                                                                                                                                    >  6:1208479/1‑149 (MQ=255)
                                                                                                                                            cACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGa                                                                                                                                                                                                                                                  <  7:897587/149‑1 (MQ=255)
                                                                                                                                               ccGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATaa                                                                                                                                                                                                                                               >  6:747261/1‑149 (MQ=255)
                                                                                                                                                cGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGa                                                                                                                                                                                                                                                  <  5:214885/145‑1 (MQ=255)
                                                                                                                                                cGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGa                                                                                                                                                                                                                                                  >  6:214885/1‑145 (MQ=255)
                                                                                                                                                    cgccgcCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGaa                                                                                                                                                                                                                                          <  8:610515/149‑1 (MQ=255)
                                                                                                                                                     gccgccGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAAc                                                                                                                                                                                                                                         <  7:193140/149‑1 (MQ=255)
                                                                                                                                                     gccgccGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAAc                                                                                                                                                                                                                                         <  8:1024176/149‑1 (MQ=255)
                                                                                                                                                           gAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAgg                                                                                                                                                                                                                                   <  5:1192196/149‑1 (MQ=255)
                                                                                                                                                           gAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAgg                                                                                                                                                                                                                                   <  6:807928/149‑1 (MQ=255)
                                                                                                                                                            aaCTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGt                                                                                                                                                                                                                                  <  7:151678/149‑1 (MQ=255)
                                                                                                                                                               tgaTCGGAAAAAGCACCCCCAACCGGCCGCGGGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAg                                                                                                 <  4:271148/273‑1 (MQ=255)
                                                                                                                                                                 aTCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCaa                                                                                                                                                                                                                             >  5:187805/1‑149 (MQ=255)
                                                                                                                                                                  tCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAg                                                                                                                                                                                                                            >  5:784256/1‑149 (MQ=255)
                                                                                                                                                                     gAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTa                                                                                                                                                                                                                         >  8:739490/1‑149 (MQ=255)
                                                                                                                                                                        aTAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAg                                                                                                                                                                                                                      <  8:1328553/149‑1 (MQ=255)
                                                                                                                                                                         tAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGc                                                                                                                                                                                                                     <  7:748098/149‑1 (MQ=255)
                                                                                                                                                                          aGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCa                                                                                                                                                                                                                    <  8:1174888/149‑1 (MQ=255)
                                                                                                                                                                           gCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGAt                                                                                                                                                                                                                                                 <  5:360499/119‑1 (MQ=255)
                                                                                                                                                                           gCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAg                                                                                                                                                                                                                   >  8:809994/1‑149 (MQ=255)
                                                                                                                                                                            caccacCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAgg                                                                                                                                                                                                                  >  6:1123757/1‑149 (MQ=255)
                                                                                                                                                                              ccaccaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCt                                                                                                                                                                                                                >  6:685965/1‑149 (MQ=255)
                                                                                                                                                                               caccaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTc                                                                                                                                                                                                               >  8:772355/1‑149 (MQ=255)
                                                                                                                                                                                accaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTca                                                                                                                                                                                                              <  6:486329/149‑1 (MQ=255)
                                                                                                                                                                                   aaCCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTcagca                                                                                                                                                                                                           >  8:1415303/1‑149 (MQ=255)
                                                                                                                                                                                    accgaccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCATGc                                                                            >  4:237798/1‑275 (MQ=255)
                                                                                                                                                                                        accgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCCGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGAGGCGATCGt                                                                                                             >  1:57109/1‑238 (MQ=255)
                                                                                                                                                                                        accgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGata                                                                                                                                          <  2:16176/209‑1 (MQ=255)
                                                                                                                                                                                           gcgcgATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAgaga                                                                                                                                                                                                   <  6:437057/149‑1 (MQ=255)
                                                                                                                                                                                            cgcgATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGAt                                                                                                                                                                                                  <  7:410264/149‑1 (MQ=255)
                                                                                                                                                                                                    cACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGt                                                                                                                                                                                          <  6:265173/149‑1 (MQ=255)
                                                                                                                                                                                                     aCCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGtt                                                                                                                                                                                         <  7:670137/149‑1 (MQ=255)
                                                                                                                                                                                                     aCCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGtt                                                                                                                                                                                         >  6:1271434/1‑149 (MQ=255)
                                                                                                                                                                                                      ccGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTg                                                                                                                                                                                        >  5:554781/1‑149 (MQ=255)
                                                                                                                                                                                                      ccGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTg                                                                                                                                                                                        <  7:172363/149‑1 (MQ=255)
                                                                                                                                                                                                       cGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCagag                                                                                                                                                                                                    <  6:860495/136‑1 (MQ=255)
                                                                                                                                                                                                       cGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCagag                                                                                                                                                                                                    >  5:860495/1‑136 (MQ=255)
                                                                                                                                                                                                       cGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGt                                                                                                                                                                                       <  8:1093520/149‑1 (MQ=255)
                                                                                                                                                                                                        gATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATaaa                                                                                                                                                                                                                                              >  6:1157276/1‑93 (MQ=255)
                                                                                                                                                                                                        gATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTa                                                                                                                                                                                      >  8:339388/1‑149 (MQ=255)
                                                                                                                                                                                                         aTAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                                                                                                                                                     <  6:659338/149‑1 (MQ=255)
                                                                                                                                                                                                         aTAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                                                                                                                                                     <  6:1367679/149‑1 (MQ=255)
                                                                                                                                                                                                               gggCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAgg                                                                                                                                                                               <  6:63430/149‑1 (MQ=255)
                                                                                                                                                                                                                 gCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTg                                                                                                                                                                                                                                 >  5:205969/1‑97 (MQ=255)
                                                                                                                                                                                                                 gCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGa                                                                                                                                                                             <  7:1073518/149‑1 (MQ=255)
                                                                                                                                                                                                                   cAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGc                                                                                                                                                                           <  7:284776/149‑1 (MQ=255)
                                                                                                                                                                                                                    aaaGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCt                                                                                                                                                                          >  6:606284/1‑149 (MQ=255)
                                                                                                                                                                                                                     aaGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAgg                                                                                                                                                                               <  6:782507/143‑1 (MQ=255)
                                                                                                                                                                                                                         aCGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCg                                                                                                                                                                     >  8:1142828/1‑149 (MQ=255)
                                                                                                                                                                                                                           gCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                                                                                                                                                     <  8:1051734/131‑1 (MQ=255)
                                                                                                                                                                                                                           gCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                                                                                                                                                     >  7:1051734/1‑131 (MQ=255)
                                                                                                                                                                                                                               ggCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGAAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAAtt                                                                                                                                                               <  7:423540/149‑1 (MQ=255)
                                                                                                                                                                                                                                 cgcAGAATATTGACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAATCCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCaa                                                                                                                                                                 >  5:1207054/1‑145 (MQ=255)
                                                                                                                                                                                                                                 cgcAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGt                                                                                                                                                                                          >  4:65992/1‑120 (MQ=255)
                                                                                                                                                                                                                                 cgcAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGa                                                                                                                                                             <  6:936337/149‑1 (MQ=255)
                                                                                                                                                                                                                                  gcAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAg                                                                                                                                                            <  8:196140/149‑1 (MQ=255)
                                                                                                                                                                                                                                     gAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTa                                                                                                                                                                                      <  2:65992/120‑1 (MQ=255)
                                                                                                                                                                                                                                        tatTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGt                                                                                                                                                      <  5:1456042/149‑1 (MQ=255)
                                                                                                                                                                                                                                               aTCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTgg                                                                                                                                               >  8:13006/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  cGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGc                                                                                                                                              >  6:120103/1‑147 (MQ=255)
                                                                                                                                                                                                                                                  cGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATTCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCATGCGGCATTCGCCAGCGACCTCACCGGAACGTCGGCCTGCGGTACGCCATTGACGccagggaagg              >  2:304885/1‑267 (MQ=255)
                                                                                                                                                                                                                                                        cTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcat                                                                                                                                      <  7:718067/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         tGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcatt                                                                                                                                     >  5:1146885/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         tGGGTTTTGGCCAGAATGGCAAAACCCGCAGAAAAACCGAGAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTAAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATAc       <  2:170842/275‑1 (MQ=255)
                                                                                                                                                                                                                                                           ggTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGC‑‑GCGTTGTAAGACAGGGAGc                                                                                                                                                                           <  8:433684/107‑1 (MQ=255)
                                                                                                                                                                                                                                                            gggttttGCCATGCTGGCAAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAg                                                                                                                                                                            >  8:151807/4‑109 (MQ=255)
                                                                                                                                                                                                                                                            gTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattcat                                                                                                                                  >  5:147105/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              tttttGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATc                                                                                                                                 <  7:1142828/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               ttttGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATcc                                                                                                                                <  8:785578/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               gttggCCAGCAAGGCAAAACCCGAAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTa  <  2:58896/271‑1 (MQ=255)
                                                                                                                                                                                                                                                                tttGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATccc                                                                                                                               <  5:517924/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aTGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATgc                                                                                                                      <  5:138752/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                         
GAACCACTTCCCAGACCCGCAGTCAGAAACAGCCCCATAAATACGGCGTAAAAGGCGATGAAATTACCGGAGCCTGTGCCCGGTAAGGTAAGGAACAGCAGGGCACTGAAAATCGCCATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTA  >  NC_000913/1542967‑1543495

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