Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,770,646 G→T 13.4% P245Q (CCG→CAG)  yibH ← inner membrane protein YibH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,770,6460GT13.4% 265.7 / 23.1 93P245Q (CCG→CAG) yibHinner membrane protein YibH
Reads supporting (aligned to +/- strand):  ref base G (33/47);  new base T (7/5);  total (41/52)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

CAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTCGCCGTTAACC  >  NC_000913/3770498‑3770795
                                                                                                                                                    |                                                                                                                                                     
cAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCgg                                                                                                                                                      >  1:1208209/1‑150 (MQ=255)
 aGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGc                                                                                                                                                     >  1:1505955/1‑150 (MQ=255)
  gTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCa                                                                                                                                                    <  2:1253477/150‑1 (MQ=255)
    ttGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAg                                                                                                                                                  <  1:352834/150‑1 (MQ=255)
    ttGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAg                                                                                                                                                  >  1:1206311/1‑150 (MQ=255)
    ttGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAg                                                                                                                                                  <  2:22526/150‑1 (MQ=255)
      gCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGc                                                                                                                                                >  1:92331/1‑150 (MQ=255)
       ccGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGGTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCa                                                                                                                                               <  2:1112548/150‑1 (MQ=255)
        cGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAg                                                                                                                                              >  2:967946/1‑150 (MQ=255)
          tGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCt                                                                                                                                            >  2:932387/1‑150 (MQ=255)
             aaCACCTGCCCAGGTAGCGCGTTAATCACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCg                                                                                                                                         <  2:1380820/150‑1 (MQ=255)
               cacCTGCCGAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTa                                                                                                                                       <  1:588197/150‑1 (MQ=255)
               cacCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCTGCCTTGCCGCTGCGtc                                                                                                                                       <  2:803817/150‑2 (MQ=38)
               cacCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTa                                                                                                                                       <  2:494648/150‑1 (MQ=255)
               cacCTGCCCAGGTAGCGCGGTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTa                                                                                                                                       <  2:435279/150‑1 (MQ=255)
                acCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTAt                                                                                                                                      <  2:205689/150‑1 (MQ=255)
                 ccTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCTGCCTTGCCGCTGCGtccg                                                                                                                                     <  2:1004115/150‑4 (MQ=255)
                  cTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCTGCCTTGCCGCTGCGtccgg                                                                                                                                    <  2:40848/150‑5 (MQ=255)
                  cTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCTGCCTTGCCGCTGCGtccgg                                                                                                                                    <  1:1178091/150‑5 (MQ=255)
                  cTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCTGCCTTGCCGCTGCGtccgg                                                                                                                                    <  1:648129/150‑5 (MQ=255)
                    gCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCACACGTAACAGCCAGTTTTGCCGAAATTGGCAGACAATTTGCCGTATTTGCTCCGGCATGAAGACCATCACCGGACGCAGCTGCAATGCAGCAGCGTATGTAc                                                                                                                                  >  1:641242/1‑150 (MQ=255)
                        aGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTgg                                                                                                                              >  2:1616469/1‑150 (MQ=255)
                            aGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGc                                                                                                                                                     <  1:988151/123‑1 (MQ=255)
                            aGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGc                                                                                                                                                     >  2:988151/1‑123 (MQ=255)
                             gcgcGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGa                                                                                                                         <  2:1206311/150‑1 (MQ=255)
                               gcgTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATc                                                                                                                       <  2:51725/150‑1 (MQ=255)
                               gcgTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATc                                                                                                                       <  1:1285296/150‑1 (MQ=255)
                                  ttAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGt                                                                                                                    <  1:707118/150‑1 (MQ=255)
                                   tAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTa                                                                                                                   <  1:1471737/150‑1 (MQ=255)
                                    aaaCACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTAc                                                                                                                  >  1:391897/1‑150 (MQ=255)
                                       caccacTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTg                                                                                                               >  2:360104/1‑150 (MQ=255)
                                       caccacTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTg                                                                                                               <  1:372490/150‑1 (MQ=255)
                                         ccacTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAg                                                                                                             >  2:685938/1‑150 (MQ=255)
                                         ccacTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTCGGCGACAATTTGCAGTTTTTGCTCGGGGATGAAGAACATCACCGGACGCAGCGGCAAGCCAGCTGCGTATCTACCTGGGCGGATCAGTACCTGAg                                                                                                             >  2:155531/1‑150 (MQ=255)
                                            cTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTaa                                                                                                          >  1:1199041/1‑150 (MQ=255)
                                               ccGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGt                                                                                                       >  2:1289862/1‑150 (MQ=255)
                                                 gCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAg                                                                                                     <  2:272245/150‑1 (MQ=255)
                                                  catcatCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGGCCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGc                                                                                                    >  2:64312/1‑150 (MQ=255)
                                                  catcatCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGc                                                                                                    >  1:1130685/1‑150 (MQ=255)
                                                  catcatCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGc                                                                                                    <  2:1470525/150‑1 (MQ=255)
                                                         gCCAGGTTTCAGACGTAACAGCGAGTTTTCCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCt                                                                                             <  1:738580/150‑1 (MQ=255)
                                                             ggTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGgc                                                                                         <  1:907373/150‑1 (MQ=255)
                                                               attaaGACGTAAAATCGATTTTTGACGAAATTGGTAGACAATTTGAAGTTTTTACTAGGGGATGAAGACCATAACAGGAATCAACATCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGgcgc                                                                                       >  2:819077/5‑150 (MQ=255)
                                                                ttCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATAGCTCGgcgcg                                                                                      >  1:1218919/1‑150 (MQ=255)
                                                                  cAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCgg                                                                                          >  1:278259/1‑144 (MQ=255)
                                                                  cAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCgg                                                                                          <  2:278259/144‑1 (MQ=255)
                                                                         aaCAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACa                                                                             >  2:1578397/1‑150 (MQ=255)
                                                                             tcGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGAGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCTCGAATGACAGTCt                                                                         >  1:125145/2‑150 (MQ=255)
                                                                                    ttGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTTTGCTGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAg                                                                  <  2:151525/150‑1 (MQ=255)
                                                                                     tGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAgg                                                                 <  1:274123/150‑1 (MQ=255)
                                                                                      gCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTTGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGt                                                                <  2:700696/150‑1 (MQ=255)
                                                                                      gCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGt                                                                <  2:899962/150‑1 (MQ=255)
                                                                                      gCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGt                                                                >  2:481387/1‑150 (MQ=255)
                                                                                         gAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTtat                                                             >  2:451489/1‑150 (MQ=255)
                                                                                         gAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTtat                                                             >  1:455535/1‑150 (MQ=255)
                                                                                              tGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATAtttt                                                        <  1:344022/150‑1 (MQ=255)
                                                                                              tGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATAtttt                                                        >  1:585400/1‑150 (MQ=255)
                                                                                                ggCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGc                                                      <  2:120407/150‑1 (MQ=255)
                                                                                                 gCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCt                                                     <  2:651005/150‑1 (MQ=255)
                                                                                                 gCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCt                                                     <  2:391897/150‑1 (MQ=255)
                                                                                                 gCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCt                                                     <  2:1199041/150‑1 (MQ=255)
                                                                                                       aTTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTa                                               <  2:139017/150‑1 (MQ=255)
                                                                                                        tttGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTaa                                              <  1:967946/150‑1 (MQ=255)
                                                                                                          tGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGt                                            >  2:1204002/1‑150 (MQ=255)
                                                                                                          tGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGt                                            >  1:610966/1‑150 (MQ=255)
                                                                                                            ccGTTTTTGCTCGGGGATGAAGACCATCACCGTACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCTAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTg                                          <  1:549531/150‑1 (MQ=255)
                                                                                                              gTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTgcg                                        <  1:267171/150‑1 (MQ=255)
                                                                                                                  ttGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTTCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGctct                                    <  1:1222128/150‑1 (MQ=255)
                                                                                                                    gCTCGGGGATGAAGGCCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTa                                  <  1:64312/150‑1 (MQ=255)
                                                                                                                    gCTCGGGGATGAAGACCATCACCGGACGGAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTa                                  <  1:155531/150‑1 (MQ=255)
                                                                                                                    gCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTa                                  <  2:1218919/150‑1 (MQ=255)
                                                                                                                    gCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCGGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGGCTGCTCAAGGTTATATTTTGCTTCAGTAATTTGCGCTCTTa                                  <  2:520298/150‑1 (MQ=255)
                                                                                                                       cGGGGATGAAGACCATCACCGGACGCAGCGTCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGc                               <  1:765660/150‑1 (MQ=255)
                                                                                                                        ggggATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCt                              <  2:1208209/150‑1 (MQ=255)
                                                                                                                          ggATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTca                            >  1:485602/1‑150 (MQ=255)
                                                                                                                               aaGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATc                       >  2:570259/1‑150 (MQ=255)
                                                                                                                                aGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCt                      <  2:1099737/150‑1 (MQ=255)
                                                                                                                                 gACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTg                     >  1:592784/1‑150 (MQ=255)
                                                                                                                                  aCCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGc                    <  2:455535/150‑1 (MQ=255)
                                                                                                                                  aCCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGc                    <  2:1607097/150‑1 (MQ=255)
                                                                                                                                   ccATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCg                   >  2:1591853/1‑150 (MQ=255)
                                                                                                                                   ccATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCg                   >  1:373931/1‑150 (MQ=255)
                                                                                                                                   ccATCACCAGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCg                   >  1:631425/1‑150 (MQ=255)
                                                                                                                                     aTCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGAt                 >  2:580921/1‑150 (MQ=255)
                                                                                                                                       cACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTg               <  2:830094/150‑1 (MQ=255)
                                                                                                                                       cACCGGACGCAGCGGCAAGGCAGCTGCATATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTg               >  1:104944/1‑150 (MQ=255)
                                                                                                                                        tacataCGCAGCTGCCTTGCCGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGc              >  2:134331/6‑150 (MQ=38)
                                                                                                                                        tacataCGCAGCTGCCTTGCCGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGc              >  1:973894/6‑150 (MQ=38)
                                                                                                                                        aCCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGc              <  2:799722/150‑1 (MQ=255)
                                                                                                                                          cataCGCAGCTGCCTTGCCGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTc            >  1:1477687/4‑150 (MQ=38)
                                                                                                                                           ataCGCAGCTGCCTTGCCGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTcc           >  2:133877/3‑149 (MQ=38)
                                                                                                                                           ataCGCAGCTGCCTTGCCGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTCg           >  1:1595517/3‑150 (MQ=38)
                                                                                                                                           ataCGCAGCTGCCTTGCCGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTCg           >  1:614886/3‑150 (MQ=38)
                                                                                                                                                 aGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTCGCCGTTa     <  1:1616469/150‑1 (MQ=255)
                                                                                                                                                    ggCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTCGCCGTTAAcc  <  1:550376/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTCGCCGTTAACC  >  NC_000913/3770498‑3770795

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: