Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,831,757 C→T 10.2% A434V (GCC→GTC)  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,7570CT10.2% 311.2 / 17.7 108A434V (GCC→GTC) yicHputative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base C (49/48);  new base T (7/4);  total (56/52)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

Only 100 of 118 total aligned reads displayed.

GAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAAC  >  NC_000913/3831608‑3831906
                                                                                                                                                     |                                                                                                                                                     
taaGAGATTTAGCGCCAGCCGCAATCTGCTCATTGATATCGACCCGGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgc                                                                                                                                                       <  1:1279793/149‑1 (MQ=255)
gAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgc                                                                                                                                                       >  1:548392/1‑150 (MQ=255)
 aaGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgcc                                                                                                                                                      >  2:72800/1‑150 (MQ=255)
  agagATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccg                                                                                                                                                     <  1:1113038/150‑1 (MQ=255)
     gaTTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccgccg                                                                                                                                                  >  1:1713264/1‑150 (MQ=255)
      aTTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccgccgc                                                                                                                                                 >  1:1497768/1‑150 (MQ=255)
         tAGCGCCAGCCGCAATCTGATCAATCATAGCGCCCCTGACTTCCCTGCGCAGCTCACCGCGCTCGATGGCTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCcgccgccgc                                                                                                                                              <  1:281237/150‑1 (MQ=255)
          aGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCCCCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccgtcgccagt                                                                                                                                             >  2:64668/1‑147 (MQ=255)
             gcCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACa                                                                                                                                          >  2:752639/1‑150 (MQ=255)
             gcCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACa                                                                                                                                          >  2:1440038/1‑150 (MQ=255)
              ccAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGGGGCAGCTCACCGCGCGCGATGGTTCCGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACAt                                                                                                                                         <  1:180534/150‑1 (MQ=255)
                aGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACTCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGCATGCAGCCGCCGCGagaatc                                                                                                                                       >  2:277221/1‑145 (MQ=255)
                 gCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCa                                                                                                                                      <  2:88225/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)
                     cAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGAGAATCTGAATGCCGCCGCCGCGACCTTCAATCg                                                                                                                                  >  2:475055/1‑150 (MQ=255)
                      aaTCTGATCATTGATCTCGACCCTGACTTCCCGTGGCAGCTCACCGCGCGCGATGGTTACGGTGCCAACCTGACTCTGGTTACCGAGCATAAATGGCGCGTCTGGAGTGTCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCgt                                                                                                                                 <  1:1255880/150‑1 (MQ=255)
                      aaTCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCGCGGCGAATCTGAATGTCGCGGCGGCGGCATTCagattc                                                                                                                                 <  2:468289/150‑6 (MQ=255)
                         cTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGTTCATAAATGGGGCGTCTGGAGTGGCTCGGCGATTCTGAATGCCGCCGCCGCGACATTCAATCGTGtt                                                                                                                              >  1:1341054/1‑150 (MQ=255)
                          tGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTTGCTCGGCGAATCTGAATGCCGCCGCCGCGCCATTCAATCGTGTTg                                                                                                                             >  2:94833/1‑150 (MQ=255)
                          tGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTg                                                                                                                             <  2:1153455/150‑1 (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)
                                            cTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGACAACCTGACGCTGGTTAACGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTAGTCGCACGTCGc                                                                                                           >  1:831240/1‑150 (MQ=255)
                                             tGACTTCGCGTGGCAGCTCACCGCGCTCGTTGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAAGGGGGCGTCTGGAGTGGCGCGGCGAATCTGAATGCCGCCGCCGCGAGATTCAATCGTGTTGATGTTCGTCGCCCGTCGCt                                                                                                          <  1:282042/150‑1 (MQ=255)
                                                 ttCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGgcg                                                                                                      >  2:1682703/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)
                                                               cACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcag                                                                                        >  2:1478217/1‑150 (MQ=255)
                                                                  cgcgCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcagcag                                                                                     <  2:613460/150‑1 (MQ=255)
                                                                  cgcgCTCGATGGGTAGGGTGCCAACCTGACGCTGGGTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCAGCGACATTCAATCGTGTTGATGGTCGTCGCCCGTCGCTGGCGGTGACCGCCAAcagcag                                                                                     <  2:413294/150‑1 (MQ=255)
                                                                      cTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGt                                                                                                              >  1:693305/1‑121 (MQ=255)
                                                                         gATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCCACAGCAGCACGGTg                                                                              >  2:1552831/1‑150 (MQ=255)
                                                                           tGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGaa                                                                            >  1:191073/1‑150 (MQ=255)
                                                                                aCGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCa                                                                       >  1:746835/1‑150 (MQ=255)
                                                                                 cGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAg                                                                      <  2:575112/150‑1 (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)
                                                                                            tGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTg                                                           >  1:411480/1‑150 (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)
                                                                                                   ggTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTAc                                                    <  1:974037/150‑1 (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)
                                                                                                          gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGATTATCAGCGAACTGAGTGCATTTACTGaaaaa                                             >  1:212218/1‑150 (MQ=255)
                                                                                                           aTCATAAATGGTGCGTCTGGAGTGGCTCGGCTAATCTGAATGCCGCCGCCTCGACATTCAATCGTGTTGATGTTCGTAGCCCGTCGCTGTCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAg                                            <  1:894939/150‑1 (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)
                                                                                                                                 tGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAGCAGCAGCACGGTGAATATCAGCGAACTGAGTGCTTTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAg                      <  1:76705/150‑1 (MQ=255)
                                                                                                                                 gggCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAg                      <  1:191101/149‑1 (MQ=255)
                                                                                                                                  ggCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAgt                     >  1:1651206/1‑150 (MQ=255)
                                                                                                                                  ggCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAgt                     >  2:240256/1‑150 (MQ=255)
                                                                                                                                    cTCGTCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAgtgt                   <  1:1129042/150‑1 (MQ=255)
                                                                                                                                    cTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAgtgt                   >  2:747225/1‑150 (MQ=255)
                                                                                                                                      cGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGttt                 >  2:722983/1‑150 (MQ=255)
                                                                                                                                        gCGAATCTGAATGCCGCCGCCTCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTca               <  1:1216915/150‑1 (MQ=255)
                                                                                                                                          acgattGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGTAAGCCACCGCCAGTGTTTcaca             >  2:405291/6‑150 (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)
                                                                                                                                                 aaTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCca      >  2:66658/1‑150 (MQ=255)
                                                                                                                                                  aTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCcac     <  2:538682/150‑1 (MQ=255)
                                                                                                                                                  aTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCcac     >  2:1113999/1‑150 (MQ=255)
                                                                                                                                                     ccgccgccgCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACACACGCCACAAc  >  1:1722888/1‑150 (MQ=255)
                                                                                                                                                     ccgccgccgCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAAc  >  1:1347097/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAAC  >  NC_000913/3831608‑3831906

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: