Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 886,189 A→G 5.5% V384A (GTC→GCC)  ybjJ ← inner membrane protein YbjJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913886,1890AG 5.5% 391.0 / 9.0 127V384A (GTC→GCC) ybjJinner membrane protein YbjJ
Reads supporting (aligned to +/- strand):  ref base A (60/60);  new base G (2/5);  total (62/65)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

Only 100 of 127 total aligned reads displayed.

AGCCATAAACCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTGCAT  >  NC_000913/886048‑886336
                                                                                                                                             |                                                                                                                                                   
aGCCATAAACCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGcc                                                                                                                                             >  2:1800991/1‑150 (MQ=255)
 gCCATAAACCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCcg                                                                                                                                            >  1:1944254/1‑150 (MQ=255)
 gCCATAAACCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCcg                                                                                                                                            >  2:2141710/1‑150 (MQ=255)
  ccATAAACCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGGTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCcgc                                                                                                                                           >  2:2345743/1‑150 (MQ=255)
   cATAAACCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCcgcg                                                                                                                                          >  1:70462/1‑150 (MQ=255)
   cATAAACCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCcgcg                                                                                                                                          <  1:1315266/150‑1 (MQ=255)
    aTAAACCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCga                                                                                                                                         <  2:317778/150‑1 (MQ=255)
     tAAACCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCgag                                                                                                                                        <  1:144278/150‑1 (MQ=255)
      aaaCCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCgaga                                                                                                                                       <  1:226660/150‑1 (MQ=255)
       aaCCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGGCTATTGTCGCgaaag                                                                                                                                      <  2:1519936/150‑4 (MQ=255)
       aaCCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGGCTATTGTCGCgaaag                                                                                                                                      <  1:2133947/150‑4 (MQ=255)
        aCCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGGCTATTGTCGCgaaagc                                                                                                                                     <  1:666472/150‑5 (MQ=255)
        aCCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGGCTATTGTCGCgaaagc                                                                                                                                     <  2:1314188/150‑5 (MQ=255)
         ccGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATa                                                                                                                                    <  2:1481118/150‑1 (MQ=255)
         ccGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATa                                                                                                                                    >  1:296046/1‑150 (MQ=255)
            cTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAAcc                                                                                                                                 >  1:974445/1‑150 (MQ=255)
              cACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAg                                                                                                                               >  1:1864755/1‑150 (MQ=255)
                 ggTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAgcgc                                                                                                                            <  1:239972/150‑1 (MQ=255)
                    tataGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAg                                                                                                                         >  2:1943393/1‑150 (MQ=255)
                     ataGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGt                                                                                                                        >  2:1852894/1‑150 (MQ=255)
                      taGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTa                                                                                                                       <  1:2310979/150‑1 (MQ=255)
                      taGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTa                                                                                                                       <  2:643269/150‑1 (MQ=255)
                        ggTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACa                                                                                                                     >  2:1066741/1‑150 (MQ=255)
                        ggTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACa                                                                                                                     >  1:975066/1‑150 (MQ=255)
                         gTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGATTAACCAGCGCCAGTACaa                                                                                                                    >  2:1626332/1‑150 (MQ=255)
                         gTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACaa                                                                                                                    >  1:1622127/1‑150 (MQ=255)
                          ttttttGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAAc                                                                                                                   >  2:1227645/1‑150 (MQ=255)
                           tttttGATCGCTTAAGAAAGTACCATCCATTTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAAcc                                                                                                                  <  2:109258/150‑1 (MQ=255)
                           tttttGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAAcc                                                                                                                  <  2:322921/150‑1 (MQ=255)
                             tttGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAg                                                                                                                >  1:1732883/1‑150 (MQ=255)
                              ttGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGc                                                                                                               <  1:2300133/150‑1 (MQ=255)
                                gATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCAt                                                                                                             >  2:1704866/1‑150 (MQ=255)
                                gATCGCTTAAGAAAGTACCATCCATGTCAACCGAAAAAAATTTAATCCAAATCAAATATTCTCCATCGACGACTGCGATATTGTATCGGGTTTCGCGACGGCTTTCGCGACAATAGCCGCGACAATAACCAACGCCAGTACAACCAGCaa                                                                                                             >  2:1704274/1‑149 (MQ=255)
                                   cGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGc                                                                                                          >  2:1542505/1‑150 (MQ=255)
                                   cGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGc                                                                                                          <  1:344424/150‑1 (MQ=255)
                                    gCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCa                                                                                                         <  2:384863/150‑1 (MQ=255)
                                    gCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCa                                                                                                         >  1:1837944/1‑150 (MQ=255)
                                     cTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCAc                                                                                                        >  2:518137/1‑150 (MQ=255)
                                     cTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCAc                                                                                                        <  2:1642000/150‑1 (MQ=255)
                                      ttAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGCCGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACt                                                                                                       >  1:1952802/1‑150 (MQ=255)
                                       tAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTa                                                                                                      >  1:164780/1‑150 (MQ=255)
                                         aGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACg                                                                                                    >  1:1219937/1‑150 (MQ=255)
                                         aGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACg                                                                                                    >  1:171969/1‑150 (MQ=255)
                                          taaaGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGt                                                                                                   <  2:249191/149‑1 (MQ=255)
                                           aaaGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTa                                                                                                  <  1:998645/150‑1 (MQ=255)
                                                  catccatGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATa                                                                                           >  2:269462/1‑150 (MQ=255)
                                                    tccatGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAAt                                                                                         <  2:1292661/150‑1 (MQ=255)
                                                    tccatGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAAt                                                                                         >  1:503453/1‑150 (MQ=255)
                                                    tccatGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAAt                                                                                         >  1:1039041/1‑150 (MQ=255)
                                                     ccatGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCAACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATg                                                                                        >  1:1743600/1‑150 (MQ=255)
                                                        tGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTc                                                                                     >  1:430/1‑150 (MQ=255)
                                                         gTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCg                                                                                    <  2:1925285/150‑1 (MQ=255)
                                                           cTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGcc                                                                                  <  2:1108082/150‑1 (MQ=255)
                                                            tACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCg                                                                                 <  1:1542505/150‑1 (MQ=255)
                                                             aCCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCga                                                                                <  1:1395899/150‑1 (MQ=255)
                                                             aCCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCga                                                                                >  2:874891/1‑150 (MQ=255)
                                                              ccGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACCACCAGCATTGCACTACGTAATCCATAATGTTCGCCgag                                                                               >  2:2152408/1‑150 (MQ=255)
                                                              ccGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCgag                                                                               >  2:2161016/1‑150 (MQ=255)
                                                                gCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCa                                                                                                                 >  2:598746/1‑114 (MQ=255)
                                                                gCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCa                                                                                                                 <  1:598746/114‑1 (MQ=255)
                                                                     taatttaatGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGccc                                                                        <  2:1847843/150‑1 (MQ=255)
                                                                     taatttaatGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGccc                                                                        <  1:1027103/150‑1 (MQ=255)
                                                                      aatttaatGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCca                                                                       <  2:974445/150‑1 (MQ=255)
                                                                      aatttaatGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCca                                                                       >  1:897121/1‑150 (MQ=255)
                                                                       atttaatGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCcag                                                                      >  1:938378/1‑150 (MQ=255)
                                                                          taatGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCcagcag                                                                   >  1:1196431/1‑150 (MQ=255)
                                                                            atGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAgcg                                                                 <  1:1920780/150‑1 (MQ=255)
                                                                             tGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAgcgg                                                                <  1:416505/150‑1 (MQ=255)
                                                                             tGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAgcgg                                                                <  2:1944254/150‑1 (MQ=255)
                                                                                tcatcaACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGcggcgg                                                             <  1:336300/150‑1 (MQ=255)
                                                                                  atcaACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGcc                                                           <  2:296046/150‑1 (MQ=255)
                                                                                    caACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCcg                                                         >  1:1032191/1‑150 (MQ=255)
                                                                                    caACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCcg                                                         >  1:988201/1‑150 (MQ=255)
                                                                                      aCTATTCTCCATCGCCGTGTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATATCAAGCGCCAGTACAACCAGCATTGCACTAGGTAATCCATAATGTTCGCCGAGTTAGGCCAGCAGCGGCGGCCcgac                                                       <  1:1704274/150‑1 (MQ=255)
                                                                                      aCTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCcgac                                                       >  2:688426/1‑150 (MQ=255)
                                                                                            ctcCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGaa                                                 <  2:11551/150‑1 (MQ=255)
                                                                                            ctcCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGaa                                                 >  1:1277575/1‑150 (MQ=255)
                                                                                             tcCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGaaa                                                >  2:1076952/1‑150 (MQ=255)
                                                                                              ccATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAg                                               >  2:639342/1‑150 (MQ=255)
                                                                                                   gCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGa                                          <  1:1943393/150‑1 (MQ=255)
                                                                                                        cTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAAcc                                     <  2:1837944/150‑1 (MQ=255)
                                                                                                                    aTCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTAc                         >  2:1299170/1‑150 (MQ=255)
                                                                                                                     tCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTAcc                        <  2:423301/150‑1 (MQ=255)
                                                                                                                      cGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACca                       <  2:1864755/150‑1 (MQ=255)
                                                                                                                             ggCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGAc                <  2:1245809/150‑1 (MQ=255)
                                                                                                                              gCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGAcg               <  2:1219937/150‑1 (MQ=255)
                                                                                                                                 aCGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTGGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCgg            <  2:1598250/150‑1 (MQ=255)
                                                                                                                                 aCGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCgg            >  2:1599430/1‑150 (MQ=255)
                                                                                                                                  ttattctCGCGGCTATTGTCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCggg           >  2:113560/7‑150 (MQ=255)
                                                                                                                                  cGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCggg           >  1:721465/1‑150 (MQ=255)
                                                                                                                                   ggCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGt          >  2:1237416/1‑150 (MQ=255)
                                                                                                                                    gCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGtt         <  1:1852894/150‑1 (MQ=255)
                                                                                                                                    gCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGtt         <  1:1089573/150‑1 (MQ=255)
                                                                                                                                     cTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTg        >  2:2025991/1‑150 (MQ=255)
                                                                                                                                      tttCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTgg       <  2:1798091/150‑1 (MQ=255)
                                                                                                                                        tCGCGGCTATTGTCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGCTAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTg     >  2:458656/1‑150 (MQ=255)
                                                                                                                                        tCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTg     <  2:503540/150‑1 (MQ=255)
                                                                                                                                        tCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTg     <  1:1472828/150‑1 (MQ=255)
                                                                                                                                         cgcgACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTGc    <  2:639105/150‑1 (MQ=255)
                                                                                                                                          gcgACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTGCa   >  1:1492108/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
AGCCATAAACCGCTCACGGTTATAGGTTTTTTGATCGCTTAAGAAAGTACCATCCATGTCTACCGCAATTAATTTAATGCTCATCAACTATTCTCCATCGCCGTCTGCGTTTTGGTATCGGGTTTGGCGACGGCTTTCGCGACAATAGCCGCGAGAATAACCAGCGCCAGTACAACCAGCATTGCACTACGTAATCCATAATGTTCGCCGAGATAGCCCAGCAGCGGCGGCCCGACGAGGAAAGCCAGATAACCGGTCGTTGCTACCACACTGACGCGGGTTGGTGCAT  >  NC_000913/886048‑886336

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