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.0% 1278.4 / ‑0.8 403F427L (TTC→CTC) ptsAPEP‑protein phosphotransferase system enzyme I/IIA
Reads supporting (aligned to +/- strand):  ref base A (4/0);  new base G (231/168);  total (235/168)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 406 total aligned reads displayed.

TATCGCGATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCG  >  NC_000913/4138067‑4138246
                                                                                          |                                                                                         
tATCGCGATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCa                                                                                         >  2:4531448/1‑93 (MQ=255)
tATCGCGATCCACCGCCAGCAAATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCa                                                                                         <  2:7549980/93‑1 (MQ=255)
  tCGCGATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATc                                                                                       <  2:4612727/93‑1 (MQ=255)
  tCGCGATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATc                                                                                       >  1:4442487/1‑93 (MQ=255)
   cgcgATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCa                                                                                      >  2:7299859/1‑93 (MQ=255)
    gcgATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCAc                                                                                     >  1:3679856/1‑93 (MQ=255)
       aTCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGa                                                                                  >  2:5882845/1‑93 (MQ=255)
       aTCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGa                                                                                  >  2:1636986/1‑93 (MQ=255)
       aTCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGa                                                                                  <  1:7863022/93‑1 (MQ=255)
        tCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGAc                                                                                 >  1:766128/1‑93 (MQ=255)
        tCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGAc                                                                                 >  1:7477827/1‑93 (MQ=255)
           cccGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGc                                                                              <  1:2952544/92‑1 (MQ=255)
           aCCGCAAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGc                                                                              >  1:1478033/1‑93 (MQ=255)
              gCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACt                                                                           >  1:1058301/1‑93 (MQ=255)
              gCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACt                                                                           <  1:7787955/93‑1 (MQ=255)
                cagcagATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTc                                                                         <  2:2657048/93‑1 (MQ=255)
                  gcagATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCa                                                                       >  1:5607669/1‑93 (MQ=255)
                  gcagATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCa                                                                       >  2:4187533/1‑93 (MQ=255)
                  gcagATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCa                                                                       <  1:4952696/93‑1 (MQ=255)
                  gcagATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCa                                                                       <  1:3594639/93‑1 (MQ=255)
                   cagATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAg                                                                      <  1:4655703/93‑1 (MQ=255)
                    agATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGc                                                                     >  2:4668757/1‑93 (MQ=255)
                     gATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCa                                                                    >  2:2099680/1‑93 (MQ=255)
                     gATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCa                                                                    <  2:2983698/93‑1 (MQ=255)
                      aTACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCat                                                                   <  1:3120935/93‑1 (MQ=255)
                      aTACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCat                                                                   >  2:7400580/1‑93 (MQ=255)
                      aTACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCat                                                                   >  1:3643484/1‑93 (MQ=255)
                       tACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCatg                                                                  <  1:4821723/93‑1 (MQ=255)
                        aCTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCatga                                                                 <  1:5751176/93‑1 (MQ=255)
                          tGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCatgatg                                                               >  2:1888762/1‑93 (MQ=255)
                           gCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGc                                                              <  1:2630603/93‑1 (MQ=255)
                             gTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCg                                                            >  1:6393626/1‑93 (MQ=255)
                               cAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAg                                                          <  1:2476450/93‑1 (MQ=255)
                                aGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGc                                                         <  1:599880/93‑1 (MQ=255)
                                 ggTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCt                                                        >  1:6822146/1‑93 (MQ=255)
                                  gTCATTACTACCAATACTAATGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTg                                                       >  1:678006/1‑93 (MQ=255)
                                  gTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTg                                                       >  1:4173710/1‑93 (MQ=255)
                                  gTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTg                                                       >  2:5455668/1‑93 (MQ=255)
                                   tCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTgg                                                      <  2:3587122/93‑1 (MQ=255)
                                    cATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGa                                                     >  2:5634438/1‑93 (MQ=255)
                                    cATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGa                                                     >  2:31153/1‑93 (MQ=255)
                                     aTTACTACCAATACTAAAGAAATCAATCTCTTCGGAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGAt                                                    >  2:4958335/1‑93 (MQ=255)
                                     aTTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGAt                                                    >  2:658527/1‑93 (MQ=255)
                                      ttACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGAtt                                                   >  2:5338750/1‑93 (MQ=255)
                                      ttACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGAtt                                                   >  2:631745/1‑93 (MQ=255)
                                      ttACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGAtt                                                   >  2:1070933/1‑93 (MQ=255)
                                       tactacCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGAttt                                                  >  1:7358856/1‑93 (MQ=255)
                                         ctacCAATCCTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTtc                                                <  2:1927708/93‑1 (MQ=255)
                                         ctacCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTtc                                                <  1:3339538/93‑1 (MQ=255)
                                          tacCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTtct                                               >  1:2064738/1‑93 (MQ=255)
                                          tacCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTtct                                               >  1:1403276/1‑93 (MQ=255)
                                            ccAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCtca                                             <  2:7865712/93‑1 (MQ=255)
                                            ccAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCtca                                             >  1:845383/1‑93 (MQ=255)
                                               aTACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCtcatca                                          >  1:3463084/1‑93 (MQ=255)
                                               aTACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCtcatca                                          <  2:4393385/93‑1 (MQ=255)
                                                 aCTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCaaa                                        <  2:5162170/93‑1 (MQ=255)
                                                   tAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAGACg                                      >  2:23855/1‑93 (MQ=255)
                                                   tAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACg                                      >  1:803445/1‑93 (MQ=255)
                                                   tAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACg                                      >  2:2218209/1‑93 (MQ=255)
                                                   tAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACg                                      >  1:1231561/1‑93 (MQ=255)
                                                     aaGAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGa                                    >  1:1034859/1‑93 (MQ=255)
                                                       gAAATCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAAt                                  >  2:6027982/1‑93 (MQ=255)
                                                        aaaTCAATCTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAAtg                                 >  1:4274020/1‑93 (MQ=255)
                                                        aaaTCAATCTATTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAAtg                                 <  1:4529961/93‑1 (MQ=255)
                                                         aatcaatcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAAtgt                                >  1:6051661/1‑93 (MQ=255)
                                                          atcaatcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAAtgtg                               >  1:2309694/1‑93 (MQ=255)
                                                          atcaatcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAAtgtg                               >  2:2596523/1‑93 (MQ=255)
                                                           tcaatcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATgtgt                              <  2:3688738/93‑1 (MQ=255)
                                                             aatcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTc                            <  1:4006681/93‑1 (MQ=255)
                                                              atcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCg                           >  1:1621973/1‑93 (MQ=255)
                                                              atcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCg                           >  2:4825976/1‑93 (MQ=255)
                                                              atcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCg                           <  2:1402626/93‑1 (MQ=255)
                                                              atcTCTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCg                           >  1:3013184/1‑93 (MQ=255)
                                                               tctcTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGt                          <  2:7054608/93‑1 (MQ=255)
                                                               tctcTTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGt                          >  1:5119921/1‑93 (MQ=255)
                                                                  ctTCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTAc                       <  1:7123054/93‑1 (MQ=255)
                                                                   ttCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACg                      >  2:331965/1‑93 (MQ=255)
                                                                    tCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGt                     <  1:3246112/93‑1 (MQ=255)
                                                                    tCGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGt                     <  2:102846/93‑1 (MQ=255)
                                                                     cGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTa                    <  1:7218748/93‑1 (MQ=255)
                                                                     cGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTa                    <  2:3980213/93‑1 (MQ=255)
                                                                     cGCAGCATTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTa                    <  1:615170/93‑1 (MQ=255)
                                                                      gcagcaTTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAg                   >  2:126513/1‑93 (MQ=255)
                                                                      gcagcaTTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAg                   >  2:2740/1‑93 (MQ=255)
                                                                       cagcaTTGATCGATGATGAGCATCACCGACGTCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGt                  <  1:3569428/93‑1 (MQ=255)
                                                                       cagcaTTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGt                  <  2:6199829/93‑1 (MQ=255)
                                                                       cagcaTTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGt                  <  2:6856419/93‑1 (MQ=255)
                                                                        agcaTTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGtt                 >  2:7344312/1‑93 (MQ=255)
                                                                         gcaTTGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTg                >  1:2158440/1‑93 (MQ=255)
                                                                         gcaTTGATCGATGATGAGCATCACCGACAGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAAAGGAATGTGTTCGGTACGTattcg                >  2:5718161/1‑89 (MQ=255)
                                                                            ttGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTg             <  2:526895/93‑1 (MQ=255)
                                                                            ttGATCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTg             <  1:6676638/93‑1 (MQ=255)
                                                                               aTCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGttt          <  1:379374/93‑1 (MQ=255)
                                                                               aTCGATAATGAGCATAACCGACGGCACATCCAGCATGATGACGAGCTGGAATTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGttt          <  2:7003176/93‑1 (MQ=255)
                                                                                tCGATGATGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTg         >  1:6557394/1‑93 (MQ=255)
                                                                                  gatgatGAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGc       <  1:5947803/93‑1 (MQ=255)
                                                                                   atgatgAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCt      <  1:2706108/93‑1 (MQ=255)
                                                                                   atgatgAGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCt      >  1:2858281/1‑93 (MQ=255)
                                                                                     gatgaGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTc    >  1:3893478/1‑93 (MQ=255)
                                                                                      atgaGCATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTcc   <  2:6502704/93‑1 (MQ=255)
                                                                                          |                                                                                         
TATCGCGATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCG  >  NC_000913/4138067‑4138246

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