Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,327,379 A→G K82R (AAA→AGA)  crfC → clamp‑binding sister replication fork colocalization protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,327,3790AG99.7% 1486.7 / NA 358K82R (AAA→AGA) crfCclamp‑binding sister replication fork colocalization protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (186/171);  total (186/172)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 358 total aligned reads displayed.

TAATATCGCCGAACGTCACGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGG  >  NC_000913/4327281‑4327476
                                                                                                  |                                                                                                 
tAATATCGCCGAACGTCACGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGAt                                                                                                 <  1:7380066/101‑1 (MQ=255)
  atatCGCCGAACGTCACGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCa                                                                                               >  2:7380697/1‑101 (MQ=255)
   tatCGCCGAACGTCACGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCa                                                                                               >  1:7380701/1‑100 (MQ=255)
     tCGCCGAACGTCACGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAacc                                                                                            <  2:7380253/101‑1 (MQ=255)
     tCGCCGAACGTCACGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAacc                                                                                            <  1:7380411/101‑1 (MQ=255)
     tCGCCGAACGTCACGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAacc                                                                                            >  2:7380706/1‑101 (MQ=255)
        ccGAACGTCACGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAaccacc                                                                                         >  1:7380713/1‑101 (MQ=255)
         cGAACGTCACGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAAccacca                                                                                        >  1:7380717/1‑101 (MQ=255)
         cGAACGTCACGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAAccacca                                                                                        <  2:7380166/101‑1 (MQ=255)
         cGAACGTCACGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAAccacca                                                                                        >  2:7380714/1‑101 (MQ=255)
            aCGTCACGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTa                                                                                     >  1:7380722/1‑101 (MQ=255)
             cGTCACGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTaa                                                                                    >  2:7380724/1‑101 (MQ=255)
              gTCACGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAAt                                                                                   <  1:7380129/101‑1 (MQ=255)
               tCACGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATg                                                                                  <  1:7379906/101‑1 (MQ=255)
                cACGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGc                                                                                 >  2:7380729/1‑101 (MQ=255)
                 aCGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGcc                                                                                >  1:7380733/1‑101 (MQ=255)
                 aCGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGcc                                                                                <  2:7380278/101‑1 (MQ=255)
                  cgcgATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATAAACCACCATTa                                                                                     >  2:7380734/1‑95 (MQ=255)
                   gcgATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCAt                                                                              <  1:7380303/101‑1 (MQ=255)
                    cgATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCAtt                                                                             <  2:7380326/101‑1 (MQ=255)
                      aTGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCAttgt                                                                           >  2:7380742/1‑101 (MQ=255)
                      aTGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCAttgt                                                                           >  1:7380743/1‑101 (MQ=255)
                       tGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCAttgtt                                                                          >  2:7380747/1‑101 (MQ=255)
                        gCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCAttgttg                                                                         <  2:7380240/101‑1 (MQ=255)
                        gCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCAttgttg                                                                         <  1:7380212/101‑1 (MQ=255)
                        gCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCAttgttg                                                                         <  1:7380120/101‑1 (MQ=255)
                         cTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTgg                                                                        >  2:7380750/1‑101 (MQ=255)
                          tGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGt                                                                       <  2:7380530/101‑1 (MQ=255)
                            aataatGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTAc                                                                     <  2:7380054/101‑1 (MQ=255)
                            aataatGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTAc                                                                     >  1:7380754/1‑101 (MQ=255)
                             ataatGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACg                                                                    <  2:7379935/101‑1 (MQ=255)
                             ataatGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACg                                                                    <  2:7380130/101‑1 (MQ=255)
                              taatGAGTTGCGTAAAATTTGCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACgg                                                                   >  2:7380764/1‑101 (MQ=255)
                              taatGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACgg                                                                   >  2:7380760/1‑101 (MQ=255)
                              taatGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACgg                                                                   >  1:7380761/1‑101 (MQ=255)
                               aatGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGa                                                                  <  2:7380092/101‑1 (MQ=255)
                                atGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAg                                                                 <  1:7380366/101‑1 (MQ=255)
                                 tGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAgg                                                                <  2:7380008/101‑1 (MQ=255)
                                  gAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGt                                                               >  2:7380771/1‑101 (MQ=255)
                                   aGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGtt                                                              >  1:7380772/1‑101 (MQ=255)
                                    gTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTc                                                             <  1:7380374/101‑1 (MQ=255)
                                     ttGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCt                                                            >  1:7380775/1‑101 (MQ=255)
                                      tGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTg                                                           <  1:7380053/101‑1 (MQ=255)
                                      tGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTg                                                           <  2:7380203/101‑1 (MQ=255)
                                       gCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGc                                                          >  1:7380780/1‑101 (MQ=255)
                                         gTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTa                                                                                     >  1:7380783/1‑72 (MQ=255)
                                          tAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCta                                                       <  1:7380034/101‑1 (MQ=255)
                                             aaTTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGtt                                                              <  1:7380195/91‑1 (MQ=255)
                                             aaTTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCtaatc                                                    >  2:7380787/1‑101 (MQ=255)
                                               tttCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCtaatcgt                                                  >  1:7380791/1‑101 (MQ=255)
                                               tttCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCtaatcgt                                                  >  1:7380796/1‑101 (MQ=255)
                                                ttCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCtaatcgta                                                 >  2:7380799/1‑101 (MQ=255)
                                                ttCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCtaatcgta                                                 <  2:7380227/101‑1 (MQ=255)
                                                 tCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCtaatcgtaa                                                >  1:7380801/1‑101 (MQ=255)
                                                   ccGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCtaatcgtaatc                                              <  1:7379705/101‑1 (MQ=255)
                                                   ccGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCtaatcgtaatc                                              >  1:7380806/1‑101 (MQ=255)
                                                    cGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCtaatcgtaatcg                                             >  2:7380809/1‑101 (MQ=255)
                                                     gACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGc                                            >  2:7380812/1‑101 (MQ=255)
                                                      aCTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGcc                                           <  2:7380257/101‑1 (MQ=255)
                                                       cTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGccc                                          >  2:7380817/1‑101 (MQ=255)
                                                       cTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGccc                                          <  1:7380150/101‑1 (MQ=255)
                                                        tGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCa                                         >  1:7380819/1‑101 (MQ=255)
                                                         ggAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCaa                                        <  1:7380319/101‑1 (MQ=255)
                                                         ggAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCaa                                        >  2:7380821/1‑101 (MQ=255)
                                                         gaaaaTGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCaa                                        <  1:7380282/99‑1 (MQ=255)
                                                          gAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAAt                                       <  2:7379834/101‑1 (MQ=255)
                                                            aaTGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGa                                     <  1:7380401/101‑1 (MQ=255)
                                                            aaTGGTGCTTGCAATTGTCGGTACCACGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGa                                     <  1:7380225/101‑1 (MQ=255)
                                                              tGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACt                                   <  2:7380139/101‑1 (MQ=255)
                                                               ggTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTg                                  >  2:7380828/1‑101 (MQ=255)
                                                                gTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTgc                                 <  2:7380305/101‑1 (MQ=255)
                                                                gTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTgc                                 >  1:7380829/1‑101 (MQ=255)
                                                                 tGCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTgcg                                <  1:7380459/101‑1 (MQ=255)
                                                                  gCTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTgcgc                               <  2:7380385/101‑1 (MQ=255)
                                                                   cTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCt                              >  2:7380835/1‑101 (MQ=255)
                                                                   cTTGCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCt                              >  1:7380836/1‑101 (MQ=255)
                                                                      gCAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGcc                           >  2:7380840/1‑101 (MQ=255)
                                                                       cAATTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCg                          <  1:7380601/101‑1 (MQ=255)
                                                                        aaTTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGa                         <  2:7380636/101‑1 (MQ=255)
                                                                         aTTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTTCCTAATCGTAATCGCCCAATGACTGCGCTGCCGAc                        <  1:7380328/101‑1 (MQ=255)
                                                                         aTTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGAc                        >  1:7380845/1‑101 (MQ=255)
                                                                         aTTGTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGAc                        <  2:7380333/101‑1 (MQ=255)
                                                                            gTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCt                     <  2:7380432/101‑1 (MQ=255)
                                                                            gTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCt                     >  1:7380850/1‑101 (MQ=255)
                                                                            gTCGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCt                     <  2:7380267/101‑1 (MQ=255)
                                                                              cGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTa                   >  2:7380852/1‑101 (MQ=255)
                                                                              cGGTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTa                   <  1:7380340/101‑1 (MQ=255)
                                                                               ggTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTAt                  >  2:7380855/1‑101 (MQ=255)
                                                                                gTACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTAtt                 <  1:7380181/101‑1 (MQ=255)
                                                                                 tACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTc                >  1:7380858/1‑101 (MQ=255)
                                                                                 tACCATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTc                <  1:7379987/101‑1 (MQ=255)
                                                                                   ccATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGc              >  1:7380863/1‑101 (MQ=255)
                                                                                    cATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGcc             <  2:7380478/101‑1 (MQ=255)
                                                                                    cATGAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGcc             >  1:7380866/1‑101 (MQ=255)
                                                                                       gAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATa          <  2:7380186/101‑1 (MQ=255)
                                                                                       gAAAGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATa          <  1:7379972/101‑1 (MQ=255)
                                                                                        aaaGCAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATAc         >  1:7380871/1‑101 (MQ=255)
                                                                                            cAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGccc     >  2:7380876/1‑101 (MQ=255)
                                                                                            cAGGGAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGccc     >  1:7380875/1‑101 (MQ=255)
                                                                                              gggAGATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCgg   >  1:7380895/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
TAATATCGCCGAACGTCACGCGATGCTGAATAATGAGTTGCGTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGG  >  NC_000913/4327281‑4327476

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: