Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,831,761 C→G 10.2% A435A (GCC→GCG 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,7610CG10.2% 315.1 / 22.5 107A435A (GCC→GCGyicHputative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base C (49/47);  new base G (7/4);  total (56/51)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

Only 100 of 119 total aligned reads displayed.

GATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCA  >  NC_000913/3831613‑3831910
                                                                                                                                                    |                                                                                                                                                     
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:1440038/1‑150 (MQ=255)
        gcCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACa                                                                                                                                              >  2:752639/1‑150 (MQ=255)
        gcCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGATACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACa                                                                                                                                              <  2:1385609/150‑1 (MQ=255)
           aGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACTCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGCATGCAGCCGCCGCGagaatc                                                                                                                                           >  2:277221/1‑145 (MQ=255)
            gCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCa                                                                                                                                          <  2:88225/150‑1 (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                                                                                                                                       >  2:1505806/1‑150 (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)
                              atatCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTtcgtc                                                                                                                        <  1:768159/150‑1 (MQ=255)
                                  cGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGccc                                                                                                                    <  1:1702978/150‑1 (MQ=255)
                                   gACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCg                                                                                                                   <  1:769517/150‑1 (MQ=255)
                                       cTGCCTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGGGCCAACCTGACGCGGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCGCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGc                                                                                                               <  2:1703617/150‑1 (MQ=255)
                                       cTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGACAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGc                                                                                                               >  2:76705/1‑150 (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)
                                                cGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGa                                                                                                      >  2:943791/1‑150 (MQ=255)
                                                  tGGCAGCTCACCGCGCTGGATGGTTTCGGTGCCAACCTGACGCTGGTTAGCGATCATCTGTGGGGCGTCTGGAGGGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTGGCTGTTCGTCGCCCGTCGCGGGCGCGGAcc                                                                                                    <  2:1163236/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)
                                                                 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)
                                                                           aCGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCa                                                                           >  1:746835/1‑150 (MQ=255)
                                                                            cGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAg                                                                          <  2:575112/150‑1 (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                                                                <  1:976086/150‑1 (MQ=255)
                                                                                      cTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGt                                                                <  2:1713264/150‑1 (MQ=255)
                                                                                       tGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTg                                                               >  1:411480/1‑150 (MQ=255)
                                                                                       tGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTg                                                               <  1:706934/150‑1 (MQ=255)
                                                                                         aCGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCa                                                             >  1:48134/1‑150 (MQ=255)
                                                                                              ggTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAAc                                                                     >  1:493908/1‑137 (MQ=255)
                                                                                              ggTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAAc                                                                     <  2:493908/137‑1 (MQ=255)
                                                                                              ggTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTAc                                                        <  1:974037/150‑1 (MQ=255)
                                                                                              ggTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTAc                                                        >  1:1674668/1‑150 (MQ=255)
                                                                                               gTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACt                                                       <  2:1338406/150‑1 (MQ=255)
                                                                                                ttACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTg                                                      <  2:1062077/150‑1 (MQ=255)
                                                                                                    cGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTg                                                                   <  1:1483392/133‑1 (MQ=255)
                                                                                                     gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGATTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGATCTGAGTGCATTTACTGaaaaa                                                 >  1:206619/1‑150 (MQ=255)
                                                                                                     gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGATTATCAGCGAACTGAGTGCATTTACTGaaaaa                                                 >  1:212218/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)
                                                                                                                               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)
                                                                                                                                         ttGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACg             >  2:176362/2‑150 (MQ=255)
                                                                                                                                         ttGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCCCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACAGTGTGCATTTAGTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGGTTCCCAAAcc             >  1:556129/2‑149 (MQ=255)
                                                                                                                                         cTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACg             >  1:1344800/1‑150 (MQ=255)
                                                                                                                                          tGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGc            >  1:1355287/1‑150 (MQ=255)
                                                                                                                                          tGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGc            >  1:64668/1‑150 (MQ=255)
                                                                                                                                           gAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCCCCGCCAGTGTTTCACAAACGcc           >  1:1733347/1‑150 (MQ=255)
                                                                                                                                            aaTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCca          >  2:66658/1‑150 (MQ=255)
                                                                                                                                             aTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCcac         >  2:1113999/1‑150 (MQ=255)
                                                                                                                                             aTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCcac         <  1:755298/150‑1 (MQ=255)
                                                                                                                                                ccgccgccgCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACACACGCCACAAc      >  1:1722888/1‑150 (MQ=255)
                                                                                                                                                ccgccgccgCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAAc      >  1:1347097/1‑150 (MQ=255)
                                                                                                                                                  gccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACgg    >  2:1241894/1‑149 (MQ=255)
                                                                                                                                                  gccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGt    <  1:1201290/150‑1 (MQ=255)
                                                                                                                                                  gccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGt    >  2:460719/1‑150 (MQ=255)
                                                                                                                                                  gccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGt    >  2:1215587/1‑150 (MQ=255)
                                                                                                                                                    cgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCa  <  2:191073/150‑1 (MQ=255)
                                                                                                                                                    cgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCAACCGCCAGTGTTTCACAAACGCCACAACGTCa  >  2:94860/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCA  >  NC_000913/3831613‑3831910

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: