Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,111,823 C→T 5.9% pseudogene (166/861 nt) yghE ← putative type II secretion system L‑type protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,111,8230CT 5.9% 473.8 / 12.0 151pseudogene (166/861 nt)yghEputative type II secretion system L‑type protein
Reads supporting (aligned to +/- strand):  ref base C (72/70);  new base T (6/3);  total (78/73)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.57e-01

Only 100 of 157 total aligned reads displayed.

GGCGGGCAGAGCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATAT  >  NC_000913/3111674‑3111972
                                                                                                                                                     |                                                                                                                                                     
ggCGGGCAGAGCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGccc                                                                                                                                                       <  2:857105/150‑1 (MQ=255)
ggCGGGCAGAGCTGTAACGGTCAGGCGTTCAGCGTCTGCGCAAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGccc                                                                                                                                                       >  2:645609/1‑150 (MQ=255)
  cGGGCAGAGCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCgg                                                                                                                                                     <  2:1097061/150‑1 (MQ=255)
  cGGGCAGAGCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCgg                                                                                                                                                     <  2:539285/150‑1 (MQ=255)
  cGGGCAGAGCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCgg                                                                                                                                                     <  2:151740/150‑1 (MQ=255)
      cAGAGCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCt                                                                                                                                                 >  2:188629/1‑150 (MQ=255)
      cAGAGCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCt                                                                                                                                                 <  2:1971120/150‑1 (MQ=255)
      cAGAGCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCt                                                                                                                                                 >  2:2124080/1‑150 (MQ=255)
        gagCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTcc                                                                                                                                               >  1:2244439/1‑150 (MQ=255)
         agCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCGGGCTGGCGCtggc                                                                                                                                              <  1:2206698/150‑4 (MQ=255)
         agCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCGGGCTGGCGCtggc                                                                                                                                              <  1:2176354/150‑4 (MQ=255)
         agCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCa                                                                                                                                              <  1:1494453/150‑1 (MQ=255)
            tGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATcc                                                                                                                                           <  1:637879/150‑1 (MQ=255)
             gTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCATCCCGGCGCTCCATccc                                                                                                                                          <  2:1934857/150‑1 (MQ=255)
                      aGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCCGCGCCAGCCCGGCGCTCCATCCCCGGCATTCa                                                                                                                                 <  2:543395/150‑1 (MQ=255)
                      aGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCa                                                                                                                                 <  2:2405282/150‑1 (MQ=255)
                       ggCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGTCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCAt                                                                                                                                >  2:490483/1‑150 (MQ=255)
                       ggCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACGGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCAt                                                                                                                                >  1:1037891/1‑150 (MQ=255)
                        gCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATc                                                                                                                               <  2:2238504/150‑1 (MQ=255)
                         cGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCg                                                                                                                              >  1:1241363/1‑150 (MQ=255)
                         cGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGGGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCg                                                                                                                              <  2:1173024/150‑1 (MQ=255)
                          nttCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTANANNCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGa                                                                                                                             >  2:742740/2‑150 (MQ=255)
                          gTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGa                                                                                                                             >  1:1553351/1‑150 (MQ=255)
                          gTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGa                                                                                                                             >  1:317584/1‑150 (MQ=255)
                          gTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGa                                                                                                                             >  1:160214/1‑150 (MQ=255)
                             cAGCGTCTGCGCCAGCCCATGCTGATAGATTTTCCGGTGCAATGACGACTAAGCTTAGCCAAAGAGGTAAATCCTGCTCATTGAGATGGCCGCGGAACATGACCCGAGCCAGCGCCCGCCCCCCCCTCCATCCCCGCCATTCATCGATGc                                                                                                                          >  2:960083/1‑150 (MQ=255)
                              aGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCAACCCCGGCATTCATCGATGCg                                                                                                                         >  2:1358982/1‑150 (MQ=255)
                                cGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGc                                                                                                                       <  2:2244439/150‑1 (MQ=255)
                                  tCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAg                                                                                                                     <  1:69275/150‑1 (MQ=255)
                                   cTGCGGCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGa                                                                                                                    <  2:1153918/150‑1 (MQ=255)
                                   cTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCCTCCCCGGCATTCATCGATGCGGCAGa                                                                                                                    >  1:1626478/1‑150 (MQ=255)
                                   cTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGa                                                                                                                    >  2:1111488/1‑150 (MQ=255)
                                        ccAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCTCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAAcc                                                                                                               <  1:1192221/150‑1 (MQ=255)
                                            cccATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGtt                                                                                                           <  1:1710577/150‑1 (MQ=255)
                                             ccATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTg                                                                                                          >  1:1001574/1‑150 (MQ=255)
                                              cATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTgc                                                                                                         <  1:886174/150‑1 (MQ=255)
                                               aTGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGGGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTgcg                                                                                                        <  1:1663490/150‑1 (MQ=255)
                                                    gATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCg                                                                                                   <  1:772678/150‑1 (MQ=255)
                                                      tAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCa                                                                                                 <  1:588039/150‑1 (MQ=255)
                                                      tAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCa                                                                                                 >  2:1495502/1‑150 (MQ=255)
                                                         cTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGc                                                                                              <  1:2036361/150‑1 (MQ=255)
                                                         cTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGc                                                                                              >  1:979249/1‑150 (MQ=255)
                                                             tCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATa                                                                                          >  1:2005202/1‑150 (MQ=255)
                                                              ccGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGc                                                                                                                                      <  1:1713991/105‑1 (MQ=255)
                                                              ccGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGc                                                                                                                                      >  2:1713991/1‑105 (MQ=255)
                                                              ccGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATAGCAAGCCATaa                                                                                         >  2:201368/1‑150 (MQ=255)
                                                                ggTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAAc                                                                                       <  1:1205621/150‑1 (MQ=255)
                                                                   gCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATAGCAAGCCATAAACCca                                                                                    >  1:617084/1‑150 (MQ=255)
                                                                      aTGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACaa                                                                                 >  2:483548/1‑150 (MQ=255)
                                                                       tGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAg                                                                                <  1:231874/150‑1 (MQ=255)
                                                                       tGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAg                                                                                <  1:576357/150‑1 (MQ=255)
                                                                           gACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAAt                                                                            <  1:869672/150‑1 (MQ=255)
                                                                              tACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGa                                                                         <  2:32301/150‑1 (MQ=255)
                                                                              tACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGa                                                                         <  1:2337570/150‑1 (MQ=255)
                                                                               aCGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGAc                                                                        >  1:1933387/1‑150 (MQ=255)
                                                                                  gTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATc                                                                     >  1:1667543/1‑150 (MQ=255)
                                                                                  gTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATc                                                                     >  1:1891329/1‑150 (MQ=255)
                                                                                   ttAGCGAAAGCGTTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATcc                                                                    <  2:583461/150‑1 (MQ=255)
                                                                                   ttAGCGAAAGCGGTAAATCCTGCGCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATcc                                                                    <  2:1807797/150‑1 (MQ=255)
                                                                                     aGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCAACCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAg                                                                  <  1:1358982/150‑1 (MQ=255)
                                                                                      gCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGt                                                                 <  2:1002147/150‑1 (MQ=255)
                                                                                          aaGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCgg                                                             >  1:687547/1‑150 (MQ=255)
                                                                                           aGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGgc                                                            >  2:584938/1‑150 (MQ=255)
                                                                                           aGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGgc                                                            >  1:336925/1‑150 (MQ=255)
                                                                                             cGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGgcgc                                                          <  1:112213/150‑1 (MQ=255)
                                                                                             cGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGgcgc                                                          <  2:1444971/150‑1 (MQ=255)
                                                                                              ggTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGcgcg                                                         >  1:452264/1‑150 (MQ=255)
                                                                                                     ccTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAg                                                  <  1:1071816/150‑1 (MQ=255)
                                                                                                       tGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAgcg                                                >  2:1397245/1‑150 (MQ=255)
                                                                                                         cTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGcgcg                                              >  1:2394927/1‑150 (MQ=255)
                                                                                                          tCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCgg                                             >  2:277867/1‑150 (MQ=255)
                                                                                                           cATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCggg                                            >  1:928576/1‑150 (MQ=255)
                                                                                                           cATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCggg                                            <  2:1335569/150‑1 (MQ=255)
                                                                                                           cATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCggg                                            <  1:2042992/150‑1 (MQ=255)
                                                                                                             ttGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCggggg                                          <  2:698756/150‑3 (MQ=255)
                                                                                                             ttGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTg                                          >  1:928034/1‑150 (MQ=255)
                                                                                                              tGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGa                                         >  2:442364/1‑150 (MQ=255)
                                                                                                               gAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGAt                                        >  2:1155292/1‑150 (MQ=255)
                                                                                                                  cTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGcc                                     >  2:2097159/1‑150 (MQ=255)
                                                                                                                   tGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCg                                    >  2:1404575/1‑150 (MQ=255)
                                                                                                                         cgGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTcc                              >  2:1829554/1‑150 (MQ=255)
                                                                                                                         cgGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTcc                              <  1:104521/150‑1 (MQ=255)
                                                                                                                          ggAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCt                             <  2:419790/150‑1 (MQ=255)
                                                                                                                          ggAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCACAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCt                             >  1:812569/1‑150 (MQ=255)
                                                                                                                            aaCATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGc                           <  1:1362865/150‑1 (MQ=255)
                                                                                                                            aaCATGAACGGAGCAAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGAGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATACAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGc                           >  2:1695892/1‑150 (MQ=255)
                                                                                                                             aCATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCa                          <  1:483548/150‑1 (MQ=255)
                                                                                                                               aTGACCGGAGCCAGCTCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGTAATGGACATCCAGTCCGGCGCTACGGTAGCGCTGGTGATGCCGTTTTCCTGCAAc                        <  2:359395/150‑1 (MQ=255)
                                                                                                                                tGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAAcc                       >  1:1452133/1‑150 (MQ=255)
                                                                                                                                 gACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCa                      >  2:886102/1‑150 (MQ=255)
                                                                                                                                  aCCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCAc                     >  2:814200/1‑150 (MQ=255)
                                                                                                                                     ggAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTgc                  >  2:254455/1‑150 (MQ=255)
                                                                                                                                      gatggagcgccgggCTGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTgcg                 >  1:1009843/6‑150 (MQ=255)
                                                                                                                                        tggagcgccgggCTGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCCCAAGGAATGGCCATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCa               >  2:481519/4‑150 (MQ=255)
                                                                                                                                        tggagcgccgggCTGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCa               >  2:1882480/4‑150 (MQ=255)
                                                                                                                                        tggagcgccgggCTGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCa               >  1:608106/4‑150 (MQ=255)
                                                                                                                                         ccagcgccagcCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCaa              >  1:2235793/1‑150 (MQ=255)
                                                                                                                                               ccagcCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTgc        <  1:769127/150‑1 (MQ=255)
                                                                                                                                                    ccGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAACATTGCGCagt   <  2:979249/150‑3 (MQ=255)
                                                                                                                                                     cGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAACATTGCGCagtg  <  1:998325/150‑4 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCGGGCAGAGCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATAT  >  NC_000913/3111674‑3111972

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