Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,138,157 A→G F427L (TTC→CTC)  ptsA ← PEP‑protein phosphotransferase system enzyme I/IIA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,138,1570AG99.4% 988.1 / ‑4.0 308F427L (TTC→CTC) ptsAPEP‑protein phosphotransferase system enzyme I/IIA
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (176/130);  total (176/132)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 313 total aligned reads displayed.

TATCGCGATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCG  >  NC_000913/4138067‑4138246
                                                                                          |                                                                                         
 aTCGCGATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCAt                                                                                        <  2:4028548/93‑1 (MQ=255)
   cgcgATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCa                                                                                      <  1:533808/93‑1 (MQ=255)
    gcgATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCAc                                                                                     <  2:1097710/93‑1 (MQ=255)
    gcgATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCAc                                                                                     >  2:3082541/1‑93 (MQ=255)
     cgATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCAcc                                                                                    <  2:4651097/93‑1 (MQ=255)
      gATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCg                                                                                   >  1:1036694/1‑93 (MQ=255)
       aTCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGa                                                                                  <  1:891071/93‑1 (MQ=255)
        tCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGAc                                                                                 <  1:1558479/93‑1 (MQ=255)
        tCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGAc                                                                                 >  2:1989335/1‑93 (MQ=255)
           aCCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGc                                                                              <  1:3307435/93‑1 (MQ=255)
           aCCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGc                                                                              <  1:622282/93‑1 (MQ=255)
           aCCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGc                                                                              >  2:2973047/1‑93 (MQ=255)
             cGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCAc                                                                            <  2:3210448/93‑1 (MQ=255)
               ccAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACtt                                                                          <  2:4296226/93‑1 (MQ=255)
                cagcagATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTc                                                                         >  2:3538377/1‑93 (MQ=255)
                cagcagATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTc                                                                         >  2:666356/1‑93 (MQ=255)
                 agcagATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTcc                                                                        >  1:1427066/1‑93 (MQ=255)
                     gATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCa                                                                    >  1:449420/1‑93 (MQ=255)
                      aTACTGCGTCAGGTCATTACTACCAATTCTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCat                                                                   <  1:1933318/93‑1 (MQ=255)
                      aTACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCat                                                                   >  1:3589101/1‑93 (MQ=255)
                       tACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCatg                                                                  >  1:1001875/1‑93 (MQ=255)
                       tACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCatg                                                                  >  1:2357887/1‑93 (MQ=255)
                        aCTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATATCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCatga                                                                 >  2:911986/1‑93 (MQ=255)
                         cTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCatgat                                                                <  1:269598/93‑1 (MQ=255)
                         cTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCatgat                                                                >  2:4227202/1‑93 (MQ=255)
                          tGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCatgatg                                                               >  2:1770874/1‑93 (MQ=255)
                            cGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGcc                                                             <  1:3870745/93‑1 (MQ=255)
                              tCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGa                                                           >  1:2056998/1‑93 (MQ=255)
                              tCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGa                                                           >  2:3817199/1‑93 (MQ=255)
                                aGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGc                                                         <  2:1889286/93‑1 (MQ=255)
                                aGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGc                                                         >  1:3465632/1‑93 (MQ=255)
                                 ggTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCt                                                        <  1:4665399/93‑1 (MQ=255)
                                   tCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTgg                                                      >  2:278465/1‑93 (MQ=255)
                                     aTTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGAt                                                    >  1:2776054/1‑93 (MQ=255)
                                     aTTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGAt                                                    >  1:2556452/1‑93 (MQ=255)
                                      ttACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGAtt                                                   <  1:1294925/93‑1 (MQ=255)
                                      ttACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGAtt                                                   >  2:2901760/1‑93 (MQ=255)
                                       tactacCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGTCACTTCCAGCATGATGCCGAGCTGGAttt                                                  >  1:2119365/1‑93 (MQ=255)
                                       tactacCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGAttt                                                  <  1:4528541/93‑1 (MQ=255)
                                       tactacCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGAttt                                                  >  1:1457863/1‑93 (MQ=255)
                                       tactacCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGAgtt                                                  >  2:2044504/1‑93 (MQ=255)
                                        actacCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGAtttt                                                 >  1:747184/1‑93 (MQ=255)
                                          tacCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTtct                                               >  1:3164396/1‑93 (MQ=255)
                                            ccAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCtca                                             <  2:220387/93‑1 (MQ=255)
                                            caaaTACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCtca                                             <  2:1229165/91‑1 (MQ=255)
                                             cAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCtcat                                            <  1:4400810/93‑1 (MQ=255)
                                             cAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCtcat                                            <  1:942974/93‑1 (MQ=255)
                                              aaTACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCtcatc                                           <  1:4218634/93‑1 (MQ=255)
                                              aaTACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCtcatc                                           >  1:3359995/1‑93 (MQ=255)
                                               aTACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCtcatca                                          <  1:1120897/93‑1 (MQ=255)
                                                 aCTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCaaa                                        >  1:68809/1‑93 (MQ=255)
                                                  cTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAAc                                       >  1:4650734/1‑93 (MQ=255)
                                                  cTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAAc                                       >  1:2597562/1‑93 (MQ=255)
                                                   tAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACg                                      <  2:335005/93‑1 (MQ=255)
                                                     aaGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGa                                    >  2:3674844/1‑93 (MQ=255)
                                                      aGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGaa                                   <  1:1897664/93‑1 (MQ=255)
                                                       gAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAAt                                  <  1:865696/93‑1 (MQ=255)
                                                        aaaTCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAAtg                                 <  2:4164174/93‑1 (MQ=255)
                                                         aatcaatcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAAtgt                                <  2:2195037/93‑1 (MQ=255)
                                                          atcaatcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAAtgtg                               <  2:3160513/93‑1 (MQ=255)
                                                          atcaatcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAAtgtg                               >  1:1890914/1‑93 (MQ=255)
                                                           tcaatcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATgtgt                              <  1:2966826/93‑1 (MQ=255)
                                                           tcaatcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATgtgt                              >  1:1702918/1‑93 (MQ=255)
                                                           tcaatcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATgtgt                              >  2:770094/1‑93 (MQ=255)
                                                            caatcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGtt                             <  2:2441163/93‑1 (MQ=255)
                                                             aatcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTc                            >  1:460988/1‑93 (MQ=255)
                                                             aatcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTc                            <  1:394048/93‑1 (MQ=255)
                                                              atcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCg                           <  1:372346/93‑1 (MQ=255)
                                                               tctcTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGt                          >  1:409673/1‑93 (MQ=255)
                                                                 tctTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTa                        <  2:1838771/93‑1 (MQ=255)
                                                                  ctTCGCAGCATTGATCGATGATTAGCACCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTAc                       <  1:4341391/93‑1 (MQ=255)
                                                                   ttCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACg                      >  2:4453796/1‑93 (MQ=255)
                                                                    tCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGt                     <  2:797510/93‑1 (MQ=255)
                                                                      gcagcaTTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAg                   >  1:996064/1‑93 (MQ=255)
                                                                       cagcaTTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGt                  >  2:1525121/1‑93 (MQ=255)
                                                                       cagcaTTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGt                  <  1:4054411/93‑1 (MQ=255)
                                                                        agcaTTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGtt                 <  2:3478464/93‑1 (MQ=255)
                                                                        agcaTTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGtt                 >  2:2270482/1‑93 (MQ=255)
                                                                         gcaTTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTg                >  2:517658/1‑93 (MQ=255)
                                                                         gcaTTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCAGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTg                >  2:3738380/1‑93 (MQ=255)
                                                                          caTTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGc               <  1:1075377/93‑1 (MQ=255)
                                                                          caTTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGc               >  2:2713842/1‑93 (MQ=255)
                                                                          caTTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGc               <  2:2518138/93‑1 (MQ=255)
                                                                            ttGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTg             >  1:1792515/1‑93 (MQ=255)
                                                                            ttGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTg             >  2:873473/1‑93 (MQ=255)
                                                                             tGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGt            >  1:3531327/1‑93 (MQ=255)
                                                                               aTCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGttt          >  1:1324574/1‑93 (MQ=255)
                                                                                tCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATgtgt                              <  1:1176968/72‑1 (MQ=255)
                                                                                tCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATgtgt                              >  2:1176968/1‑72 (MQ=255)
                                                                                tCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTg         <  1:1222084/93‑1 (MQ=255)
                                                                                tCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTg         <  2:413856/93‑1 (MQ=255)
                                                                                tCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTACGTTACGTAGTTGCTGTTTg         >  1:3125824/1‑93 (MQ=255)
                                                                                  gatgatGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGc       >  2:628879/1‑93 (MQ=255)
                                                                                  gatgatGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGc       <  1:3072035/93‑1 (MQ=255)
                                                                                   atgatgAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCt      >  2:1277415/1‑93 (MQ=255)
                                                                                   atgatgAGCATCACCGACGGCACTTCCAGCATGAGGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCt      <  2:1593457/93‑1 (MQ=255)
                                                                                    tgatgaGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCtt     >  1:583535/1‑93 (MQ=255)
                                                                                      atgaGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTcc   <  1:2169084/93‑1 (MQ=255)
                                                                                      atgaGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTcc   >  1:3556851/1‑93 (MQ=255)
                                                                                       tgaGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCg  >  2:3956764/1‑93 (MQ=255)
                                                                                          |                                                                                         
TATCGCGATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCG  >  NC_000913/4138067‑4138246

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: