Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,882,550 G→A L394F (CTC→TTC)  dnaA ← chromosomal replication initiator protein DnaA, DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,882,5500GA98.9% 608.9 / NA 175L394F (CTC→TTC) dnaAchromosomal replication initiator protein DnaA, DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (70/103);  minor base C (1/1);  total (71/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 176 total aligned reads displayed.

CAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGC  >  NC_000913/3882454‑3882646
                                                                                                |                                                                                                
    ctcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCcgcg                                                                                           <  2:2299129/100‑1 (MQ=255)
    ctcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCcgcg                                                                                           <  2:1918971/100‑1 (MQ=255)
     tcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGa                                                                                          >  2:2167489/1‑100 (MQ=255)
     tcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGa                                                                                          >  2:3869789/1‑100 (MQ=255)
      ccGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGCGCGGGATCGACGCTTGGAAAGGAAATCCGCGAc                                                                                         <  2:2318943/100‑1 (MQ=255)
      ccGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGAc                                                                                         <  2:2690012/100‑1 (MQ=255)
      ccGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGAc                                                                                         <  1:798750/100‑1 (MQ=255)
      ccGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGAc                                                                                         <  2:3814948/100‑1 (MQ=255)
       cGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACt                                                                                        <  1:3208994/100‑1 (MQ=255)
         gCAGACTGTGGTTAGTCAGCTCTTTCTCCAGCGCCAGCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACttt                                                                                      <  2:3452624/100‑1 (MQ=255)
          cAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTg                                                                                     <  2:3068592/100‑1 (MQ=255)
            gACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGAt                                                                                   >  1:1535867/1‑100 (MQ=255)
             aCTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATc                                                                                  <  1:1478520/100‑1 (MQ=255)
              cTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCt                                                                                 <  2:1435036/100‑1 (MQ=255)
              cTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCt                                                                                 <  1:935146/100‑1 (MQ=255)
                gtgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTg                                                                               <  2:598421/100‑1 (MQ=255)
                 tgattAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTgt                                                                              <  2:1881995/97‑1 (MQ=255)
                 tgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTgt                                                                              <  2:3862131/100‑1 (MQ=255)
                  ggTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTgta                                                                             >  1:3698728/1‑100 (MQ=255)
                   gTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTgtag                                                                            <  1:814400/100‑1 (MQ=255)
                   gTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTgtag                                                                            <  2:1751099/100‑1 (MQ=255)
                    ttAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTgtagt                                                                           <  1:1311568/100‑1 (MQ=255)
                    ttAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCCCCACCGAGTGGGATCGATGCTTGGAAAGGAAAGCCGCGACTTTGATCTTgtag                                                                            <  2:1604973/99‑1 (MQ=255)
                     tAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTgtagta                                                                          >  1:682446/1‑100 (MQ=255)
                      aGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTAc                                                                         >  2:890135/1‑100 (MQ=255)
                      aGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTAc                                                                         >  2:2019573/1‑100 (MQ=255)
                       gTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGACATCCGCGACTTTGATCTTGTAGTACt                                                                        >  2:25952/1‑100 (MQ=255)
                       gTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACt                                                                        >  2:908986/1‑100 (MQ=255)
                       gTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACt                                                                        >  1:191142/1‑100 (MQ=255)
                            ctctTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGcc                                                                   <  1:300924/100‑1 (MQ=255)
                             tctTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCa                                                                  <  2:1365754/100‑1 (MQ=255)
                              ctTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCAc                                                                 >  1:1740448/1‑100 (MQ=255)
                              ctTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCAc                                                                 >  2:582965/1‑100 (MQ=255)
                                 tCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGt                                                              >  1:3433911/1‑100 (MQ=255)
                                 tCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGt                                                              >  2:3572940/1‑100 (MQ=255)
                                   ggcAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtct                                                            >  1:2607709/3‑100 (MQ=255)
                                   gCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtct                                                            >  1:97042/1‑100 (MQ=255)
                                   gCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtct                                                            <  1:1013945/100‑1 (MQ=255)
                                    ccAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtctt                                                           >  1:2388510/1‑100 (MQ=255)
                                     cAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtcttc                                                          >  2:213829/1‑100 (MQ=255)
                                     cAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtcttc                                                          <  2:3444219/100‑1 (MQ=255)
                                       gcgcCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTg                                                        <  1:524186/100‑1 (MQ=255)
                                       gcgcCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTg                                                        <  2:3424213/100‑1 (MQ=255)
                                       gcgcCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTg                                                        <  2:2204830/100‑1 (MQ=255)
                                          ccatcgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAat                                                     >  2:2878313/1‑100 (MQ=255)
                                          ccatcgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAat                                                     <  1:742695/100‑1 (MQ=255)
                                          ccatcgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAat                                                     >  1:3861077/1‑100 (MQ=255)
                                          ccatcgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAat                                                     <  1:3383766/100‑1 (MQ=255)
                                           catcgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAata                                                    <  1:3007288/100‑1 (MQ=255)
                                            atcgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAatat                                                   >  2:3669425/1‑100 (MQ=255)
                                            atcgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAatat                                                   >  1:1633867/1‑100 (MQ=255)
                                              cgccatcTTGCCCGGACGCGCCACCGAGCGGGATCTACGCCTGTAAAGGAAATCCGCGACTTTGATCTTGTACTACTCCGCCACCGTCTTCTGAATATTg                                                 <  2:686368/100‑1 (MQ=255)
                                              cgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTg                                                 <  2:1535862/100‑1 (MQ=255)
                                               gccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGt                                                <  1:3300493/100‑1 (MQ=255)
                                                 catcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCg                                              <  2:889924/100‑1 (MQ=255)
                                                  atcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGa                                             <  1:3587653/100‑1 (MQ=255)
                                                   tcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCTGCCACCGTCTTCTGAATATTGTCGAt                                            <  1:2354791/100‑1 (MQ=255)
                                                   tcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGAt                                            <  1:3861148/100‑1 (MQ=255)
                                                    cTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATg                                           <  2:3751229/100‑1 (MQ=255)
                                                     tGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATgg                                          >  1:1224633/1‑100 (MQ=255)
                                                     tGGCGCGGACGCGCCACCGACCGGGATCGACTCCTGGAAAGGCAATCCGCGACTTTGACCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATgg                                          <  2:1682288/100‑1 (MQ=255)
                                                      ggCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGt                                         >  1:2686342/1‑100 (MQ=255)
                                                       gcgcGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTg                                        <  2:3097954/100‑1 (MQ=255)
                                                        cgcgGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGa                                       >  1:1269434/1‑100 (MQ=255)
                                                        cgcgGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGCCTTCTGAATATTGTCGATGGTg                                        >  2:735705/1‑99 (MQ=255)
                                                         gcgGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTACTACTCCGCCACCGTCTTCTGAATACTGTCGATGGTGAc                                      <  2:2388501/100‑1 (MQ=255)
                                                          cgGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGAcc                                     >  1:2385425/1‑100 (MQ=255)
                                                           ggACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCa                                    >  1:1164920/1‑100 (MQ=255)
                                                           ggACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCa                                    >  1:3232866/1‑100 (MQ=255)
                                                            gACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAg                                   <  1:582966/100‑1 (MQ=255)
                                                               gcgcCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGttt                                <  1:179266/100‑1 (MQ=255)
                                                                cgcCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGtttt                               <  2:97041/100‑1 (MQ=255)
                                                                cgcCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGtttt                               >  2:1194462/1‑100 (MQ=255)
                                                                cgcCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGtttt                               <  2:116755/100‑1 (MQ=255)
                                                                  ccACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGttttt                              <  1:1666691/99‑1 (MQ=255)
                                                                  ccACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTc                             <  1:1504424/100‑1 (MQ=255)
                                                                  ccACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTc                             >  1:3187688/1‑100 (MQ=255)
                                                                  ccACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTc                             >  1:3556098/1‑100 (MQ=255)
                                                                   cACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTcc                            <  2:3590808/100‑1 (MQ=255)
                                                                    aCCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCt                           <  1:2818737/100‑1 (MQ=255)
                                                                     ccGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTg                          <  1:2167498/100‑1 (MQ=255)
                                                                        aGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCaa                       >  1:999359/1‑100 (MQ=255)
                                                                          cGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATg                     <  1:114222/100‑1 (MQ=255)
                                                                          cGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATg                     <  1:3669442/100‑1 (MQ=255)
                                                                            ggATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGcc                   >  2:2454242/1‑100 (MQ=255)
                                                                              aTCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAg                 >  1:639942/1‑100 (MQ=255)
                                                                               tCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGc                <  1:2930003/100‑1 (MQ=255)
                                                                               tCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTCTAGTGCTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGc                <  1:2721327/100‑1 (MQ=255)
                                                                                cGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCa               <  2:3861060/100‑1 (MQ=255)
                                                                                cGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCa               <  2:3013842/100‑1 (MQ=255)
                                                                                 gACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCaa              >  2:423202/1‑100 (MQ=255)
                                                                                  aCGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAg             <  2:1258002/100‑1 (MQ=255)
                                                                                  aCGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAg             <  2:935914/100‑1 (MQ=255)
                                                                                   cGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGt            >  1:2457420/1‑100 (MQ=255)
                                                                                     cTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTcg          >  1:876889/1‑100 (MQ=255)
                                                                                     cTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTcg          <  1:908988/100‑1 (MQ=255)
                                                                                        ggaaaggaaaTCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCgcgc       <  2:191140/100‑1 (MQ=255)
                                                                                         gAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCgcgc       <  2:3338173/99‑1 (MQ=255)
                                                                                           aaGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAgc    <  1:407304/100‑1 (MQ=255)
                                                                                            aGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAgcg   <  1:2673340/100‑1 (MQ=255)
                                                                                                |                                                                                                
CAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGC  >  NC_000913/3882454‑3882646

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