Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 216,244 A→G Y254H (TAT→CAT)  yaeF ← peptidase C92 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913216,2440AG87.7% 1213.4 / 134.3 383Y254H (TAT→CAT) yaeFpeptidase C92 family protein
Reads supporting (aligned to +/- strand):  ref base A (23/24);  new base G (186/150);  total (209/174)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 383 total aligned reads displayed.

CTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATACTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCAGCCCGATTGCGCCAGTGTAAGTGG  >  NC_000913/216145‑216343
                                                                                                   |                                                                                                   
 tCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCt                                                                                                   <  1:404966/101‑1 (MQ=255)
 tCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCt                                                                                                   >  2:405528/1‑101 (MQ=255)
  cATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTg                                                                                                  >  2:405533/1‑101 (MQ=255)
   aTACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGt                                                                                                 <  1:404882/101‑1 (MQ=255)
   aTACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGt                                                                                                 <  2:404912/101‑1 (MQ=255)
   aTACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATACTGt                                                                                                 >  1:405537/1‑101 (MQ=255)
    tACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATACTGTa                                                                                                >  1:405540/1‑101 (MQ=255)
     aCCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTaa                                                                                               >  2:405545/1‑101 (MQ=255)
      ccTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAAc                                                                                              <  1:404785/101‑1 (MQ=255)
      ccTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAAc                                                                                              >  1:405547/1‑101 (MQ=255)
       cTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACt                                                                                             <  2:405020/101‑1 (MQ=255)
        ttAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTg                                                                                            <  1:404754/101‑1 (MQ=255)
        ttAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTg                                                                                            >  1:405552/1‑101 (MQ=255)
         tAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGc                                                                                           >  2:405555/1‑101 (MQ=255)
         tAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGc                                                                                           >  1:405558/1‑101 (MQ=255)
          aaGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCg                                                                                          >  2:405563/1‑101 (MQ=255)
           aGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGt                                                                                         >  1:405567/1‑101 (MQ=255)
               ggATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATACTGTAACTGCGTTTCt                                                                                     >  2:405573/1‑101 (MQ=255)
                gATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTg                                                                                    <  1:404925/101‑1 (MQ=255)
                gATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTg                                                                                    <  1:405246/101‑1 (MQ=255)
                 aTAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTgg                                                                                   >  1:405577/1‑101 (MQ=255)
                  tAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGc                                                                                  >  2:405579/1‑101 (MQ=255)
                   aaGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCt                                                                                 <  1:404851/101‑1 (MQ=255)
                   aaGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATACTGTAACTGCGTTTCTGGCt                                                                                 <  2:404934/101‑1 (MQ=255)
                    aGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCtt                                                                                >  1:405584/1‑101 (MQ=255)
                      gCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTaa                                                                              >  2:405586/1‑101 (MQ=255)
                       cATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTaaa                                                                             >  1:405589/1‑101 (MQ=255)
                          ttACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGc                                                                          <  1:404944/101‑1 (MQ=255)
                          ttACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGc                                                                          >  2:405595/1‑101 (MQ=255)
                            aCTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCag                                                                        >  2:405601/1‑101 (MQ=255)
                            aCTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATACTGTAACTGCGTTTCTGGCTTAAACGCag                                                                        >  1:405602/1‑101 (MQ=255)
                               ttATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAgaga                                                                     >  2:405607/1‑101 (MQ=255)
                                tATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGAc                                                                    >  1:405609/1‑101 (MQ=255)
                                  tCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACAt                                                                  >  2:405616/1‑101 (MQ=255)
                                   cACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATc                                                                 >  1:405617/1‑101 (MQ=255)
                                    aacGGGTCAAACCAACAAAGCGACCCACCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCa                                                                >  1:405629/3‑101 (MQ=255)
                                    aCCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCa                                                                >  2:405625/1‑101 (MQ=255)
                                    aCCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATACTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCa                                                                >  1:405624/1‑101 (MQ=255)
                                     ccGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCa                                                                <  2:405121/100‑1 (MQ=255)
                                     ccGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATACTGTAACTGCGTTTCTgg                                                                                   <  2:405007/81‑1 (MQ=255)
                                      cGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCAcc                                                              <  1:405088/101‑1 (MQ=255)
                                        ggTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCaa                                                            <  2:405148/101‑1 (MQ=255)
                                         gTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTaa                                                                                               <  2:404759/65‑1 (MQ=255)
                                         gTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAAt                                                           >  2:405635/1‑101 (MQ=255)
                                          tCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATa                                                          >  1:405638/1‑101 (MQ=255)
                                           cAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATACTGTAACTGCGtt                                                                                        <  2:404831/70‑1 (MQ=255)
                                            aaaCCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACg                                                        >  2:405641/1‑101 (MQ=255)
                                            aaaCCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACg                                                        <  1:405386/101‑1 (MQ=255)
                                               ccAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCat                                                     >  1:405646/1‑101 (MQ=255)
                                               ccAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCat                                                     <  1:405036/101‑1 (MQ=255)
                                                caacaaAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCata                                                    <  2:404717/101‑1 (MQ=255)
                                                 aacaaAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCatat                                                   >  1:405651/1‑101 (MQ=255)
                                                 aacaaAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCatat                                                   <  1:404536/101‑1 (MQ=255)
                                                  acaaAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATg                                                  <  1:405098/101‑1 (MQ=255)
                                                   caaAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGc                                                 >  2:405655/1‑101 (MQ=255)
                                                   caaAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGc                                                 <  1:404996/101‑1 (MQ=255)
                                                    aaaGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCa                                                >  1:405659/1‑101 (MQ=255)
                                                     aaGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCAt                                               <  2:405221/101‑1 (MQ=255)
                                                     aaGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATACTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCAt                                               >  2:405660/1‑101 (MQ=255)
                                                      aGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATc                                              <  1:405147/101‑1 (MQ=255)
                                                      aGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATc                                              <  2:404883/101‑1 (MQ=255)
                                                         gACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGa                                           >  1:405667/1‑101 (MQ=255)
                                                          aCCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGAt                                          >  2:405668/1‑101 (MQ=255)
                                                           cccGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATc                                         <  1:405006/101‑1 (MQ=255)
                                                           cccGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATc                                         <  1:405169/101‑1 (MQ=255)
                                                            ccGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCg                                        <  2:405099/101‑1 (MQ=255)
                                                             cGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCgg                                       <  2:405354/101‑1 (MQ=255)
                                                             cGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCgg                                       <  2:404842/101‑1 (MQ=255)
                                                             cGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCgg                                       <  1:405199/101‑1 (MQ=255)
                                                              gCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGc                                      <  2:405044/101‑1 (MQ=255)
                                                              gCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGc                                      >  2:405673/1‑101 (MQ=255)
                                                              gCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGc                                      <  1:404982/101‑1 (MQ=255)
                                                              gCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGc                                      <  2:404970/101‑1 (MQ=255)
                                                               ccTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCa                                     >  1:405675/1‑101 (MQ=255)
                                                               ccTTGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCa                                     <  1:405109/101‑1 (MQ=255)
                                                                 ttGATATAAATCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAgg                                   >  1:405682/1‑101 (MQ=255)
                                                                 ttGATATAAATCCCAGGCTTCAGGTGCCCAACATACTGGAACTGCGtt                                                                                        <  2:404814/48‑1 (MQ=255)
                                                                     tataAATCCCAGGCTTCAGGTGCCCAACATACTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTa                               >  2:405686/1‑101 (MQ=255)
                                                                      ataAATCCCAGGCTTCAGGTGCCCAACATACTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTaa                              >  1:405689/1‑101 (MQ=255)
                                                                        aaaTCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATc                            >  2:405692/1‑101 (MQ=255)
                                                                        aaaTCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATc                            <  1:405269/101‑1 (MQ=255)
                                                                         aaTCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATcc                           <  2:404981/101‑1 (MQ=255)
                                                                          aTCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCa                          <  2:404899/101‑1 (MQ=255)
                                                                          aTCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCa                          >  2:405698/1‑101 (MQ=255)
                                                                          aTCCCAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCa                          <  1:405026/101‑1 (MQ=255)
                                                                            cccAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCAGc                        <  2:405033/101‑1 (MQ=255)
                                                                              cAGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCAGccc                      >  2:405701/1‑101 (MQ=255)
                                                                               aGGCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCAGCCCg                     <  2:405065/101‑1 (MQ=255)
                                                                                ggCTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCAGCCCGa                    >  1:405703/1‑101 (MQ=255)
                                                                                  cTTCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCAGCCCGAtt                  >  2:405706/1‑101 (MQ=255)
                                                                                   ttCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCAGCCCGATTg                 >  1:405714/1‑101 (MQ=255)
                                                                                   ttCAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCAGCCCGATTg                 >  2:405712/1‑101 (MQ=255)
                                                                                     cAGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCAGCCCGATTgcg               <  2:404368/101‑1 (MQ=255)
                                                                                      aGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCAGcc                       >  2:405722/1‑92 (MQ=255)
                                                                                      aGGTGCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCAGCCCGATTgcgc              <  1:405061/101‑1 (MQ=255)
                                                                                      aGGTGCCCAACATACTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCAGCCCGATTgcgc              >  1:405721/1‑101 (MQ=255)
                                                                                          gCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCAGCCCGATTGCGCCAgt          <  2:405414/101‑1 (MQ=255)
                                                                                          gCCCAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCAGCCCGATTGCGCCAgt          >  1:405727/1‑101 (MQ=255)
                                                                                            ccAACATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCAGccc                      >  2:405731/1‑87 (MQ=255)
                                                                                              aaCATGCTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCAGCCCGATTGCGCCAGTGTaa      >  1:405734/1‑101 (MQ=255)
                                                                                                   |                                                                                                   
CTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGCTTCAGGTGCCCAACATACTGTAACTGCGTTTCTGGCTTAAACGCAGAGACATCACCAATACGCATATGCATCAGATCGGCAGGGCTAATCCAGCCCGATTGCGCCAGTGTAAGTGG  >  NC_000913/216145‑216343

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: