Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,881,913 G→A intergenic (‑556/‑52) 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,9130GA88.4% 892.2 / 66.1 301intergenic (‑556/‑52)dnaA/rpmHDnaA/50S ribosomal subunit protein L34
Reads supporting (aligned to +/- strand):  ref base G (16/19);  new base A (132/134);  total (148/153)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-03

Only 100 of 311 total aligned reads displayed.

AAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCG  >  NC_000913/3881822‑3882005
                                                                                           |                                                                                            
 aGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCCCCCATGACTTCGGTGTCCATCGTTTCATTTTTcgacg                                                                                            <  1:4502154/93‑1 (MQ=255)
  ggAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTcggcgg                                                                                           >  1:519129/1‑93 (MQ=255)
   gAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGa                                                                                          <  1:4697322/93‑1 (MQ=255)
    aaaGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGat                                                                                         >  2:837288/1‑93 (MQ=255)
    aaaGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGat                                                                                         <  2:1136568/93‑1 (MQ=255)
     aaGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGata                                                                                        <  1:817266/93‑1 (MQ=255)
     aaGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGata                                                                                        >  2:336360/1‑93 (MQ=255)
      agagAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGatat                                                                                       >  2:6480921/1‑93 (MQ=255)
       gagaATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATc                                                                                      >  1:2989874/1‑93 (MQ=255)
         gaATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCa                                                                                    >  2:6438297/1‑93 (MQ=255)
          aaTTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCaa                                                                                   >  1:5998396/1‑93 (MQ=255)
          aaTTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCaa                                                                                   >  1:5492274/1‑93 (MQ=255)
           aTTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAAt                                                                                  >  2:5369485/1‑93 (MQ=255)
            ttGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATa                                                                                 <  1:3575498/93‑1 (MQ=255)
              gACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATaaa                                                                               <  2:3641390/93‑1 (MQ=255)
              gACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATaaa                                                                               <  1:5296443/93‑1 (MQ=255)
              gACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAACCCCCCAGGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATaaa                                                                               <  1:1161845/93‑1 (MQ=255)
                cTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGc                                                                             >  2:4440993/1‑93 (MQ=255)
                  ccGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCCCCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCa                                                                           <  1:3195966/93‑1 (MQ=255)
                  ccGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCa                                                                           <  1:109772/93‑1 (MQ=255)
                  ccGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCa                                                                           <  1:1005473/93‑1 (MQ=255)
                    ggAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCAtt                                                                         >  2:449717/1‑93 (MQ=255)
                    ggAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCAtt                                                                         <  1:7089010/93‑1 (MQ=255)
                     gAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTg                                                                        >  1:6904864/1‑93 (MQ=255)
                      aGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGa                                                                       <  1:6478826/93‑1 (MQ=255)
                       gtgtACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGaa                                                                      <  2:4212818/93‑1 (MQ=255)
                       gtgtACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGaa                                                                      >  2:3797541/1‑93 (MQ=255)
                        tgtACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAAt                                                                     >  1:3304035/1‑93 (MQ=255)
                         gtACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAAtt                                                                    <  2:4144486/93‑1 (MQ=255)
                         gtACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAAtt                                                                    <  2:7211513/93‑1 (MQ=255)
                         gtACAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAAtt                                                                    >  1:197611/1‑93 (MQ=255)
                          tACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAAttt                                                                   >  1:6729068/1‑93 (MQ=255)
                           aCAATTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTa                                                                  <  2:5488024/93‑1 (MQ=255)
                             aaTTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTAtt                                                                <  2:5889372/93‑1 (MQ=255)
                             aaTTATTACAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTAtt                                                                >  2:3948492/1‑93 (MQ=255)
                               ttattaCAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCa                                                              <  2:2433479/93‑1 (MQ=255)
                               ttattaCAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCa                                                              >  1:5383889/1‑93 (MQ=255)
                                 attaCAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAg                                                            >  2:2620636/1‑93 (MQ=255)
                                  ttaCAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGt                                                           <  1:5635677/93‑1 (MQ=255)
                                   taCAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGtt                                                          <  2:4054060/93‑1 (MQ=255)
                                     cAATCCGGCCTCTTTAATCACCCATGCCTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTa                                                        <  1:6183135/93‑1 (MQ=255)
                                     cAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTa                                                        <  2:4896198/93‑1 (MQ=255)
                                     cAATCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTa                                                        <  1:4901335/93‑1 (MQ=255)
                                       aTCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAgg                                                      >  1:5872019/1‑93 (MQ=255)
                                       aTCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAgg                                                      <  2:7313444/93‑1 (MQ=255)
                                        tCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGt                                                     >  1:1615067/1‑93 (MQ=255)
                                        tCCGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGt                                                     >  1:4026545/1‑93 (MQ=255)
                                         ccGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTa                                                    <  2:3198555/93‑1 (MQ=255)
                                         ccGGCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTa                                                    <  1:3901940/93‑1 (MQ=255)
                                            gCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGaa                                                 >  1:568941/1‑93 (MQ=255)
                                            gCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGaa                                                 <  1:2479405/93‑1 (MQ=255)
                                            gCCTCTTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGaa                                                 >  2:2172397/1‑93 (MQ=255)
                                              ctctTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAAt                                               >  1:5975080/1‑93 (MQ=255)
                                               tctTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATc                                              <  1:2860610/93‑1 (MQ=255)
                                                ctTTAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCg                                             >  1:145066/1‑93 (MQ=255)
                                                 tttAATCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGc                                            <  1:3057642/93‑1 (MQ=255)
                                                    aaTCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCAt                                         <  2:2304888/93‑1 (MQ=255)
                                                    aaTCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCAt                                         <  1:3095900/93‑1 (MQ=255)
                                                     aTCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATg                                        >  1:1370262/1‑93 (MQ=255)
                                                     aTCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATg                                        >  1:6570689/1‑93 (MQ=255)
                                                      tCACCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGa                                       <  1:3433674/93‑1 (MQ=255)
                                                        aCCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGaaa                                     >  2:5220344/1‑93 (MQ=255)
                                                        aCCCATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGaaa                                     <  2:6861632/93‑1 (MQ=255)
                                                         cccATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAAc                                    <  1:4165264/93‑1 (MQ=255)
                                                         cccATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAAc                                    >  1:2416252/1‑93 (MQ=255)
                                                          ccATGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACg                                   >  2:7437378/1‑93 (MQ=255)
                                                            aTGACTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCa                                 <  2:3095352/93‑1 (MQ=255)
                                                             tgacTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCAc                                >  1:1723852/4‑93 (MQ=255)
                                                              gacTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACt                               <  1:7417884/91‑1 (MQ=255)
                                                              gacTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACt                               >  2:4727361/3‑93 (MQ=255)
                                                               acTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACtt                              >  1:3826483/2‑93 (MQ=255)
                                                                cTTCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACttt                             >  2:5650795/1‑93 (MQ=255)
                                                                 ttCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACtttt                            <  2:205295/93‑1 (MQ=255)
                                                                  tCGGTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTc                           <  1:4329278/93‑1 (MQ=255)
                                                                    ggTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCaa                         >  1:4415384/1‑93 (MQ=255)
                                                                    ggTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCaa                         <  2:2890434/93‑1 (MQ=255)
                                                                    ggTGTCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCaa                         <  2:5701109/93‑1 (MQ=255)
                                                                     gtgtCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAAc                        <  2:171280/93‑1 (MQ=255)
                                                                      tgtCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAAcc                       <  1:4817586/93‑1 (MQ=255)
                                                                      tgtCCATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAAcc                       <  2:6689931/93‑1 (MQ=255)
                                                                         ccATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTc                    >  1:6342751/1‑93 (MQ=255)
                                                                         ccATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTc                    <  1:963246/93‑1 (MQ=255)
                                                                         ccATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTc                    <  2:5970296/93‑1 (MQ=255)
                                                                         ccATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTc                    <  2:4638971/93‑1 (MQ=255)
                                                                          cATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCt                   >  1:6064894/1‑93 (MQ=255)
                                                                          cATCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCt                   >  2:3271120/1‑93 (MQ=255)
                                                                           aTCGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTg                  <  1:2168435/93‑1 (MQ=255)
                                                                             cGTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTa                <  2:626220/93‑1 (MQ=255)
                                                                              gTTTCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTAc               <  2:3583553/93‑1 (MQ=255)
                                                                               tttCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACt              >  1:7341899/1‑93 (MQ=255)
                                                                               tttCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACt              >  1:7296601/1‑93 (MQ=255)
                                                                                ttCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTg             >  2:1353506/1‑93 (MQ=255)
                                                                                ttCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTg             >  2:618699/1‑93 (MQ=255)
                                                                                ttCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTg             >  1:2059891/1‑93 (MQ=255)
                                                                                 tCATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGa            <  2:2533617/93‑1 (MQ=255)
                                                                                  cATTTTTCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGaa           <  2:4569207/93‑1 (MQ=255)
                                                                                      tttCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAgcgc       <  2:1588618/93‑1 (MQ=255)
                                                                                      tttCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAgcgc       <  2:2714490/93‑1 (MQ=255)
                                                                                       ttCGACGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCa      >  2:4283763/1‑93 (MQ=255)
                                                                                           gcggATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCg  >  1:2656627/1‑93 (MQ=255)
                                                                                           |                                                                                            
AAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCG  >  NC_000913/3881822‑3882005

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