Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 552,371 G→A 6.1% P75L (CCG→CTG)  purK ← N5‑carboxyaminoimidazole ribonucleotide synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913552,3710GA 6.1% 207.5 / 2.5 65P75L (CCG→CTG) purKN5‑carboxyaminoimidazole ribonucleotide synthase
Reads supporting (aligned to +/- strand):  ref base G (32/29);  new base A (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

CCAGTGCGACGCTTAACAATCGCCAGCTCACCTAAACGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGTCCAGCCCGACTGGCCAGACAG  >  NC_000913/552224‑552518
                                                                                                                                                   |                                                                                                                                                   
ccaGTGCGACGCTTAACAATCGCCAGCTCACCTAAACGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCgg                                                                                                                                                    <  6:101857/149‑1 (MQ=255)
ccaGTGCGACGCTTAACAATCGCCAGCTCACCTAAACGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCgg                                                                                                                                                    >  7:363450/1‑149 (MQ=255)
   gTGCGACGCTTAACAATCGCCAGCTCACCTAAACGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGaa                                                                                                                                                 <  7:53845/149‑1 (MQ=255)
       gACGCTTAACAATCGCCAGCTCACCTAAACGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAAcaca                                                                                                                                             >  5:16362/1‑149 (MQ=255)
        aCGCTTAACAATCGCCAGCTCACCTAAACGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACAt                                                                                                                                            <  4:136084/149‑1 (MQ=255)
         cGCTTAACAATCGCCAGCTCACCTAAACGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATc                                                                                                                                           >  6:370929/1‑149 (MQ=255)
           cTTAACAATCGCCAGCTCACCTAAACGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATcgc                                                                                                                                         >  3:210498/1‑149 (MQ=255)
            ttAACAATCGCCAGCTCACCTAAACGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATcgcg                                                                                                                                        >  5:88158/1‑149 (MQ=255)
              aaCAATCGCCAGCTCACCTAAACGATC‑AACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGt                                                                                                                                      <  8:324234/148‑1 (MQ=255)
               aCAATCGCCAGCTCACCTAAACGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGtt                                                                                                                                     >  7:32950/1‑149 (MQ=255)
                   tCGCCAGCTCACCTAAACGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTcaca                                                                                                                                 <  4:237403/149‑1 (MQ=255)
                        aGCTCACCTAAACGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGc                                                                                                                            <  2:118790/149‑1 (MQ=255)
                         gCTCACCTAAACGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGcc                                                                                                                           <  5:34587/149‑1 (MQ=255)
                             aCCTAAACGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGa                                                                                                                        <  6:357427/148‑1 (MQ=255)
                               cTAAACGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATgg                                                                                                                     >  7:273232/1‑149 (MQ=255)
                                  aaCGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGcc                                                                                                                           >  3:241692/1‑140 (MQ=255)
                                  aaCGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGcc                                                                                                                           <  4:241692/140‑1 (MQ=255)
                                  aaCGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGgcgc                                                                                                                  <  8:363450/149‑1 (MQ=255)
                                  aaCGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGgcgc                                                                                                                  >  5:115052/1‑149 (MQ=255)
                                  aaCGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGgcgc                                                                                                                  <  8:216277/149‑1 (MQ=255)
                                     gATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGgcgcgc                                                                                                                >  1:59506/1‑148 (MQ=255)
                                       tCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGAt                                                                                                                       <  3:7371/139‑1 (MQ=255)
                                       tCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGAt                                                                                                                       >  4:7371/1‑139 (MQ=255)
                                            cacCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTcgc                                                                                                        <  2:59506/149‑1 (MQ=255)
                                                gCAGGCCACTCGCTGCGTTCGGCAAGTCACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGt                                                                                                    <  3:408030/149‑1 (MQ=255)
                                                      cACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGCCGGTCAGGAATAATCGGGAGCACATCGCGGTTCGCATAGGCCGGATGGCGCGGCAGCTGGGGGGTTAATGc                                                                                              >  8:362068/1‑149 (MQ=255)
                                                       aCTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCg                                                                                             >  6:58228/1‑149 (MQ=255)
                                                         tCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGt                                                                                           >  4:299164/1‑149 (MQ=255)
                                                           gCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGtt                                                                                                                                     <  2:295421/105‑1 (MQ=255)
                                                           gCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGtt                                                                                                                                     >  1:295421/1‑105 (MQ=255)
                                                           gCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGtt                                                                                          <  4:210498/148‑1 (MQ=255)
                                                            cTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTc                                                                                        <  7:359106/149‑1 (MQ=255)
                                                                gTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGCAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGc                                                                                    >  4:374189/1‑149 (MQ=255)
                                                                 ttCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGc                                                                                    >  5:29615/1‑148 (MQ=255)
                                                                   cGGCAAGTAACTGCCACGGTGCAGTCGTCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACCTCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAg                                                                                 >  5:74018/1‑149 (MQ=255)
                                                                     gCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAaca                                                                                                                                               >  6:368160/1‑85 (MQ=255)
                                                                     gCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAaca                                                                                                                                               <  5:368160/85‑1 (MQ=255)
                                                                      cAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAATAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAAGACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGt                                                                              <  1:191417/149‑1 (MQ=255)
                                                                            aCTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGCAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTAtc                                                                        <  3:374189/149‑1 (MQ=255)
                                                                             cTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTAtct                                                                       >  7:14545/1‑149 (MQ=255)
                                                                                      tGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCTGCAATAATCGGGAACACATCGCGGTTCACAAAGGCGGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTaa                                                              <  6:88158/149‑1 (MQ=255)
                                                                                      tGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCGGCGTTCTATCTCAGCGGTaa                                                              <  8:32950/149‑1 (MQ=255)
                                                                                      tGCAGGCGGCAGGTGGAGCTTATCGAACAGCTGCTTCTGAGTCAGACGGGCTGCAATAATCGGGAACACATCGCGGGTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGGTTCCGGCCAGCGTTCTATCTCAGCGGTaa                                                              <  5:234541/149‑1 (MQ=255)
                                                                                        cAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGgcgcgc                                                                                                                <  6:233327/97‑1 (MQ=255)
                                                                                        cAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGgcgcgc                                                                                                                >  5:233327/1‑97 (MQ=255)
                                                                                               cAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCttt                                                     <  8:273232/149‑1 (MQ=255)
                                                                                                  gTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGACCACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTttgtt                                                  >  5:103988/1‑149 (MQ=255)
                                                                                                   tGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTttgtt                                                  >  4:69448/1‑148 (MQ=255)
                                                                                                           aTCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGc                                          <  6:16362/148‑1 (MQ=255)
                                                                                                           aTCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCa                                         >  5:157658/1‑149 (MQ=255)
                                                                                                             cGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACTTCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCAcc                                       >  5:287198/1‑149 (MQ=255)
                                                                                                                 aaGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCAccgccg                                   >  4:245675/1‑149 (MQ=255)
                                                                                                                             aGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGt                       >  3:153219/1‑149 (MQ=255)
                                                                                                                             aGTCAGACGGTCAGCAATAATCAGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTaa                                                                                                 <  8:14384/75‑1 (MQ=255)
                                                                                                                             aGTCAGACGGTCAGCAATAATCAGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTaa                                                                                                 <  8:7935/75‑1 (MQ=255)
                                                                                                                             aGTCAGACGGTCAGCAATAATCAGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTaa                                                                                                 >  7:14384/1‑75 (MQ=255)
                                                                                                                             aGTCAGACGGTCAGCAATAATCAGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTaa                                                                                                 >  7:7935/1‑75 (MQ=255)
                                                                                                                                 aGACGGTCAGCAATAATCGTGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCCGCGTTCTATCTCAGCGGTAAGCCCGCTTTGGTGAAAAGGCACCGCCGCCGGTTCAGCGTCCAg                   >  5:23445/1‑149 (MQ=255)
                                                                                                                                 aGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGTCCa                    <  6:115052/148‑1 (MQ=255)
                                                                                                                                  gACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGTCCCGc                  >  3:131218/1‑149 (MQ=255)
                                                                                                                                       tCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGTCCAGCCCGAc             <  3:299164/149‑1 (MQ=255)
                                                                                                                                         aGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGTCCAGCCCGACTg           >  2:299404/1‑149 (MQ=255)
                                                                                                                                             ataatCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCTGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGTCCAGCCCGACTGGCCa       <  8:35439/149‑1 (MQ=255)
                                                                                                                                             ataatCGGGAACACATAGCGGTTCACAAAGGCCGGATAGCGCGCCACCTCGCGGGTTAATGCGGTTTCCGGCAAGCATTCTATCTCAGCGGTAATCACGCTTTCTTGAAAAGGAACAGCCGACGGTTCAGCGTCCAGCCCGACTGGc         >  3:395671/1‑147 (MQ=255)
                                                                                                                                              taatCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGTCCAGCCCGACTGGc         >  5:353623/1‑146 (MQ=255)
                                                                                                                                                atCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGTCCAGCCCGACTGGCCAGa     <  5:370929/148‑1 (MQ=255)
                                                                                                                                                 tCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGTCCAGCCCGACTGGCCAGa     >  7:298944/1‑147 (MQ=255)
                                                                                                                                                   gggAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGTCCAGCCCGACTGGCCAGACAg  >  3:98486/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCAGTGCGACGCTTAACAATCGCCAGCTCACCTAAACGATCAAACACCGCAGGCCACTCGCTGCGTTCGGCAAGTAACTGCCACGGTGCAGTCGGCAGGTGGAGCTTATCGAAAAGCTGCTTCTGAGTCAGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGTCCAGCCCGACTGGCCAGACAG  >  NC_000913/552224‑552518

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