Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,845,824 T→C E448E (GAA→GAG uhpT ← hexose phosphate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,845,8240TC92.9% 290.7 / 3.9 99E448E (GAA→GAGuhpThexose phosphate transporter
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base C (77/14);  total (81/18)
Fisher's exact test for biased strand distribution p-value = 2.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 102 total aligned reads displayed.

GGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCAAAGCTGTC  >  NC_000913/3845676‑3845970
                                                                                                                                                       |                                                                                                                                                  
ggggGTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGttc                                                                                                                                                  <  4:1567297/151‑1 (MQ=255)
 gggCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGctct                                                                                                                                                 >  1:1150998/1‑151 (MQ=255)
 aggCTATTTGATAGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGctct                                                                                                                                                 >  4:1859058/2‑151 (MQ=255)
    cTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTcc                                                                                                                                              >  4:800629/1‑151 (MQ=255)
     tATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTcc                                                                                                                                              >  4:399535/1‑150 (MQ=255)
      aTTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCAt                                                                                                                                            >  1:1008473/1‑151 (MQ=255)
        ttGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCAtt                                                                                                                                           <  4:91374/150‑1 (MQ=255)
        ttGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCAtt                                                                                                                                           >  1:156305/1‑150 (MQ=255)
        ttGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCAtt                                                                                                                                           >  4:265941/1‑150 (MQ=255)
         tGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTAc                                                                                                                                         >  1:647781/1‑151 (MQ=255)
          gATGGCGAGAAGTTTGCCTTCACTACGTTGGAAGTCACGGAATAATAAAAAAAGCCAGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGATCTTCCATTa                                                                                                                                          >  4:1770174/1‑149 (MQ=255)
           aTGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCg                                                                                                                                       <  2:44622/151‑1 (MQ=255)
              gCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCa                                                                                                                                    <  4:1041950/151‑1 (MQ=255)
               cGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGc                                                                                                                                      >  1:1437464/1‑148 (MQ=255)
                  gaAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTAt                                                                                                                                >  4:338559/1‑151 (MQ=255)
                    aGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCg                                                                                                                              >  2:545447/1‑151 (MQ=255)
                      tttGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGcc                                                                                                                            >  1:467075/1‑151 (MQ=255)
                          ccTTCCCTGCGCGGAAAGTCACGAAATAAAAAAAAACGCCCGGCGTCATGCCGGAAAAAAGTCCCCATTTATTTTTATGCCACTTTCAACTGCTGATTTTTTTTCCCGCGGGGTTTTTG‑‑GCGCTCTTCCATTACCGCCACTATCGCCAt                                                                                                                          <  1:568712/149‑1 (MQ=255)
                              cACTACGCTGAAAGTCACGAAATAAAAAAAAAAGCCCGGCGTCAGGCCGGAAAAAAGCCCCCAGTTACTTTTATCCCCCTGTCAACTGTTAATTTTTTTTCCCGCGGGGGTTTTT‑‑GCCCTCTCCCATTACCGCCACTATCGCCATCAGGc                                                                                                                     <  4:54380/150‑1 (MQ=255)
                                cTACGCGGGAACTAACGAAAAAAAAAAAAAGCCCCGGCGTCATGCCGGGCAAAAGCCCCCAGTTACTTTTATCCCCCTGTCAATTGCTAATTTTTTTTCCCGCGGCGTTTTTT‑‑GCCCTTTCCCATTACCGCCACTATCGCCATCAGGCa                                                                                                                    <  4:895119/149‑1 (MQ=255)
                                  aCGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGAt                                                                                                                 >  4:29555/1‑150 (MQ=255)
                                   cGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATa                                                                                                                >  1:1223030/1‑150 (MQ=255)
                                   cGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATAc                                                                                                               >  3:103869/1‑151 (MQ=255)
                                    gCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATAc                                                                                                               >  4:1797001/1‑150 (MQ=255)
                                      tGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAAc                                                                                                            >  1:333304/1‑151 (MQ=255)
                                             cACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATcgc                                                                                                     >  1:316100/1‑151 (MQ=255)
                                               cGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCgcgg                                                                                                   >  1:1493529/1‑151 (MQ=255)
                                                 aaaaaaaaaaaaagGCCCGGCGCCAGGCCGGAAAAAAGCCCCCATTTTGTTTTACCCCATTACCACCTCTTATTTTTTTTCCCGGGGGCGTTTTTC‑‑CCCTTCTTCCATTTCCGCCACTATCGCCATCAGGCAGATACAACCAATCgcggcg                                                                                                 <  1:1812587/144‑1 (MQ=255)
                                                  aataatAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCgcggcg                                                                                                 >  4:294640/1‑150 (MQ=255)
                                                  aaaaaaaaaaaaaGCCGGGCGTCATGCGGGGAAAAAGCCCCCAGTTCCTTTTATCCCACTGTCACCTGCTAATTTTTTTTCCGGGGGCGTTTTTT‑‑GCGCTCTCCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCgcggcg                                                                                                 <  1:491818/144‑1 (MQ=255)
                                                    taataaAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGat                                                                                               >  1:604869/1‑150 (MQ=255)
                                                    taataaAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGat                                                                                               >  1:142069/1‑150 (MQ=255)
                                                     aaaaaaaaaaGCCGGGCTTCAGGCGGGGAAAAAGCCACCATTTCCTTTTACCCCACTGTAAACTCTTTATTTTTTTCCCCGGGGCGATTTTC‑‑CCCCTTTCCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGatat                                                                                             <  4:1712283/148‑1 (MQ=255)
                                                      ataaAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑TCGCTCTTCCATTACCGCCACTATCGTCATCAGGCAGATACAACCAATCGCGGCGATATc                                                                                            >  4:29207/1‑151 (MQ=255)
                                                         aaaaagCCCCGGCGCAAGGCGGGGAAAAGGCCCCCATTTACTTTTATCCCACTGTCAACTGTTGTTTTTTTTTCGCGGGGGGTTTTTC‑‑GCGCTCTCCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAg                                                                                         <  1:940106/151‑1 (MQ=255)
                                                           aaaaGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAgcg                                                                                       >  2:184502/1‑151 (MQ=255)
                                                            aaaGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAgcgc                                                                                      >  1:1089513/1‑151 (MQ=255)
                                                             aaGCCCGGCGTCATGCCGGGCAAAAGTCCCCAGTTACGCTTTTGCCACTGTCCACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GAGCTCTTCCATTACCTCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAgc                                                                                        >  1:1544809/1‑148 (MQ=255)
                                                             aaGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGcgcg                                                                                     >  1:1462771/1‑151 (MQ=255)
                                                                cccGGCTTCACGCCGGGAAAACTCCCCCTTTTCCTTTTATCCCCCTGTCACCTGCTAATTTTTTTCCCCGCGGGGGTTTTT‑‑GCCCTTTTCCATTGCCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCgcgg                                                                                    <  4:292322/149‑1 (MQ=255)
                                                                cccGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTGCCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCgcggcg                                                                                  >  4:273695/1‑151 (MQ=255)
                                                                 ccGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCgcggcg                                                                                  >  2:504523/1‑150 (MQ=255)
                                                                    gCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAgg                                                                              >  1:405760/1‑151 (MQ=255)
                                                                      gTCAGGCCGGAAAAAAGCCCCCAGTTCCTTTTACGCCCCTGTCAACTGCTAATTTTTTTTCCCGGGGGGTTTTTT‑‑GCCCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAg                                                                               <  1:1196780/148‑1 (MQ=255)
                                                                         aTGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGtgcct                                                                         >  1:902558/1‑151 (MQ=255)
                                                                         aTGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGtgcct                                                                         >  3:422433/1‑151 (MQ=255)
                                                                           gcgggaaaaaatcccccATTTCCGTTTATCCCACTGTCAACTGCTGATTTTTTTTCCCGGGGGGGTTTTT‑‑GCCCTCTTCCATTCCCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGtgcctg                                                                        <  4:710917/135‑1 (MQ=255)
                                                                           gCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGtgcctgc                                                                       >  1:1887719/1‑151 (MQ=255)
                                                                             cGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGccc                                                                     >  4:973365/1‑151 (MQ=255)
                                                                             cGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGccc                                                                     >  4:285627/1‑151 (MQ=255)
                                                                                gCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAg                                                                   >  4:483441/1‑150 (MQ=255)
                                                                                  aaaaGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCg                                                                >  1:744569/1‑151 (MQ=255)
                                                                                  aaaaGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCg                                                                >  2:313906/1‑151 (MQ=255)
                                                                                   aaaGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCgg                                                               >  1:850766/1‑151 (MQ=255)
                                                                                     aGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGt                                                              >  4:653903/1‑150 (MQ=255)
                                                                                       tCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTa                                                             >  1:1456855/1‑149 (MQ=255)
                                                                                        cACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAgg                                                          >  1:1764289/1‑151 (MQ=255)
                                                                                           cAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCg                                                       >  4:1422445/1‑151 (MQ=255)
                                                                                               tACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAAt                                                    >  4:875669/1‑150 (MQ=255)
                                                                                                aCGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATa                                                   >  4:1300171/1‑150 (MQ=255)
                                                                                                aCGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATAc                                                  >  1:958107/1‑151 (MQ=255)
                                                                                                  gTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCg                                                >  4:395176/1‑151 (MQ=255)
                                                                                                   tttATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCgg                                               >  4:1720768/1‑151 (MQ=255)
                                                                                                     tATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGc                                              >  4:1304218/1‑150 (MQ=255)
                                                                                                     tATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCg                                             >  4:1863011/1‑151 (MQ=255)
                                                                                                     tATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCg                                             >  1:151782/1‑151 (MQ=255)
                                                                                                      aTGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGt                                            >  4:41108/1‑151 (MQ=255)
                                                                                                        gCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTcc                                          >  2:358560/1‑151 (MQ=255)
                                                                                                          cACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTccc                                         >  3:517726/1‑150 (MQ=255)
                                                                                                          cACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTccc                                         >  1:1219910/1‑150 (MQ=255)
                                                                                                          cACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTccc                                         >  3:269361/1‑150 (MQ=255)
                                                                                                           aCTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCAt                                       >  1:352610/1‑151 (MQ=255)
                                                                                                             tGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCg                                     >  1:214861/1‑151 (MQ=255)
                                                                                                                    tgctgAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCTATc                               >  4:1299650/1‑150 (MQ=255)
                                                                                                                     tgttattttttttCGCGGGGCGGTTTTG‑‑CCCCTCTCCCATCACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCaa                                 <  3:79046/141‑1 (MQ=255)
                                                                                                                     tctttaTTTTTTTCCGCGGGGGGTTTTTGGCCTTCTTCTTTTCCCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAAt                                <  4:1169655/144‑1 (MQ=255)
                                                                                                                        gAATTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTcc                          >  2:6027/1‑151 (MQ=255)
                                                                                                                         aaTTTTTTTCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTccc                         >  1:273966/1‑151 (MQ=255)
                                                                                                                           tttttttCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTccc                         >  1:358183/1‑149 (MQ=255)
                                                                                                                           tttttttCCCGGGGGGATTTTT‑‑CCCTGCTCCCTTAACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATc                               <  1:1631939/143‑1 (MQ=255)
                                                                                                                            ttttttCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGATACCGGCCGCCCATTGGGAATTCTTccc                         >  1:1785058/1‑147 (MQ=255)
                                                                                                                             tttttCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGAc                     >  3:115649/1‑151 (MQ=255)
                                                                                                                             tttttCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGAc                     >  4:1099360/1‑151 (MQ=255)
                                                                                                                             tttttCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGAc                     >  4:521104/1‑151 (MQ=255)
                                                                                                                               tttcccggcggggGTTTTT‑‑GCCCTTTCCTATTCCCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGAc                     <  4:186715/143‑1 (MQ=255)
                                                                                                                               tttCTCGCGGCGGATTTT‑‑GCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCt                   >  4:1859163/1‑151 (MQ=255)
                                                                                                                                ttCTCGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCt                   >  2:137448/1‑150 (MQ=255)
                                                                                                                                 tctcGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTa                  >  1:148803/1‑150 (MQ=255)
                                                                                                                                  ctcGCGGCGGATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAAc                >  1:535015/1‑151 (MQ=255)
                                                                                                                                        ggcggATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGc           >  4:11967/1‑151 (MQ=255)
                                                                                                                                           ggATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCaa         >  4:4179/1‑150 (MQ=255)
                                                                                                                                           ggATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCaa         >  2:435312/1‑150 (MQ=255)
                                                                                                                                           ggATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCaa         >  1:1862876/1‑150 (MQ=255)
                                                                                                                                           ggATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCaa         >  1:1110121/1‑150 (MQ=255)
                                                                                                                                           ggATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCa          >  1:377172/1‑149 (MQ=255)
                                                                                                                                            gATTTT‑‑GCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCaa         >  1:271435/1‑149 (MQ=255)
                                                                                                                                            gATTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCaa         >  4:1222178/1‑149 (MQ=255)
                                                                                                                                             aTTTT‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCAAAGc      >  2:388228/1‑151 (MQ=255)
                                                                                                                                                 t‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCAAAGCTGTc  >  1:358168/1‑151 (MQ=255)
                                                                                                                                                 t‑‑GCGCTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCAAAGCTGTc  >  2:41702/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
GGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT‑‑GCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCAAAGCTGTC  >  NC_000913/3845676‑3845970

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