Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,944,744 T→C 5.9% S205G (AGT→GGT)  ruvB ← ATP‑dependent DNA helicase, component of RuvABC resolvasome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,944,7440TC 5.9% 197.9 / 3.0 69S205G (AGT→GGT) ruvBATP‑dependent DNA helicase, component of RuvABC resolvasome
Reads supporting (aligned to +/- strand):  ref base T (34/30);  new base C (2/2);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

AACATATCCAGCGCCTGAGCAGCGATATCTGCCGAGATGGTGCCATCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATACCAAAACGGTCGCGCAACGGTGATGTCAGCGAACCTGCGCGCGTGGTTGCACCAATCAGG  >  NC_000913/1944596‑1944889
                                                                                                                                                    |                                                                                                                                                 
aaCATATCCAGCGCCTGAGCAGCGCTATCTGCCGAGATGGTGCCATCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGCCGGTTGGCACTGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACt                                                                                                                                                   >  2:387921/1‑149 (MQ=255)
  cATATCCAGCGCCTGAGCAGCGATATCTGCCGAGATGGTGCCATCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCCTCACCCa                                                                                                                                                 >  1:1071732/1‑149 (MQ=255)
  cATATCCAGCGCCTGAGCAGCGATATCTGCCGAGATGGTGCCATCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACCCa                                                                                                                                                 >  1:310022/1‑149 (MQ=255)
    tatCCAGCGCCTGAGCAGCGATATCTGCCGAGATGGTGCCATCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGTCGAGCAACTTCCAGCGCGCCGTCATCACTCAtc                                                                                                                                               <  1:614895/149‑1 (MQ=255)
     atCCAGCGCCTGAGCAGCGATATCTGCCGAGATGGTGCCATCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCAtct                                                                                                                                              >  1:342972/1‑149 (MQ=255)
     atCCAGCGCCTGAGCAGCGATATCTGCCGAGATGGTGCCATCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCAtct                                                                                                                                              >  1:452050/1‑149 (MQ=255)
      tCCAGCGCCTGAGCAGCGATATCTGCCGAGATGGTGCCATCGGGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCAtctc                                                                                                                                             <  2:224180/149‑1 (MQ=255)
          gcgcCTGAGCAGCGATATCTGCCGAGATGGTGCCATCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGc                                                                                                                                         <  1:837576/149‑1 (MQ=255)
             ccTGAGCAGCGATATCTGCCGAGATGGTGCCATCGTGCTTCACTTCGGCGAAATCACGCACTCGCCGCAGCAGACGGTTGGCAATGCGCGGCGAACCGCGAGCGCGACGAGAAACTTCAAGAGCGCCCTCATCACTCATCTCAAGccgc                                                                                                                                      >  2:842531/1‑147 (MQ=255)
               tGAGCAGCGATATCTGCCGAGATGGTGCCATCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCAt                                                                                                                                    <  2:1064387/149‑1 (MQ=255)
                gAGCAGCGATATCTGCCGAGATGGTGCCATCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATa                                                                                                                                   >  2:714369/1‑149 (MQ=255)
                 tgcagcGATATCTGCCGAGAAGGTGCCATCGTGCTTCACTTCTGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATaa                                                                                                                                  <  2:193619/148‑1 (MQ=255)
                    agcGATATCTGCCGAGATGGTGCCATCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGc                                                                                                                               >  1:709468/1‑149 (MQ=255)
                    agcGATATCTGCCGAGATGGTGCCATCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGc                                                                                                                               <  2:165194/149‑1 (MQ=255)
                      cgttatCTGCAGTGATGGTGCCCTCGTGCTTCACTTCGGCGTAATAACGCACTCGACGCTGCAGACGGTTGTCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGt                                                                                                                             <  2:474437/146‑1 (MQ=255)
                      cGATATCTGCCGAGATGGTGCCATCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCTCTCATCTCAAGCCCCATAAAGCGt                                                                                                                             <  1:464990/149‑1 (MQ=255)
                           tCTGCCGAGATGGTGCCATCGTGCTTCACTTCGGCGAAATCACGCTCTCGACGCAGCAGACGGTTGGCTATGCGCGGCGTACCGAGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTGAAGCCCCATAAAGCGTGCGCt                                                                                                                        <  2:900291/149‑1 (MQ=255)
                                cGAGATGGTGCCATCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGGTGGCAATGCGCGGCGGACCGCGCGCGCGCCGAGCAACTTCCAGCGCGACGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGAc                                                                                                                   >  1:573723/1‑149 (MQ=255)
                                    aTGGGGCCATCGTGCTTCACGTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGAc                                                                                                               <  2:1173637/149‑1 (MQ=255)
                                       gTGCCATCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGat                                                                                                            <  1:518373/149‑1 (MQ=255)
                                           cATCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCACGCCCCATAAAGCGTGCGCTGCGACTGACGctatat                                                                                                        >  1:437595/1‑149 (MQ=255)
                                           cATCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGatatat                                                                                                        <  1:1169398/149‑1 (MQ=255)
                                            aTCGTGCTTCACTTCGGCGAAATCACGCACTCGCCGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGTGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATAtt                                                                                                       >  2:549719/1‑149 (MQ=255)
                                            aTCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCCTCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATAtt                                                                                                       >  2:723568/1‑149 (MQ=255)
                                             tCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCAGCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTg                                                                                                      >  1:16726/1‑149 (MQ=255)
                                               gTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCa                                                                                                    <  1:575178/149‑1 (MQ=255)
                                                  cTTCACTTCGGCGAAATCACGCTCTCGACGCCGCAGCCGGTTGGCAATGCGCGGCGTACAGCGATCGCGACGAGCAACTTCCAGTGCGCCGTCATCACTCATCTCAAGCCACATAAAGCGTGCGCTGCGACTGACGATATATTGCAGAt                                                                                                 <  2:816469/149‑1 (MQ=255)
                                                   ttCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGAAGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCTAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATc                                                                                                <  2:9294/149‑1 (MQ=255)
                                                       cTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTGGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGc                                                                                            <  2:1037022/149‑1 (MQ=255)
                                                        ttCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCa                                                                                           <  2:852099/149‑1 (MQ=255)
                                                           ggCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGCCGATATATTGCAGATCCGGCACCt                                                                                        >  1:168172/1‑149 (MQ=255)
                                                           ggCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCt                                                                                        >  2:846590/1‑149 (MQ=255)
                                                              gAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCACGCCCCATAAAGCGTGCGCTGCGAGTGACGATATATTGCAGATCCGGCACCTGAt                                                                                     >  2:129381/1‑149 (MQ=255)
                                                                  tCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGCCTGACGATATATTGCAGATCCGGCACCTGATaaaa                                                                                 >  1:361480/1‑149 (MQ=255)
                                                                  tCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATaaaa                                                                                 >  1:998046/1‑149 (MQ=255)
                                                                    ccGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTACAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGATGCGACTGACGCTATATTGCAGATCCGGCACCTGATAAAAca                                                                               >  2:1049351/2‑148 (MQ=255)
                                                                    aCGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGCCGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGAGCGCTGCGACTGCCGATATATTGCAGATCCGGCACCTGATaaaccg                                                                               >  1:858700/1‑146 (MQ=255)
                                                                     cGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCCGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTc                                                                              >  2:116356/1‑149 (MQ=255)
                                                                      gCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGAGACTGACGATATGTTGCAGATCCGGCACCTGATAAAACTcc                                                                             >  1:541550/1‑149 (MQ=255)
                                                                        aCTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGCGCAACTTCCAGCGCGCCGTCATCACTCATCTCCAGCCGCATAAAGCGTGCGCTGCGCCTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAg                                                                           >  1:833687/1‑149 (MQ=255)
                                                                          tCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGAc                                                                         >  2:15913/1‑149 (MQ=255)
                                                                             tcGCAGCAGACGGTTGGCAATGCTCGGCGTACAGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAATCCCCATAAAGCGTTCGCTGCGTCTTACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGtt                                                                      <  1:284450/148‑1 (MQ=255)
                                                                             aCGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGtt                                                                      <  2:839788/149‑1 (MQ=255)
                                                                              cggagcagACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCTTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTg                                                                     <  1:652141/146‑1 (MQ=255)
                                                                               gcagccTACTGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGc                                                                    <  1:676580/149‑1 (MQ=255)
                                                                                 agcagACGGTTGGCAATGCGCGGCGTACCTCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGcac                                                                  <  1:714369/149‑1 (MQ=255)
                                                                                 agcagACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGCGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGAGATATTGCAGATCCGGCACCTGATAAACCTCCAGACGTTGccc                                                                  >  1:86707/1‑147 (MQ=255)
                                                                                   cagACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCATGCCCTGCGACTGCAGATATATTGCAGATACGGCACCTGATAACACTCCAGACGTTGCACaa                                                                >  2:690712/1‑149 (MQ=255)
                                                                                     gACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATa                                                              >  2:598154/1‑149 (MQ=255)
                                                                                      aCGGTTGGCAATGCGCGGCGTACCGCGCGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACAGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATaa                                                             >  1:118239/1‑148 (MQ=255)
                                                                                      aCGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATTTATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATAc                                                             >  1:19544/1‑149 (MQ=255)
                                                                                            ggCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCAACACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATACCAAAAc                                                       <  2:341478/149‑1 (MQ=255)
                                                                                              cAATGCGCGGCGTACCGCGAGCGCGCCGAGCAACTTCCAGCGCGCCGTCAACACACATCTCACGCCCCATACTGCGTGCGCTGCGACTGACGATATATTGACGATCCGGCACCGGATAAATCTCCAGCCGTTGCACAATACAAAAAAgg                                                     >  1:906596/1‑149 (MQ=255)
                                                                                                  gcgcGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATACCAAAACGGTcgc                                                 >  2:488589/1‑149 (MQ=255)
                                                                                                        cgTCCCGCGAGCGCGACGAGCAACTTCCAGGGCGCCGTCATCACTCATCTCAACCCCCATCTTGCGTGCGCTGCGACTGTCGATATATTGCAGATCCGGGACCTGATAAATCTGCAGACGTTGCTCAATACCAAAACGGTGGCGCAACg                                           <  2:1014002/149‑1 (MQ=255)
                                                                                                          tACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATACCAAAACGGTCGCGCAACGGt                                         <  2:1091998/149‑1 (MQ=255)
                                                                                                           aCCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATACCAAACCGGTCGCGCAACGGTg                                        >  2:1026262/1‑149 (MQ=255)
                                                                                                             cgcgAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGAGCCGGCACCTGATAAAACTCCAGACGTTGCACAATACCAAAACGGTCGCGCAACGGTGAt                                      >  1:487759/1‑149 (MQ=255)
                                                                                                               cgAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATACCAAAACGGTCGCGCAACGGTGATGt                                    >  1:582414/1‑149 (MQ=255)
                                                                                                               cgAGCGCGACCAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATACCAAAACGGTCGCGCAACGGTGATGt                                    >  1:547899/1‑149 (MQ=255)
                                                                                                                      gacgaGCAACTTCCAGCGCGCCGTCATCCCTCATCTCAAGCCCCATAAAGCGTGAGCTGCGACTGACGATATATTGAAGCTCCGGCACCTGATACAACTCCAGACGTGGCACAATCCCAAACCGGTCGCGCTACGGTGATGTCcgcgca                             >  1:385101/1‑147 (MQ=255)
                                                                                                                       actaGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGATGCTACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATACCAAAACGGTCGCGCAACGGTGATGTCAGCGAAc                            <  2:79067/146‑1 (MQ=255)
                                                                                                                             aaCTTCCAGCGCGCCGTCCTCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATACCAAAACGGTCGCGCAACGGTGATGTCAGCGAACCTgcgc                      <  2:1034006/149‑1 (MQ=255)
                                                                                                                               cTTCCACCGCGCCGGCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATACCAAAACGGTCGCGCAACGGTGATGTCAGCGAACCTgcgcgc                    <  2:385804/149‑1 (MQ=255)
                                                                                                                                     gcgcgcCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATACCAAAACGGTCGCGCAACGGTGATGTCAGCGAACCTGCGCGCGTGGtt              >  2:616876/1‑149 (MQ=255)
                                                                                                                                     gcgcgcCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATACCAAAACGGTCGCGCAACGGTGATGTCAGCGAACCTGCGCGCGTGGtt              <  1:387921/149‑1 (MQ=255)
                                                                                                                                      cgcgcCGTCATCACCCATCTCAAGCCCCGTAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATACCAAAACGGTCGCGCAACGGTGATGTCAGCGAACCTGCGCGCGTGGTTg             <  2:310022/149‑1 (MQ=255)
                                                                                                                                      cgcgcCGTCATCACCCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATACCAAAACGGTCGCGCAACGGTGATGTCAGCGAACCTGCGCGCGTGGTTg             <  2:1071732/149‑1 (MQ=255)
                                                                                                                                       gcgcCTTCATCACTCCTCTCAATCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATACCAAAACGGTCGCGCAACGGTGATGTCAGCGAACCTGCGCGCGTGGTTGc            <  2:573723/149‑1 (MQ=255)
                                                                                                                                                 caCTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTCATAAAACTCCAGACGTTGCACAATACCAAAACGGTCGCGCAACGGTGATGTCAGCGAACCTGCGCGCGTGGTTGCACCAATCAgg  >  1:977055/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AACATATCCAGCGCCTGAGCAGCGATATCTGCCGAGATGGTGCCATCGTGCTTCACTTCGGCGAAATCACGCACTCGACGCAGCAGACGGTTGGCAATGCGCGGCGTACCGCGAGCGCGACGAGCAACTTCCAGCGCGCCGTCATCACTCATCTCAAGCCCCATAAAGCGTGCGCTGCGACTGACGATATATTGCAGATCCGGCACCTGATAAAACTCCAGACGTTGCACAATACCAAAACGGTCGCGCAACGGTGATGTCAGCGAACCTGCGCGCGTGGTTGCACCAATCAGG  >  NC_000913/1944596‑1944889

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