Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,882,550 G→A 100% 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,5500GA72.1% 214.1 / 102.8 147L394F (CTC→TTC) dnaAchromosomal replication initiator protein DnaA, DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base G (22/19);  new base A (57/48);  total (80/67)
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 150 total aligned reads displayed.

CAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCC  >  NC_000913/3882454‑3882647
                                                                                                |                                                                                                 
cAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATc                                                                                                <  2:1113318/100‑1 (MQ=255)
   tctcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCcgc                                                                                             >  1:1720185/1‑100 (MQ=255)
   tctcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCcgc                                                                                             >  1:625560/1‑100 (MQ=255)
   tctcCGGCAGACTGTGATTAGTCAGCTCTTTCGCCAGCGCCATCGCGATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCcgc                                                                                             <  2:1823780/100‑1 (MQ=255)
     tcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGa                                                                                           >  1:1802461/1‑100 (MQ=255)
     tcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGa                                                                                           >  1:768283/1‑100 (MQ=255)
      ccGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGAc                                                                                          >  1:2152205/1‑100 (MQ=255)
       cGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACt                                                                                         >  1:2807295/1‑100 (MQ=255)
       cGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACt                                                                                         >  2:163794/1‑100 (MQ=255)
        ggCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACtt                                                                                        >  1:857343/1‑100 (MQ=255)
         gCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACttt                                                                                       >  1:72894/1‑100 (MQ=255)
          cAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTg                                                                                      >  1:2506950/1‑100 (MQ=255)
            gACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGAt                                                                                    >  2:2926840/1‑100 (MQ=255)
             aCTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATc                                                                                   <  1:1514271/100‑1 (MQ=255)
              cTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCt                                                                                  >  1:128342/1‑100 (MQ=255)
              cTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCt                                                                                  <  1:21739/100‑1 (MQ=255)
              cTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCt                                                                                  <  1:21779/100‑1 (MQ=255)
                gttgTTAGTCAGCTCTTGCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTg                                                                                <  1:2127974/97‑1 (MQ=255)
                 tgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTgt                                                                               <  2:2506067/100‑1 (MQ=255)
                 tgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTgt                                                                               <  2:1455176/100‑1 (MQ=255)
                  ggTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgta                                                                              >  1:1250901/1‑100 (MQ=255)
                  ggTTAGTCAGCCCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGCCGCTTGGAAACGAAATCCGCGACCTTGATccagta                                                                              >  2:1626401/1‑95 (MQ=255)
                   gTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgtag                                                                             <  2:2870144/100‑1 (MQ=255)
                   gTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTgtag                                                                             <  2:2518223/100‑1 (MQ=255)
                   gTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTgtag                                                                             <  2:2937543/100‑1 (MQ=255)
                       gTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACt                                                                         >  1:1959702/1‑100 (MQ=255)
                         cttccgatctgCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTcc                                                                       <  1:1964086/90‑1 (MQ=255)
                         caactctTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTcc                                                                       >  1:547524/4‑100 (MQ=255)
                         cAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGGTCGACGCTTGGAAGGGGGATCCGCGACTTTGATCTTgtagtagt                                                                         >  2:2334018/1‑96 (MQ=255)
                           gCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGc                                                                     >  2:1556207/1‑100 (MQ=255)
                           gCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGc                                                                     <  1:250111/100‑1 (MQ=255)
                           gCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGc                                                                     >  1:1825052/1‑100 (MQ=255)
                            ctctTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCTCCACCGAGCGGGATCGTCGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGcc                                                                    <  2:2923024/100‑1 (MQ=255)
                             tctTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCa                                                                   >  2:323861/1‑100 (MQ=255)
                                ttCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGTTCTTGTAGTACTCCGCCACCg                                                                <  2:669607/100‑1 (MQ=255)
                                  cGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtc                                                              >  1:718323/1‑100 (MQ=255)
                                   gCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtct                                                             >  2:2600777/1‑100 (MQ=255)
                                   gCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtct                                                             >  2:796479/1‑100 (MQ=255)
                                   gCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTccagatcggaag                                                             >  2:1964069/1‑90 (MQ=255)
                                   gCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtct                                                             >  1:241911/1‑100 (MQ=255)
                                   gCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtct                                                             >  2:683234/1‑100 (MQ=255)
                                   gCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtct                                                             <  2:241909/100‑1 (MQ=255)
                                    ccAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGAAGATCCGCGACTTTGATCTTGTAGTAGTCAGcaccagcctt                                                            >  2:2419283/1‑91 (MQ=255)
                                     cAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtcttc                                                           >  2:288586/1‑100 (MQ=255)
                                      aGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtcttct                                                          <  2:3022969/100‑1 (MQ=255)
                                      aGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtcttct                                                          <  1:401188/100‑1 (MQ=255)
                                        cgccatcgccatCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGa                                                        <  1:307530/100‑1 (MQ=255)
                                          ccatcgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAat                                                      <  2:180459/100‑1 (MQ=255)
                                          ccatcgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAat                                                      <  1:491711/100‑1 (MQ=255)
                                           tcccgctaggtGGCGCGTCCGCGCCACCGAGCGGGATCGAAGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAata                                                     >  2:2128788/11‑100 (MQ=255)
                                               gccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGt                                                 >  1:1617000/1‑100 (MQ=255)
                                               gccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACGTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGt                                                 <  2:1575145/100‑1 (MQ=255)
                                                   ggtGGCGCGTCCGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGCCTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGAt                                             >  1:1273119/3‑100 (MQ=255)
                                                   ggtGGCGCGTCCGCGCCAACGAGCGGGATCGACGCTTGGAAAGGGAATCCGCGACTTTGGTGTTGTAGTACTCCGACACCGTCTTCTGAATATTGTCGAt                                             >  2:367170/3‑100 (MQ=255)
                                                    cTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATg                                            <  2:1276452/100‑1 (MQ=255)
                                                     tGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATgg                                           <  2:806929/100‑1 (MQ=255)
                                                      ggCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGt                                          >  2:2328867/1‑100 (MQ=255)
                                                      ggCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGt                                          >  2:2724144/1‑100 (MQ=255)
                                                        cgcgGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGa                                        >  1:2560661/1‑100 (MQ=255)
                                                        cgcgGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTAGATCTTGTCGTATGCGGCGACCGTCATCTGAATAGTGTCgcaggt                                          >  2:1848970/1‑93 (MQ=255)
                                                         gcgGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGAc                                       <  1:42431/100‑1 (MQ=255)
                                                           ggACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTg                                         <  2:3068350/96‑1 (MQ=255)
                                                                cgcCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGtttt                                <  1:612036/100‑1 (MQ=255)
                                                                   cACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTcc                             <  1:2671433/100‑1 (MQ=255)
                                                                   cACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTcc                             >  1:76963/1‑100 (MQ=255)
                                                                   cACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTcc                             <  2:1608513/100‑1 (MQ=255)
                                                                   cACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTcc                             >  2:1161134/1‑100 (MQ=255)
                                                                    aCCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCt                            <  1:1036359/100‑1 (MQ=255)
                                                                     ccGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCt                            <  1:829325/99‑1 (MQ=255)
                                                                      cGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGc                          >  1:1147928/1‑100 (MQ=255)
                                                                       gAGCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCa                         <  1:1412909/100‑1 (MQ=255)
                                                                         gCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAAt                       >  2:1303852/1‑100 (MQ=255)
                                                                         gCGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAAt                       <  1:2722424/100‑1 (MQ=255)
                                                                          cGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATg                      <  2:2193821/100‑1 (MQ=255)
                                                                          cGGGATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATg                      <  2:2627961/100‑1 (MQ=255)
                                                                           gggATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGc                     >  1:3070064/1‑100 (MQ=255)
                                                                            ggATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGcc                    <  1:1765174/100‑1 (MQ=255)
                                                                            ggATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGcc                    >  1:166240/1‑100 (MQ=255)
                                                                            ggATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGcc                    <  1:1290583/100‑1 (MQ=255)
                                                                             gATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCa                   <  2:180126/100‑1 (MQ=255)
                                                                             gATCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCa                   >  2:1266803/1‑100 (MQ=255)
                                                                              aTCGACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAg                  <  2:1340936/100‑1 (MQ=255)
                                                                               tCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGc                 >  2:1693573/1‑100 (MQ=255)
                                                                                 gACGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCaa               >  1:921027/1‑100 (MQ=255)
                                                                                   cGCTTGGAAAGGCAATCCGCGACTTTGATCTTGTCGTACTCCGCCACCGTCGTCTGAATATTGTCGATGGTACACGGTTTTTCCTGCCATGCCGGCaa               >  2:2009588/1‑98 (MQ=255)
                                                                                   cGCTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGt             >  1:1426208/1‑100 (MQ=255)
                                                                                     cTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTcg           >  1:1586524/1‑100 (MQ=255)
                                                                                     cTTGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTcg           <  2:2829615/100‑1 (MQ=255)
                                                                                      ttGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTcgc          <  1:2328889/100‑1 (MQ=255)
                                                                                      ttGGAAAGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATACCAGCAAGTcgc          <  2:3029032/100‑1 (MQ=255)
                                                                                        ggaaaggaaaTCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCgcgc        <  2:2476532/100‑1 (MQ=255)
                                                                                        ggaaaggaaaTCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCgcgc        >  1:1070003/1‑100 (MQ=255)
                                                                                        ggAAAGGAGATCCGCGCCTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCgcgc        <  2:148339/100‑1 (MQ=255)
                                                                                          aaaGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAg      >  1:1721776/1‑100 (MQ=255)
                                                                                          aaaGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAg      >  1:117804/1‑100 (MQ=255)
                                                                                           aaGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAgc     >  2:392362/1‑100 (MQ=255)
                                                                                            aGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGAGCAGTTTTTCCTGCAGTGCGAGCAAGTCGCGCAg      >  2:1024310/1‑98 (MQ=255)
                                                                                            aGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAgcg    <  1:2725155/100‑1 (MQ=255)
                                                                                            aGGAAATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAgcg    >  1:2495310/1‑100 (MQ=255)
                                                                                              gagaTCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGcc  <  1:1879847/100‑1 (MQ=255)
                                                                                                |                                                                                                 
CAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCC  >  NC_000913/3882454‑3882647

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