Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,831,768 A→G 8.9% T438A (ACA→GCA)  yicH → putative inner membrane‑anchored periplasmic AsmA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,831,7680AG 8.9% 339.3 / 19.3 124T438A (ACA→GCA) yicHputative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base A (54/57);  new base G (7/4);  total (63/61)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

Only 100 of 125 total aligned reads displayed.

GCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACA  >  NC_000913/3831621‑3831916
                                                                                                                                                   |                                                                                                                                                    
gcCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACa                                                                                                                                                    >  2:1440038/1‑150 (MQ=255)
gcCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACa                                                                                                                                                    >  2:752639/1‑150 (MQ=255)
gcCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGATACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACa                                                                                                                                                    <  2:1385609/150‑1 (MQ=255)
 ccAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGGGGCAGCTCACCGCGCGCGATGGTTCCGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACAt                                                                                                                                                   <  1:180534/150‑1 (MQ=255)
   aGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACTCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGCATGCAGCCGCCGCGagaatc                                                                                                                                                 >  2:277221/1‑145 (MQ=255)
    gCCGCAATCTGAGCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCGCGATGGCTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGCGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCa                                                                                                                                                <  1:836217/150‑1 (MQ=255)
    cccGCAATCTGATCATTGATCTCGACCCTGACTTCCCGGGGCAGCTCACCGCGCTCGAGGGGTACGGTGCCAACCTGACGCTGGCTACCGATCATAAATGGGGCGGCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCa                                                                                                                                                <  1:800424/149‑1 (MQ=255)
     ccGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGCCATTCaa                                                                                                                                               >  1:695721/1‑150 (MQ=255)
     ccGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCaa                                                                                                                                               >  2:141381/1‑150 (MQ=255)
     ccGCAATCTGATCATTGATATCGACCCTGACTTCCCCTGGCAGCTCACCGCGCTCGATGGGTACGGTGCCAACCTGCCGCTGGTTTCCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCag                                                                                                                                               <  1:315296/150‑2 (MQ=255)
       gCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATc                                                                                                                                             <  1:387650/150‑1 (MQ=255)
       gCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATc                                                                                                                                             >  2:1505806/1‑150 (MQ=255)
        cAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGAGAATCTGAATGCCGCCGCCGCGACCTTCAATCg                                                                                                                                            >  2:475055/1‑150 (MQ=255)
            cTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGTTCATAAATGGGGCGTCTGGAGTGGCTCGGCGATTCTGAATGCCGCCGCCGCGACATTCAATCGTGtt                                                                                                                                        >  1:1341054/1‑150 (MQ=255)
             tGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTTGCTCGGCGAATCTGAATGCCGCCGCCGCGCCATTCAATCGTGTTg                                                                                                                                       >  2:94833/1‑150 (MQ=255)
             tGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTg                                                                                                                                       >  2:90572/1‑150 (MQ=255)
             tGATCATTGATATCGACCCTGACATCCCGGTGCAGCTCCCCGCGCTCGATGGTTACGGTGCCAACCTGACGCTTGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTg                                                                                                                                       <  2:1210085/150‑1 (MQ=255)
                     gATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTtcgt                                                                                                                               >  2:755298/1‑150 (MQ=255)
                      atatCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATACATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACCTTCAATCGTGTTGATGTtcgtc                                                                                                                              >  1:263240/1‑150 (MQ=255)
                      atatCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTtcgtc                                                                                                                              <  1:768159/150‑1 (MQ=255)
                          cGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGccc                                                                                                                          <  1:1702978/150‑1 (MQ=255)
                           gACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCg                                                                                                                         <  1:769517/150‑1 (MQ=255)
                             cccTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGGTACCGATCATAAATGGGGCGGCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTc                                                                                                                       >  2:463713/1‑150 (MQ=255)
                               cTGCCTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGGGCCAACCTGACGCGGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCGCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGc                                                                                                                     <  2:1703617/150‑1 (MQ=255)
                               cTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGACAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGc                                                                                                                     >  2:76705/1‑150 (MQ=255)
                               cTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGACAACCTGACGCTGGTTAACGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTAGTCGCACGTCGc                                                                                                                     >  1:831240/1‑150 (MQ=255)
                                     tCCCGTGGCAGCGCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCGCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGgcgc                                                                                                               <  1:1235134/150‑1 (MQ=255)
                                      cccGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCt                                                                                                              >  2:499659/1‑150 (MQ=255)
                                       ccGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTg                                                                                                             >  2:834402/1‑150 (MQ=255)
                                        cGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGa                                                                                                            >  2:943791/1‑150 (MQ=255)
                                          tGGCAGCTCACCGCGCTGGATGGTTTCGGTGCCAACCTGACGCTGGTTAGCGATCATCTGTGGGGCGTCTGGAGGGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTGGCTGTTCGTCGCCCGTCGCGGGCGCGGAcc                                                                                                          <  2:1163236/150‑1 (MQ=255)
                                          tGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAcc                                                                                                          <  1:850632/150‑1 (MQ=255)
                                              aGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCa                                                                                                      <  1:33355/150‑1 (MQ=255)
                                                cTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAc                                                                                                    <  1:1005558/150‑1 (MQ=255)
                                                cTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAc                                                                                                    <  1:1539585/150‑1 (MQ=255)
                                                     cgcgCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcagcag                                                                                               <  2:613460/150‑1 (MQ=255)
                                                     cgcgCTCGATGGGTAGGGTGCCAACCTGACGCTGGGTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCAGCGACATTCAATCGTGTTGATGGTCGTCGCCCGTCGCTGGCGGTGACCGCCAAcagcag                                                                                               <  2:413294/150‑1 (MQ=255)
                                                         cTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGt                                                                                                                        <  2:693305/121‑1 (MQ=255)
                                                            gATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCCACAGCAGCACGGTg                                                                                        >  2:1552831/1‑150 (MQ=255)
                                                              tGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGaa                                                                                      >  1:191073/1‑150 (MQ=255)
                                                               ggTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAat                                                                                     <  2:532140/150‑1 (MQ=255)
                                                                  tACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACAGTGAATATc                                                                                  >  1:1254536/1‑150 (MQ=255)
                                                                   aCGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCa                                                                                 >  1:746835/1‑150 (MQ=255)
                                                                    cGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAg                                                                                <  2:575112/150‑1 (MQ=255)
                                                                      gTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGTGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAAACGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACTGCAGCACGGTGAATATCAGCg                                                                              >  2:1086235/1‑150 (MQ=255)
                                                                        gCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGaa                                                                            >  2:20509/1‑150 (MQ=255)
                                                                         ccAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAAc                                                                           <  2:481750/150‑1 (MQ=255)
                                                                          cAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACt                                                                          <  1:1440038/150‑1 (MQ=255)
                                                                           aaCCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTg                                                                         <  1:141381/150‑1 (MQ=255)
                                                                              cTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGt                                                                      <  2:1713264/150‑1 (MQ=255)
                                                                              cTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGt                                                                      <  1:976086/150‑1 (MQ=255)
                                                                               tGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTg                                                                     <  1:706934/150‑1 (MQ=255)
                                                                                 aCGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCa                                                                   >  1:48134/1‑150 (MQ=255)
                                                                                      ggTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAAc                                                                           <  2:493908/137‑1 (MQ=255)
                                                                                      ggTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAAc                                                                           >  1:493908/1‑137 (MQ=255)
                                                                                      ggTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTAc                                                              >  1:1674668/1‑150 (MQ=255)
                                                                                       gTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACt                                                             <  2:1338406/150‑1 (MQ=255)
                                                                                        ttACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTg                                                            <  2:1062077/150‑1 (MQ=255)
                                                                                          aCCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGaa                                                          >  2:705644/1‑150 (MQ=255)
                                                                                            cGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTg                                                                         >  2:1483392/1‑133 (MQ=255)
                                                                                            cGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTg                                                                         <  1:1483392/133‑1 (MQ=255)
                                                                                             gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGATTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGATCTGAGTGCATTTACTGaaaaa                                                       >  1:206619/1‑150 (MQ=255)
                                                                                             gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGaaaaa                                                       <  2:441598/150‑1 (MQ=255)
                                                                                              aTCATAAATGGTGCGTCTGGAGTGGCTCGGCTAATCTGAATGCCGCCGCCTCGACATTCAATCGTGTTGATGTTCGTAGCCCGTCGCTGTCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAg                                                      <  1:894939/150‑1 (MQ=255)
                                                                                                  tAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCACGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATTTCAGCGAACTGAGTGCATTTACTGAAAACGGCAt                                                  >  1:249006/1‑150 (MQ=255)
                                                                                                      tGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCTATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAGCGCCAACAGCAGCACGGGGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTg                                              <  1:1148052/150‑1 (MQ=255)
                                                                                                        gggAGTCTGGAGTGGCTCTGCGAATCTGAATGCCGCCGCCGCGACTTTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCTCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGa                                            <  1:168225/150‑1 (MQ=255)
                                                                                                          gCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAg                                          <  1:943791/150‑1 (MQ=255)
                                                                                                                  aGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGATGCCCGCCAACAGCAGAACGGTGAATATCAGCGAACTGAGTGCATTTACAGAAAAACGCATTCTGGAAGCCACCGcc                                  >  2:216483/1‑150 (MQ=255)
                                                                                                                    gggCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAg                                <  1:191101/149‑1 (MQ=255)
                                                                                                                     ggCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAgt                               >  2:240256/1‑150 (MQ=255)
                                                                                                                       cTCGTCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAgtgt                             <  1:1129042/150‑1 (MQ=255)
                                                                                                                         cGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGttt                           >  2:722983/1‑150 (MQ=255)
                                                                                                                           gCGAATCTGAATGCCGCCGCCTCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTca                         <  1:1216915/150‑1 (MQ=255)
                                                                                                                               gattGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACaac                     >  1:439703/4‑149 (MQ=255)
                                                                                                                                tCTGACTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAAc                    <  2:1497768/150‑1 (MQ=255)
                                                                                                                                tCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAAc                    >  1:380024/1‑150 (MQ=255)
                                                                                                                                tCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAAc                    >  2:1220013/1‑150 (MQ=255)
                                                                                                                                 ttGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACg                   >  2:176362/2‑150 (MQ=255)
                                                                                                                                 cTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACg                   >  1:1344800/1‑150 (MQ=255)
                                                                                                                                  tGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGc                  >  1:1355287/1‑150 (MQ=255)
                                                                                                                                  tGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGc                  >  1:64668/1‑150 (MQ=255)
                                                                                                                                   gAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGcc                 >  2:170761/1‑150 (MQ=255)
                                                                                                                                   gAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGcc                 >  1:595937/1‑150 (MQ=255)
                                                                                                                                     aTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCcac               >  2:1113999/1‑150 (MQ=255)
                                                                                                                                     aTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCcac               <  2:538682/150‑1 (MQ=255)
                                                                                                                                     aTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCcac               <  1:755298/150‑1 (MQ=255)
                                                                                                                                        ccgccgccgCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACACACGCCACAAc            >  1:1722888/1‑150 (MQ=255)
                                                                                                                                        ccgccgccgCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAAc            >  1:1347097/1‑150 (MQ=255)
                                                                                                                                          gccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGt          >  2:460719/1‑150 (MQ=255)
                                                                                                                                          gccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGt          >  2:1215587/1‑150 (MQ=255)
                                                                                                                                          gccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGt          <  1:1201290/150‑1 (MQ=255)
                                                                                                                                            cgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCAACCGCCAGTGTTTCACAAACGCCACAACGTCa        >  2:94860/1‑150 (MQ=255)
                                                                                                                                             gccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAg       <  1:100047/150‑1 (MQ=255)
                                                                                                                                             gccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAg       <  1:569736/150‑1 (MQ=255)
                                                                                                                                               cgcgACATTCAATCGTGTTGATGTTCGTCGCCCGTCTCTGGCGCTGACCGCCAACAGCATCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGac     <  1:475055/150‑1 (MQ=255)
                                                                                                                                                gcgACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCCTTTACTGATAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACACACGCCACAACGTCAGacc    >  1:272713/1‑149 (MQ=255)
                                                                                                                                                gcgACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGaca    >  1:577027/1‑150 (MQ=255)
                                                                                                                                                 cgACATTCTAACGTGTTGATGTTAGTCGCCCGTCGCTGGCGCTGACATCCAACAGCAGCACGGTGAATATCAGCGATCTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGacac   <  1:277221/150‑1 (MQ=255)
                                                                                                                                                 cgACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGacac   <  2:746835/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACA  >  NC_000913/3831621‑3831916

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