Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,897,558 G→A 7.5% A104V (GCC→GTC)  ygcS ← putative MFS sugar transporter; membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,897,5580GA 7.5% 272.7 / 3.6 104A104V (GCC→GTC) ygcSputative MFS sugar transporter; membrane protein
Reads supporting (aligned to +/- strand):  ref base G (50/43);  new base A (2/6);  total (55/49)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 108 total aligned reads displayed.

TCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAGCGCCGAGCCGCCGATCATTCCCGCCAT  >  NC_000913/2897412‑2897706
                                                                                                                                                   |                                                                                                                                                    
tCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAACTTTA‑AAACCTGAGTTTGAACCTACTGAATAATCGACTCCCAGACCAATGCCAATCAAAAGCCGCAGGCCAATAAGATTCACTGGCGtggtgggc                                                                                                                                                    >  1:794061/1‑147 (MQ=255)
     accccGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCaaaaaa                                                                                                                                               <  1:113947/149‑1 (MQ=255)
       cacGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAAtt                                                                                                                                             >  2:672029/1‑149 (MQ=255)
        acGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATCCACGCCGGCCAAACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTg                                                                                                                                            <  1:346339/149‑1 (MQ=255)
         cGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAACCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGt                                                                                                                                           <  1:1081722/149‑1 (MQ=255)
          gCTAAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAAGTCAGCCAGCAAGGTGTG‑ACCTACTGGATAATCGCCTCCCAGACCAATGCCCACAACAACAACCCGATCGGCAAGATGCTCTGGCGTGATGCAAAAAAATTGTa                                                                                                                                          >  2:1732894/1‑149 (MQ=255)
                   gcCCAGCAAAATACCGCGATGGCGGCGCGGGGCACATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAGAAATTGTAAAAACGaag                                                                                                                                 >  2:835876/1‑149 (MQ=255)
                   gcCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAACATTGTAAAAACGaag                                                                                                                                 >  1:408568/1‑149 (MQ=255)
                     ccAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGgagca                                                                                                                               >  2:162274/1‑149 (MQ=255)
                      cAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGaagcaa                                                                                                                              >  1:1223906/1‑149 (MQ=255)
                       aGCAAAATGCCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGAGCG‑ACCTACTGAATAATCGCCTCCCATACCAATGCCAATCAAAATGCGCAGTCCACTAAGATCCCCTGGCCGGGTTGCAAAAAATGGTGAAAACGaagcagg                                                                                                                             >  2:981020/1‑147 (MQ=255)
                        gCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGTTGGCAAAAAATTGTAAAAACGaagcaagc                                                                                                                            >  1:109411/1‑149 (MQ=255)
                        gCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGaagcaagc                                                                                                                            >  1:541310/1‑149 (MQ=255)
                          aaaaTACCGCGATGGCGGCGCGGGGGAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCACTCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGaagcaagcat                                                                                                                          >  1:489353/1‑147 (MQ=255)
                          aaaaTACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGt                                                                                                                          >  1:311914/1‑149 (MQ=255)
                          aaaaTACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGt                                                                                                                          >  1:939845/1‑149 (MQ=255)
                          aaaaTACCGCGATGGCGGCACGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGACAAGCAAAATGCGCAGTCCAAGTAGATGCTCAGGTGTGGTGACAAAAAAATGCAAAAACAaaacaagaga                                                                                                                          >  2:499491/1‑146 (MQ=255)
                            aagaCCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAAGGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAACATTGTAAAAACGAAGCAAGCGTaa                                                                                                                        <  2:850755/146‑1 (MQ=255)
                              tACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTG‑ATCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCATAAAATTGTAAAAACGAAGCAAGCGTAATc                                                                                                                      <  1:1274968/149‑1 (MQ=255)
                              tACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAACATTGTAAAAACGAAGACAGCGTAATc                                                                                                                      >  2:1245674/1‑149 (MQ=255)
                               aCCGCGATGGCGGCGCGGGGAAAGTGCAGCCAGCAAGGGGTG‑ACCTACTGAAGAATCGCCTCCCAGACCGATGCCAATCAAAAGGCGCTGTCCAATAAGATGCTCTGGGGTGGTGGCAAAAAATTTTTAAAACGAAGCAAGCGTAATCa                                                                                                                     <  2:490184/149‑1 (MQ=255)
                               aCCGCGATGGCGGCGCGGGGAAAATTCATCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCa                                                                                                                     <  2:1728191/149‑1 (MQ=255)
                                ccGCGATGGCGGCGGTGGGAAAATTCAGCCCGTAAGGGGTG‑ACCTACTGAAGAAGCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCATAGAAGAAAAAAATTGTAAAAACGAAGCAAGCGTAATCAg                                                                                                                    <  2:1250787/149‑1 (MQ=255)
                                ccGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGAGGTGCCACAAAATTGTAAAAACGAAGCAAGCGGAATCAg                                                                                                                    >  2:1444657/1‑149 (MQ=255)
                                 cgcgATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGACCAAGAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAACCGCAATCAGc                                                                                                                   >  1:870305/1‑149 (MQ=255)
                                      tGGCGGCGCGGGGAAAATTCCGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCaaaaa                                                                                                              <  2:1470726/149‑1 (MQ=255)
                                          ggcgCGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAACTg                                                                                                          >  2:1431533/1‑149 (MQ=255)
                                             gcgGAGAAAATTCAACAAGGCAGGGGGG‑ACCCACAGAATAATCGCCTCCCAGACCAATCCCCATCAAAATTCGCAGTCCAATTACTTGCTCTGAAAAAAAAAAAAAAAATTGTAAAAGCGAAGCAAGCGTAATCAGCAAAAAGCTGAAg                                                                                                       <  2:1324332/149‑1 (MQ=255)
                                              cgGGAAATATGCAGCCAGCAAGGTGTG‑ACCTACTGACTAACCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAgg                                                                                                      <  2:818940/149‑1 (MQ=255)
                                                gggAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAgggg                                                                                                    >  1:153379/1‑147 (MQ=255)
                                                 ggAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGa                                                                                                   <  1:279571/149‑1 (MQ=255)
                                                 cgAAAATAAAGCCAGAAAGGTATG‑ATCCGCAGAATAACCGGCCCCCAGCCCAATGCCAATGAAAATGCGCAGTCCAATAAGATGCTTTGAAAAAAAAAAAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGa                                                                                                   <  2:1503078/148‑1 (MQ=255)
                                                  gtgaaTTCCGCCGCCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGCCCGATGTCGATCCAGATGCGCAGCCCAATACCACGCGCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGaa                                                                                                  <  2:489165/146‑1 (MQ=255)
                                                  gagatttCAGCCAGCATGGTGTG‑GCCTACTGAAGAACCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCTTAAAGCTGAAGGTGaa                                                                                                  <  2:667391/144‑1 (MQ=255)
                                                    aaagacaggcagcAGGGGGGG‑ACCTACAGAATAAACGCCTCCCAGACGAGTGCCAATCAAAATGCGCAGGCGAGTTAGATGCTGTGACGAGGAGAAAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGa                                                                                                <  2:1467152/144‑1 (MQ=255)
                                                         cagccagcAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGGAATCAGCAAAAAGCTGAAGGTGAAGAtttct                                                                                           >  2:626050/1‑147 (MQ=255)
                                                           gccagcAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAATAACGAAGCAAGCGTAATCGGCCAAAAGCTGAAGGTGAAGATTTTTTg                                                                                         <  1:1014161/149‑1 (MQ=255)
                                                            ccagcAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCGAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGCTTTTTTGc                                                                                        >  1:1178458/1‑149 (MQ=255)
                                                            ccagcAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATGAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGGGAAGATTTTTTGc                                                                                        >  1:1093484/1‑149 (MQ=255)
                                                            ccagcAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATCAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGc                                                                                        >  2:872211/1‑149 (MQ=255)
                                                            ccagcAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGc                                                                                        <  1:405893/149‑1 (MQ=255)
                                                             cagcAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTACTCAGCAAAAAGCTGAAGGTGAAGATTTTTTGcc                                                                                       >  1:649448/1‑149 (MQ=255)
                                                               gcAAGGTGTG‑ACGTACGGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGACGCTCTGGAGAAAAAAAAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGa                                                                                     <  2:1293725/149‑1 (MQ=255)
                                                                    gggtg‑ACCTACTGACTAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAat                                                                                <  2:1321115/147‑1 (MQ=255)
                                                                     tgtg‑ACCTACTGAATACTCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGACAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAata                                                                               <  1:1402438/149‑1 (MQ=255)
                                                                     tgtg‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAata                                                                               <  1:261326/149‑1 (MQ=255)
                                                                       tg‑ACCTACTGAATAATCGCCTCCCAGACCAGTGAGGTACAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATGAGCAAAAAGATGAAGGTGAAGATTTTTTGCGGACCAATATg                                                                             >  2:1747150/1‑149 (MQ=255)
                                                                       tg‑ACCTACTGAATAATCGCCTCCCAGACCACTGACGCCCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATGGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGAGAAGATCATTTGCCGACCAatacg                                                                             >  2:1747627/1‑147 (MQ=255)
                                                                       tg‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAGAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATg                                                                             >  1:1400769/1‑149 (MQ=255)
                                                                       tg‑ACCTACTCAATCATCGCCCCCCAGAGCAAGGCCAAGCACAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTCTTTGCCGACCAATATg                                                                             <  2:1367453/149‑1 (MQ=255)
                                                                          aCCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGt                                                                           <  2:1035671/149‑1 (MQ=255)
                                                                          aCCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGt                                                                           >  1:1362341/1‑149 (MQ=255)
                                                                            cTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCGAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATAAGGTCg                                                                         >  1:797646/1‑149 (MQ=255)
                                                                            cTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCg                                                                         >  2:1315675/1‑149 (MQ=255)
                                                                            cTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCg                                                                         >  1:1038319/1‑149 (MQ=255)
                                                                            cTACTGAATAATCCCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCg                                                                         <  2:297173/149‑1 (MQ=255)
                                                                             tACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTAGCCGACCAAGATGGTCgg                                                                        >  1:1131450/1‑149 (MQ=255)
                                                                              aCTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGga                                                                       <  1:1186452/149‑1 (MQ=255)
                                                                              aCTGAACAATCGCCTCCCAGAACAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGga                                                                       <  1:733235/149‑1 (MQ=255)
                                                                                      atCGCCTCCCAGACCAATGTCAATCAAAATGCGCCGTCCAATAAGATGCTCTGGTGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCGGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCAcc                                                               <  2:386328/149‑1 (MQ=255)
                                                                                       tCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAACAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCAccc                                                              <  1:1565146/149‑1 (MQ=255)
                                                                                       gcGCCTCCCAGACCAATGCCAAGCAAAATGCGCCGTCCAATGAGCTGCTCAGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCAccc                                                              <  1:141176/148‑1 (MQ=255)
                                                                                        cGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATCGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATAGTCGGAGATCCACTCa                                                             >  2:1219255/1‑149 (MQ=255)
                                                                                         gCCTCCCAGACCCATGGCAATCCAAATGCGCAGTCCAATAAGATGCTTTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCACGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCaa                                                            <  2:703842/149‑1 (MQ=255)
                                                                                          ccTTCCAGACCAACGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAg                                                           <  1:1493835/149‑1 (MQ=255)
                                                                                             cccAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAAc                                                        <  1:384257/149‑1 (MQ=255)
                                                                                              ccAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCACAAAGCTGAAGGTGAAGA‑TTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCa                                                      <  1:957694/149‑1 (MQ=255)
                                                                                               ccgagcAAGGCCAATCAAAAGGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGAAAAAAAAAAAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCa                                                      <  2:1623248/144‑1 (MQ=255)
                                                                                               cAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCGCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCa                                                      <  1:313010/149‑1 (MQ=255)
                                                                                                aGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAAcccg                                                     >  1:1081785/1‑147 (MQ=255)
                                                                                                       tggcaaacaaaATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCa                                              >  2:311770/4‑149 (MQ=255)
                                                                                                       tGCCAATCAAAATGCGCAGCCCAATAAGATGCTCTGGCGTGGTGACAAAAAATTGTAAAAACCACGCCAGCGTCCTCAGCAAAAAGATGAAGGTGAAGATTTCTTGCCGACCAATATGGTCGGAGACACACCCAAGAACCAGGATGCCa                                              >  2:316938/1‑149 (MQ=255)
                                                                                                        gCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAAATCGACGCAAGAACCAGCCTGcaag                                             >  2:1575211/1‑146 (MQ=255)
                                                                                                        gCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGGGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCaa                                             <  1:726602/149‑1 (MQ=255)
                                                                                                        gCCAATCAAAATGCGCAGTCCAATAAGATGCGCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCaa                                             <  1:1487077/149‑1 (MQ=255)
                                                                                                           aaTCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGa                                          >  2:1457331/1‑149 (MQ=255)
                                                                                                           aaTCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGa                                          >  1:239271/1‑149 (MQ=255)
                                                                                                           aaTCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGa                                          >  2:1417231/1‑149 (MQ=255)
                                                                                                            aTCAAAATGCGCAGTCCAATGAGACGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCACGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGaa                                         <  2:1222860/149‑1 (MQ=255)
                                                                                                            aTCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGaa                                         <  2:91699/149‑1 (MQ=255)
                                                                                                             tCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTCAAAACGAAGCAAGCGTAGTCAGCAAAGAGCTGAAGGTGAAGATGATTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGACAAGGGAc                                        >  1:349406/1‑149 (MQ=255)
                                                                                                             tCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCGAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAAc                                        >  1:1078924/1‑149 (MQ=255)
                                                                                                                aaaTGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATAGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACaaa                                     >  1:1676939/1‑149 (MQ=255)
                                                                                                                   tGCGCAGTCCAATCAGATGCCTTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCGAGCGTAATCAGCAAAAAGCCGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAACCCa                                  <  1:1999/149‑1 (MQ=255)
                                                                                                                     cgcAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAg                                <  1:9363/149‑1 (MQ=255)
                                                                                                                         gTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAgc                            >  1:1648870/1‑149 (MQ=255)
                                                                                                                          tCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTGGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGCACAAACCAAGGAgcg                           >  1:1256656/1‑149 (MQ=255)
                                                                                                                           ccAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAgcgc                          <  1:302517/149‑1 (MQ=255)
                                                                                                                            cAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAGCGcc                         >  1:765296/1‑149 (MQ=255)
                                                                                                                             aaTAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTTTCGTTTAATATGGTCGGAGATCCACCCAAGAACCAGGCTACCAAGGAACAAACAAAGGAGCGCCg                        >  2:1766454/1‑149 (MQ=255)
                                                                                                                               tAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAGCGCCGAg                      >  1:720893/1‑149 (MQ=255)
                                                                                                                                   aggccctGGCGTGGTGACAAAAAATGGTAAAAACGAAGCAAGAAAAAAAAAAAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAGCGCCGAgccgc                  <  2:700126/144‑1 (MQ=255)
                                                                                                                                   aTGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTTCCAAGGAACAAACCAAGGAGCGCAGAgccgc                  >  1:179321/1‑149 (MQ=255)
                                                                                                                                   aTGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAGCGCCGAgccgc                  <  2:715384/149‑1 (MQ=255)
                                                                                                                                   aTGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAGCGCCGAgccgc                  >  1:1623209/1‑149 (MQ=255)
                                                                                                                                      ctctGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAGCGCCGAGCCGCCGa               >  1:826256/1‑149 (MQ=255)
                                                                                                                                           gCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAGCGCCGAGCCGCCGATCAtt          >  1:47775/1‑149 (MQ=255)
                                                                                                                                                  ggCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAGCGCCGAGCCGCCGATCATTCCCGCCa   >  2:1606755/1‑149 (MQ=255)
                                                                                                                                                   gCAAAAAATTGTAAAACCGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAGCGCCGAGCCGCCGATCATTCCCGCCAt  <  2:541310/149‑1 (MQ=255)
                                                                                                                                                   gCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATACGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAGCGCCGAGCCGCCGATCATTCCCGCCAt  >  1:969608/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTG‑ACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAGCGCCGAGCCGCCGATCATTCCCGCCAT  >  NC_000913/2897412‑2897706

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