Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,881,885 G→A intergenic (‑528/‑80) dnaA ← / → rpmH DnaA/50S ribosomal subunit protein L34

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,881,8850GA86.1% 863.5 / 85.4 312intergenic (‑528/‑80)dnaA/rpmHDnaA/50S ribosomal subunit protein L34
Reads supporting (aligned to +/- strand):  ref base G (22/21);  new base A (145/122);  total (168/144)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 313 total aligned reads displayed.

GCGAAAAGGTGCGGAAAAGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCA  >  NC_000913/3881793‑3881974
                                                                                            |                                                                                         
gCGAAAAGGTGCGGAAAAGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCGTCTTTAATCACCCATgg                                                                                           >  1:4777707/1‑93 (MQ=255)
  gAAAAGGTGCGGAAAAGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACt                                                                                         >  1:562139/1‑93 (MQ=255)
  gAAAAGGTGCGGAAAAGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACt                                                                                         <  1:2821855/93‑1 (MQ=255)
   aaaaGGTGCGGAAAAGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACtt                                                                                        >  1:5274279/1‑93 (MQ=255)
   aaaaGGTGCGGAAAAGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACtt                                                                                        >  1:800986/1‑93 (MQ=255)
    aaaGGTGCGGAAAAGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTc                                                                                       <  2:6545477/93‑1 (MQ=255)
      aGGTGCGGAAAAGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCgg                                                                                     >  1:6237936/1‑93 (MQ=255)
      aGGTGCGGAAAAGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCgg                                                                                     >  2:6900668/1‑93 (MQ=255)
       ggTGCGGAAAAGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGgt                                                                                    <  2:3258356/93‑1 (MQ=255)
         tGCGGAAAAGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGgtgt                                                                                  >  2:6045739/1‑93 (MQ=255)
             gAAAAGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCAt                                                                              >  1:3110334/1‑93 (MQ=255)
             gAAAAGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCAt                                                                              >  2:2169505/1‑93 (MQ=255)
             gAAAAGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCAt                                                                              >  2:6155349/1‑93 (MQ=255)
               aaaGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCg                                                                            >  2:3642491/1‑93 (MQ=255)
               aaaGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCg                                                                            >  2:2739336/1‑93 (MQ=255)
                aaGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGt                                                                           >  1:6947857/1‑93 (MQ=255)
                 aGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGtt                                                                          >  1:3935481/1‑93 (MQ=255)
                  gcgcGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGttt                                                                         >  1:4968777/1‑93 (MQ=255)
                  gcgcGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGttt                                                                         >  2:5444749/1‑93 (MQ=255)
                   cgcgGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTc                                                                        >  1:6893787/1‑93 (MQ=255)
                    gcgGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCa                                                                       <  2:5343293/93‑1 (MQ=255)
                    gcgGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCa                                                                       <  1:4568968/93‑1 (MQ=255)
                     cgGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCAt                                                                      <  2:7159274/93‑1 (MQ=255)
                      ggTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCAtt                                                                     >  2:1187749/1‑93 (MQ=255)
                       gTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAGTCACCCATGACTTCGGTGTCCATCGTTTCAttt                                                                    >  2:7260813/1‑93 (MQ=255)
                       gTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCAttt                                                                    >  1:5108823/1‑93 (MQ=255)
                        tAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCAtttt                                                                   <  2:5381606/93‑1 (MQ=255)
                        tAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCAtttt                                                                   <  1:1673219/93‑1 (MQ=255)
                          aaTAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTc                                                                 >  2:6284452/1‑93 (MQ=255)
                          aaTAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTc                                                                 >  2:3583739/1‑93 (MQ=255)
                          aaTAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTc                                                                 >  1:1203318/1‑93 (MQ=255)
                          aaTAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTc                                                                 >  2:7441155/1‑93 (MQ=255)
                           aTAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTcg                                                                >  1:1092180/1‑93 (MQ=255)
                            tAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTcga                                                               <  2:1437026/93‑2 (MQ=255)
                            tAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTcga                                                               >  2:753492/1‑92 (MQ=255)
                            tAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTcga                                                               >  1:3012318/1‑92 (MQ=255)
                             aaGGAAAGAGAATTGTCTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCTCCCATGGCTTCGGTGTCCATCGTTTCATTTTTcggc                                                              <  1:6436449/93‑1 (MQ=255)
                             aaGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTcgac                                                              <  2:2503868/93‑3 (MQ=255)
                              aGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCCCCCATGACTTCGGTGTCCATCGTTTCATTTTTcgacg                                                             <  1:4502154/93‑1 (MQ=255)
                              aGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTcgacg                                                             >  2:3360836/1‑93 (MQ=255)
                              aGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATACGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTcgacg                                                             >  1:2021409/1‑93 (MQ=255)
                               ggAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTcggcgg                                                            >  1:519129/1‑93 (MQ=255)
                                  aaGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGata                                                         >  2:6473116/1‑93 (MQ=255)
                                  aaGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGata                                                         <  1:2631479/93‑1 (MQ=255)
                                   agagAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGatat                                                        >  2:5600928/1‑93 (MQ=255)
                                   agagAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGatat                                                        <  2:6575284/93‑1 (MQ=255)
                                    gagaATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATc                                                       >  1:3713107/1‑93 (MQ=255)
                                       aaTTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCaa                                                    <  1:1423903/93‑1 (MQ=255)
                                        aTTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAAt                                                   <  1:1462107/93‑1 (MQ=255)
                                         ttGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATa                                                  >  2:3990379/1‑93 (MQ=255)
                                           gACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATaaa                                                >  2:4948301/1‑93 (MQ=255)
                                            aCTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAg                                               >  1:6572612/1‑93 (MQ=255)
                                             cTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGc                                              >  2:4440993/1‑93 (MQ=255)
                                             cTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGc                                              <  1:3210998/93‑1 (MQ=255)
                                              tCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGcc                                             >  1:6863934/1‑93 (MQ=255)
                                               ccGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCa                                            <  1:1005473/93‑1 (MQ=255)
                                                 ggAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCAtt                                          >  1:5448972/1‑93 (MQ=255)
                                                 ggAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCAtt                                          <  1:4727984/93‑1 (MQ=255)
                                                 ggAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCAtt                                          >  2:449717/1‑93 (MQ=255)
                                                   aGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGa                                        >  2:1958304/1‑93 (MQ=255)
                                                    gtgtACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGaa                                       <  1:288504/93‑1 (MQ=255)
                                                    gtgtACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGaa                                       <  2:4212818/93‑1 (MQ=255)
                                                     tgtACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAAt                                      >  1:3304035/1‑93 (MQ=255)
                                                     tgtACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAAt                                      >  2:4621719/1‑93 (MQ=255)
                                                     tgtACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAAt                                      >  2:2247603/1‑93 (MQ=255)
                                                      gtACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAAtt                                     <  2:4144486/93‑1 (MQ=255)
                                                       tACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAAttt                                    >  1:7318810/1‑93 (MQ=255)
                                                          aaTTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTAtt                                 >  2:685981/1‑93 (MQ=255)
                                                          aaTTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTAtt                                 <  2:5889372/93‑1 (MQ=255)
                                                          aaTTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTAtt                                 <  1:6486559/93‑1 (MQ=255)
                                                            ttattaCAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCa                               >  1:6174513/1‑93 (MQ=255)
                                                            ttattaCAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCa                               <  2:1081481/93‑1 (MQ=255)
                                                             tattaCAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCaa                              <  1:5328323/93‑1 (MQ=255)
                                                               ttaCAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGt                            >  1:2453761/1‑93 (MQ=255)
                                                               ttaCAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGt                            >  2:7384605/1‑93 (MQ=255)
                                                                taCAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGtt                           <  1:6006083/93‑1 (MQ=255)
                                                                taCAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGtt                           <  2:3867803/93‑1 (MQ=255)
                                                                taCAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGtt                           >  2:4064849/1‑93 (MQ=255)
                                                                  cAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTa                         <  2:4896198/93‑1 (MQ=255)
                                                                  cAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTa                         >  1:1152665/1‑93 (MQ=255)
                                                                    aTCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAgg                       >  1:5872019/1‑93 (MQ=255)
                                                                     tCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGt                      >  1:4026545/1‑93 (MQ=255)
                                                                     tCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGt                      <  2:2846380/93‑1 (MQ=255)
                                                                      ccGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTa                     <  2:3198555/93‑1 (MQ=255)
                                                                         gCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGaa                  >  1:568941/1‑93 (MQ=255)
                                                                         gCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGaa                  <  1:3517197/93‑1 (MQ=255)
                                                                         gCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGaa                  >  1:4331488/1‑93 (MQ=255)
                                                                          ccTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGaaa                 >  2:1694058/1‑93 (MQ=255)
                                                                           ctctTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAAt                >  1:5975080/1‑93 (MQ=255)
                                                                            tctTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATc               >  2:6399381/1‑93 (MQ=255)
                                                                            tctTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATc               <  2:2398814/93‑1 (MQ=255)
                                                                              tttAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGc             <  1:3057642/93‑1 (MQ=255)
                                                                                tAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCa           >  1:7103523/1‑93 (MQ=255)
                                                                                     aCCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGaaa      <  1:3419321/93‑1 (MQ=255)
                                                                                     aCCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGaaa      >  2:5220344/1‑93 (MQ=255)
                                                                                     aCCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGaaa      <  1:2227032/93‑1 (MQ=255)
                                                                                     aCCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGaaa      >  1:908898/1‑93 (MQ=255)
                                                                                      cccATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAAc     <  1:4165264/93‑1 (MQ=255)
                                                                                       ccATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACg    >  1:1859988/1‑93 (MQ=255)
                                                                                       ccATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACg    <  1:5578553/93‑1 (MQ=255)
                                                                                            |                                                                                         
GCGAAAAGGTGCGGAAAAGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCA  >  NC_000913/3881793‑3881974

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