Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,859,354 T→C 46.6% N255D (AAC→GAC)  yidK ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,859,3540TC46.6% 474.4 / inf 2001N255D (AAC→GAC) yidKputative transporter
Reads supporting (aligned to +/- strand):  ref base T (564/504);  new base C (486/445);  total (1052/949)
Fisher's exact test for biased strand distribution p-value = 7.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

Only 100 of 2043 total aligned reads displayed.

AAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGGCCAAATACCG  >  NC_000913/3859205‑3859503
                                                                                                                                                     |                                                                                                                                                     
  aCGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCa                                                                                                                                                     >  2:1834711/1‑150 (MQ=255)
  aCGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCa                                                                                                                                                     >  1:11238587/1‑150 (MQ=255)
  aCGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCa                                                                                                                                                     >  2:6277702/1‑150 (MQ=255)
    gCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCAcc                                                                                                                                                   >  2:7754647/1‑150 (MQ=255)
    gCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCAcc                                                                                                                                                   <  1:7220045/150‑1 (MQ=255)
    gCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCAcc                                                                                                                                                   >  1:2730984/1‑150 (MQ=255)
     cAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTccacca                                                                                                                                                  >  1:1166762/1‑150 (MQ=255)
       aTCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGa                                                                                                                                                >  1:3135119/1‑150 (MQ=255)
        tCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGaa                                                                                                                                               <  2:5637152/150‑1 (MQ=255)
         cAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAAt                                                                                                                                              <  2:7531288/150‑1 (MQ=255)
           aCCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCGTTTGCCCTTCCGCCAGGCGTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCa                                                                                                                                            >  2:4547538/1‑150 (MQ=255)
              cTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAAc                                                                                                                                         >  1:5108157/1‑150 (MQ=255)
                      aCCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAAt                                                                                                                                 >  2:9100784/1‑150 (MQ=255)
                        cAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGc                                                                                                                               >  2:10342368/1‑150 (MQ=255)
                        cAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGc                                                                                                                               >  1:5554565/1‑150 (MQ=255)
                        cAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGc                                                                                                                               >  1:6308324/1‑150 (MQ=255)
                         aGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGcc                                                                                                                              >  1:111024/1‑150 (MQ=255)
                             ccAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCgcgc                                                                                                                          >  1:12072935/1‑150 (MQ=255)
                              cAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGcc                                                                                                                         >  2:9330799/1‑150 (MQ=255)
                                   ggTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCg                                                                                                                    <  2:9497307/150‑1 (MQ=255)
                                   ggTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCgg                                                                                                                   >  1:7912882/1‑150 (MQ=255)
                                      cGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCa                                                                                                                 <  1:9764528/150‑1 (MQ=255)
                                       gAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCGTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCaa                                                                                                                >  2:4709437/1‑150 (MQ=255)
                                         gCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAgg                                                                                                              <  1:12231833/150‑1 (MQ=255)
                                            ttttCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGAt                                                                                                           <  1:4848531/150‑1 (MQ=255)
                                             tttCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTTACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGAtc                                                                                                          >  1:11877931/1‑150 (MQ=255)
                                                cAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGAtcggt                                                                                                       <  1:9073016/150‑1 (MQ=255)
                                                 aGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACGCCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGAtcggtc                                                                                                      <  2:4387980/150‑1 (MQ=255)
                                                 aGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCCGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGAtcggtc                                                                                                      >  2:3794633/1‑150 (MQ=255)
                                                 aGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGAtcggtc                                                                                                      <  2:11548812/150‑1 (MQ=255)
                                                 aGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGAtcggtc                                                                                                      >  1:1733524/1‑150 (MQ=255)
                                                 aGCACCGCCGTTAACAGCGCCCCCGTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGAtcggtc                                                                                                      >  1:1414785/1‑150 (MQ=255)
                                                   cACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGAtcggtcgg                                                                                                    <  1:6131593/150‑1 (MQ=255)
                                                    aCCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCGTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGAGTTGTACACCAGTAAAAGGTGTCCACCAGAATCACACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCggg                                                                                                   >  2:11381485/1‑150 (MQ=255)
                                                     ccgccgTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGc                                                                                                  >  2:4196662/1‑150 (MQ=255)
                                                      cgccgTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGcc                                                                                                 >  2:240965/1‑150 (MQ=255)
                                                       gccgTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCa                                                                                                <  1:2580978/150‑1 (MQ=255)
                                                        ccgTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCAc                                                                                               >  1:2161135/1‑150 (MQ=255)
                                                        ccgTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCAc                                                                                               <  1:3668620/150‑1 (MQ=255)
                                                        ccgTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCAc                                                                                               <  2:1402847/150‑1 (MQ=255)
                                                         cgTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCAcc                                                                                              >  2:11783828/1‑150 (MQ=255)
                                                         cgTTAACAGCGCCCCCTTGTGCCCGTCCGCCAGGCTTTTGGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCAcc                                                                                              <  2:8250071/150‑1 (MQ=255)
                                                          gTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCg                                                                                             >  1:2410514/1‑150 (MQ=255)
                                                                aGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAg                                                                                       >  1:8478415/1‑150 (MQ=255)
                                                                aGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAg                                                                                       >  1:3499223/1‑150 (MQ=255)
                                                                  cgcCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGtt                                                                                     <  1:8841172/150‑1 (MQ=255)
                                                                      cccTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAAt                                                                                 <  1:387979/150‑1 (MQ=255)
                                                                         ttttGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTtc                                                                              <  2:8858582/150‑1 (MQ=255)
                                                                          tttGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTtct                                                                             <  1:411584/150‑1 (MQ=255)
                                                                            tGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCg                                                                           <  2:3170083/150‑1 (MQ=255)
                                                                            tGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCg                                                                           <  1:6520197/150‑1 (MQ=255)
                                                                               ccgtCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAGGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCg                                                                        <  2:499257/147‑1 (MQ=255)
                                                                               ccTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCg                                                                        >  1:8198043/1‑150 (MQ=255)
                                                                               ccTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCg                                                                        >  2:8474731/1‑150 (MQ=255)
                                                                               ccTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCg                                                                        <  2:5291941/150‑1 (MQ=255)
                                                                                cTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGt                                                                       >  1:9542690/1‑150 (MQ=255)
                                                                                cTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGt                                                                       >  2:8654939/1‑150 (MQ=255)
                                                                                 ttCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTg                                                                      >  1:2978273/1‑150 (MQ=255)
                                                                                  tCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGCCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGa                                                                     <  1:10881257/150‑1 (MQ=255)
                                                                                    cGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAAc                                                                   <  2:2854081/150‑1 (MQ=255)
                                                                                            tttttGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGa                                                           <  1:6698440/150‑1 (MQ=255)
                                                                                             ttttGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAg                                                          <  1:8624500/150‑1 (MQ=255)
                                                                                               ttGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGtt                                                        <  2:12147882/150‑1 (MQ=255)
                                                                                                tGACGCCAGCGTGCGTTGCACGATGCCCTGATGTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTg                                                       <  2:10152213/150‑1 (MQ=255)
                                                                                                  aCGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCt                                                     >  2:11970791/1‑150 (MQ=255)
                                                                                                  aCGCCAGCGTGCGTTGCACGATGCCCTGATTTGGACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGTACGGTGGTGAGTTGCt                                                     <  1:6880623/150‑1 (MQ=255)
                                                                                                   cGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTc                                                    <  2:3519848/150‑1 (MQ=255)
                                                                                                     ccAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCaa                                                  <  2:302897/150‑1 (MQ=255)
                                                                                                      cAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAAt                                                 <  2:9686811/150‑1 (MQ=255)
                                                                                                      cAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAAt                                                 >  2:2586998/1‑150 (MQ=255)
                                                                                                       aGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATg                                                <  2:334236/150‑1 (MQ=255)
                                                                                                         cgtgcgtTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGcc                                              >  1:10415985/1‑150 (MQ=255)
                                                                                                              gtTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCa                                         <  2:6572432/150‑1 (MQ=255)
                                                                                                              gtTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCa                                         <  1:7107020/150‑1 (MQ=255)
                                                                                                                   aCGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAg                                    >  1:10027092/1‑150 (MQ=255)
                                                                                                                   aCGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAg                                    <  1:1984137/150‑1 (MQ=255)
                                                                                                                   aCGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAg                                    <  1:156252/150‑1 (MQ=255)
                                                                                                                    cGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGc                                   <  2:5906700/150‑1 (MQ=255)
                                                                                                                     gATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCt                                  <  2:1625976/150‑1 (MQ=255)
                                                                                                                     gATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCt                                  >  1:9332841/1‑150 (MQ=255)
                                                                                                                      aTGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTg                                 <  1:11538867/150‑1 (MQ=255)
                                                                                                                        gCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGcc                               <  1:10220739/150‑1 (MQ=255)
                                                                                                                         cccTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGccc                              >  1:10660294/1‑150 (MQ=255)
                                                                                                                          ccTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCt                             >  2:5402806/1‑150 (MQ=255)
                                                                                                                          ccTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCt                             >  1:7439670/1‑150 (MQ=255)
                                                                                                                           cTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCtt                            >  2:7974307/1‑150 (MQ=255)
                                                                                                                            tGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTg                           <  2:2161135/150‑1 (MQ=255)
                                                                                                                            tGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTg                           <  2:108990/150‑1 (MQ=255)
                                                                                                                             gATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGc                          <  2:10615379/150‑1 (MQ=255)
                                                                                                                             gATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGc                          <  2:11187612/150‑1 (MQ=255)
                                                                                                                                     cacCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATcgcg                  >  2:723767/1‑150 (MQ=255)
                                                                                                                                      acCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGa                 >  1:7681672/1‑150 (MQ=255)
                                                                                                                                        cAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGAtt               >  1:4328343/1‑150 (MQ=255)
                                                                                                                                        cAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGAtt               <  1:5857396/150‑1 (MQ=255)
                                                                                                                                         aGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTa              <  1:459179/150‑1 (MQ=255)
                                                                                                                                         aGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTa              <  1:3324198/150‑1 (MQ=255)
                                                                                                                                         aGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTa              <  1:5368240/150‑1 (MQ=255)
                                                                                                                                         aGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTa              >  2:3997640/1‑150 (MQ=255)
                                                                                                                                            aaaaGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGGc           <  2:7033278/150‑1 (MQ=255)
                                                                                                                                             aaaGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGGcc          <  2:6823668/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGGCCAAATACCG  >  NC_000913/3859205‑3859503

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