Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,894,532 T→G S113A (TCA→GCA)  ygcP → putative antiterminator regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,894,5320TG96.9% 367.4 / ‑2.2 96S113A (TCA→GCA) ygcPputative antiterminator regulatory protein
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (47/46);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TGGAAGGCGCGTCAAATAAAGAAGTCGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGA  >  NC_000913/2894386‑2894677
                                                                                                                                                  |                                                                                                                                                 
tGGAAGGCCAGTCAAATAAAGAAGTCCTTATTCAGTTTTTGAAACTGGGGACCGAGGCGGACGCCATTATCAGTACCAAAGCCTCAATGTTGCAAACGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCCGTTTATTGTTGATGCAAt                                                                                                                                               >  2:480790/1‑151 (MQ=255)
 ggAAGGCGCGTCAAATAAAGAAGTCGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAAtt                                                                                                                                              <  1:680464/151‑1 (MQ=255)
   aaGGCGCGTCAAATAAAGAAGTCGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTc                                                                                                                                            >  2:1093599/1‑151 (MQ=255)
   aaGGCGCGTCAAATAAAGAAGTCGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTc                                                                                                                                            <  2:793613/151‑1 (MQ=255)
          gTCAAATAAAGAAGTCGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTcac                                                                                                                                     >  2:277942/1‑151 (MQ=255)
             aaaTAAAGAAGTCGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCacaac                                                                                                                                  >  1:24794/1‑151 (MQ=255)
               aTAAAGAAGTCGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACAt                                                                                                                                >  2:280275/1‑151 (MQ=255)
               aTAAAGAAGTCGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACAt                                                                                                                                >  1:1040260/1‑151 (MQ=255)
                 aaaGAAGTCGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTg                                                                                                                              <  1:702036/151‑1 (MQ=255)
                 aaaGAAGTCGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTg                                                                                                                              <  1:1146191/151‑1 (MQ=255)
                    gaagTCGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATa                                                                                                                           <  1:424208/151‑1 (MQ=255)
                    gaagTCGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATa                                                                                                                           <  2:302102/151‑1 (MQ=255)
                      agTCGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGAt                                                                                                                            >  1:938886/1‑148 (MQ=255)
                      agTCGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATa                                                                                                                           >  2:224180/1‑149 (MQ=255)
                       gTCGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGc                                                                                                                        >  2:350785/1‑151 (MQ=255)
                       gTCGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGc                                                                                                                        >  2:32209/1‑151 (MQ=255)
                        tCGTTATTCAGTTTTTGAAACTGGGGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCCCAACATTGATa                                                                                                                           <  2:238211/147‑1 (MQ=255)
                         cGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCa                                                                                                                       <  2:249380/150‑1 (MQ=255)
                         cGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCa                                                                                                                       <  1:316025/150‑1 (MQ=255)
                         cGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCa                                                                                                                       >  1:59531/1‑150 (MQ=255)
                             aTTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGtt                                                                                                                   <  1:1164388/150‑1 (MQ=255)
                              ttCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTgc                                                                                                                 >  1:838099/1‑151 (MQ=255)
                                cAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTgcgc                                                                                                               <  2:458908/151‑1 (MQ=255)
                                cAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTgcgc                                                                                                               <  2:1159850/151‑1 (MQ=255)
                                  gTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCaa                                                                                                             <  1:1045319/151‑1 (MQ=255)
                                    ttttGAAACTGGTGACCGAGGCGGACGGCGTTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAAt                                                                                                            <  1:683051/150‑1 (MQ=255)
                                     tttGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATc                                                                                                           >  2:52794/1‑150 (MQ=255)
                                      ttGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCg                                                                                                          <  2:899325/150‑1 (MQ=255)
                                        gAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAAt                                                                                                       <  1:452993/151‑1 (MQ=255)
                                          aaCTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATcc                                                                                                     <  1:391104/151‑1 (MQ=255)
                                            cTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCgg                                                                                                   >  1:461272/1‑151 (MQ=255)
                                            cTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCgg                                                                                                   >  1:821537/1‑151 (MQ=255)
                                             tGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCgg                                                                                                   >  2:321557/1‑150 (MQ=255)
                                              ggTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGAt                                                                                                 <  1:771836/151‑1 (MQ=255)
                                               gTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGAtt                                                                                                <  2:936785/151‑1 (MQ=255)
                                                 gACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGt                                                                                              <  2:923042/151‑1 (MQ=255)
                                                     gAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTg                                                                                          >  1:89934/1‑151 (MQ=255)
                                                      aGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTg                                                                                          >  1:1195088/1‑150 (MQ=255)
                                                           gACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATcc                                                                                    <  1:187373/151‑1 (MQ=255)
                                                                cATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCa                                                                               <  1:960759/151‑1 (MQ=255)
                                                                cATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCa                                                                               <  2:21626/151‑1 (MQ=255)
                                                                 aTTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAg                                                                              <  2:468288/151‑1 (MQ=255)
                                                                  ttATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAgg                                                                             >  2:193449/1‑151 (MQ=255)
                                                                  ttATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAgg                                                                             >  2:272670/1‑151 (MQ=255)
                                                                   tATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGc                                                                            >  1:861257/1‑151 (MQ=255)
                                                                   tATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGc                                                                            >  1:263108/1‑151 (MQ=255)
                                                                    aTCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCt                                                                           >  1:443253/1‑151 (MQ=255)
                                                                    aTCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCt                                                                           <  2:105748/151‑1 (MQ=255)
                                                                        gTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTAt                                                                       >  2:221978/1‑151 (MQ=255)
                                                                                 ccTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGt                                                              <  2:534728/151‑1 (MQ=255)
                                                                                   tCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTgc                                                            >  2:114368/1‑151 (MQ=255)
                                                                                      aTGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTgctgg                                                         <  2:966266/151‑1 (MQ=255)
                                                                                          taaaaGCGGCAAAAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTAATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTgctgggct                                                      <  2:76285/148‑1 (MQ=255)
                                                                                           gAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTgctgggctgg                                                    <  1:1058704/151‑1 (MQ=255)
                                                                                             aaGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGctgggctggg                                                   >  1:222926/1‑150 (MQ=255)
                                                                                             aaGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGt                                                  <  2:1042232/151‑1 (MQ=255)
                                                                                                   cAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGAc                                               >  1:861945/1‑148 (MQ=255)
                                                                                                   cAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACaga                                            <  1:26826/151‑1 (MQ=255)
                                                                                                    aaGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAAACGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACaga                                            <  2:784164/150‑1 (MQ=255)
                                                                                                    aaGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACag                                             >  2:373356/1‑149 (MQ=255)
                                                                                                       agCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGaaa                                        <  1:863431/151‑1 (MQ=255)
                                                                                                           gAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATcc                                    >  1:337134/1‑151 (MQ=255)
                                                                                                            aaGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATcc                                    >  1:768875/1‑150 (MQ=255)
                                                                                                            aaGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCg                                   <  2:913603/151‑1 (MQ=255)
                                                                                                            aaGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCg                                   >  1:361231/1‑151 (MQ=255)
                                                                                                             aGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGc                                  >  1:25576/1‑151 (MQ=255)
                                                                                                             aGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCATGCCAGGCTGTATGCCCAAAGTGCCGGGCTGGGTGACAGAGAAAATCCGc                                  >  2:1044342/1‑151 (MQ=255)
                                                                                                               gTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGcc                                 >  2:668647/1‑150 (MQ=255)
                                                                                                                  ttttCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAAc                              <  1:634796/150‑1 (MQ=255)
                                                                                                                     tCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGTCTGGGTGACAGAGAAAATCCGCCAACCGCt                          >  1:661636/1‑151 (MQ=255)
                                                                                                                     tCTGTATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCt                          >  1:411891/1‑151 (MQ=255)
                                                                                                                         tATTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGAtt                      >  1:350948/1‑151 (MQ=255)
                                                                                                                          aTTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTg                     >  1:267863/1‑151 (MQ=255)
                                                                                                                          aTTCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTg                     >  1:256598/1‑151 (MQ=255)
                                                                                                                           ttCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGc                    >  2:859725/1‑151 (MQ=255)
                                                                                                                            tCATCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGcc                   >  1:307719/1‑151 (MQ=255)
                                                                                                                              aTCGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGCGACAGAGAAAATCCGCCAACCGCTGATTGCCgg                 >  2:207011/1‑151 (MQ=255)
                                                                                                                               tCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCgg                 >  2:680944/1‑150 (MQ=255)
                                                                                                                                cGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGt                <  1:273024/150‑1 (MQ=255)
                                                                                                                                cGCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTg               >  2:397169/1‑151 (MQ=255)
                                                                                                                                 gCCTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTgg              >  2:1204835/1‑151 (MQ=255)
                                                                                                                                  ccTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTggg             <  2:330194/151‑1 (MQ=255)
                                                                                                                                  ccTGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTggg             <  2:642565/151‑1 (MQ=255)
                                                                                                                                    tGTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCt           >  2:511669/1‑151 (MQ=255)
                                                                                                                                     gTTTATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCt           <  1:257935/150‑1 (MQ=255)
                                                                                                                                      tttATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAATTCCGCCAACCGCTGATTGCCGGTGGGCTgg         >  1:672362/1‑151 (MQ=255)
                                                                                                                                       ttATTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTgg         <  1:845293/150‑1 (MQ=255)
                                                                                                                                         aTTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGgtgt      <  2:7781/151‑1 (MQ=255)
                                                                                                                                         aTTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGgtg       <  2:895295/150‑1 (MQ=255)
                                                                                                                                         aTTGTTGATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGgtg       <  1:1139433/150‑1 (MQ=255)
                                                                                                                                          ttgttgATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGtgtg     >  2:500638/1‑151 (MQ=255)
                                                                                                                                            gttgATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGc    <  2:838056/150‑1 (MQ=255)
                                                                                                                                            gttgATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCg   <  2:509238/151‑1 (MQ=255)
                                                                                                                                            gttgATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCg   <  1:277953/151‑1 (MQ=255)
                                                                                                                                             ttgATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGTTGGTGTGCg   <  2:1071721/150‑1 (MQ=255)
                                                                                                                                              tgATGCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGa  <  2:567839/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGGAAGGCGCGTCAAATAAAGAAGTCGTTATTCAGTTTTTGAAACTGGTGACCGAGGCGGACGGCATTATCAGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGA  >  NC_000913/2894386‑2894677

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