Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,921,585 G→T 10.3% R665L (CGG→CTG)  PP_4330 → toxin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,921,5850GT10.3% 192.9 / 10.6 68R665L (CGG→CTG) PP_4330toxin‑like protein
Reads supporting (aligned to +/- strand):  ref base G (32/29);  new base T (3/4);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

GTCGGCGACAAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGC  >  NC_002947/4921437‑4921733
                                                                                                                                                    |                                                                                                                                                    
gTCGGCGACAAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCg                                                                                                                                                      <  8:138564/149‑1 (MQ=255)
    gCGACAAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAATGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACACAAGGCTATTGGCAGCCCGGGcc                                                                                                                                                  >  1:222534/1‑149 (MQ=255)
      gACAAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGccc                                                                                                                                                 >  4:167908/1‑148 (MQ=255)
         aaaGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAgg                                                                                                                                             >  5:192076/1‑149 (MQ=255)
         aaaGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAgg                                                                                                                                             >  1:49612/1‑149 (MQ=255)
          aaGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAggg                                                                                                                                            <  7:192960/149‑1 (MQ=255)
             tCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCCCGGGCTg                                                                                                                                         <  1:9142/149‑1 (MQ=255)
               aTTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGc                                                                                                                                        <  6:23639/148‑1 (MQ=255)
                ttAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAg                                                                                                                                      <  1:22321/149‑1 (MQ=255)
                ttAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAg                                                                                                                                      <  5:85793/149‑1 (MQ=255)
                   gCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGAGAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTTCGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGt                                                                                                                                    <  4:3297/148‑1 (MQ=255)
                   gCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCACTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTa                                                                                                                                   <  4:98641/149‑1 (MQ=255)
                    cGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCCGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCCCGGGCTGccaata                                                                                                                                   <  3:195491/148‑6 (MQ=255)
                    cGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCCGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCCCGGGCTGccaa                                                                                                                                     >  4:195491/1‑143 (MQ=255)
                       gCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGGTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGc                                                                                                                               >  4:43703/1‑149 (MQ=255)
                       gCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCTCGGGCTGc                                                                                                                                        <  1:76740/140‑1 (MQ=255)
                       gCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCTCGGGCTGc                                                                                                                                        >  2:76740/1‑140 (MQ=255)
                          aTTGCTGGCCAACCTTGCCGCCAACGCCGCGCGACTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTAGGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGc                                                                                                                            <  2:102901/149‑1 (MQ=255)
                           ttGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCt                                                                                                                           <  1:65860/149‑1 (MQ=255)
                             gCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCtt                                                                                                                          >  2:153752/1‑148 (MQ=255)
                             gCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCtt                                                                                                                          <  3:167908/148‑1 (MQ=255)
                               tGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGcccc                                                                                                                                                <  3:1978/124‑1 (MQ=255)
                               tGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGcccc                                                                                                                                                >  4:1978/1‑124 (MQ=255)
                                      ccTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAg                                                                                                                <  7:176905/149‑1 (MQ=255)
                                            cgccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCACTTGAAGAGCAAGGcttct                                                                                                          <  6:192076/149‑1 (MQ=255)
                                              ccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGcttctt                                                                                                         >  2:172387/1‑148 (MQ=255)
                                              ccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCttcttc                                                                                                        >  3:270315/1‑149 (MQ=255)
                                              ccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGCTATTGCCGCTTGAAGAGCAAGGcttctt                                                                                                         <  1:172387/148‑1 (MQ=255)
                                                 aggccccGCGCCTGCCCTACAGTGAACACCGCATGCGCCGGTTTCTGTCGAAGGCTGCGCCTACGGCTTACGACCAGACCCCAAGGCTATTGGCAGCCCGGGCCCCTGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCgg                                                                                                     <  6:184287/147‑1 (MQ=255)
                                                   gccgcGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAaga                                                                                                                     >  2:8728/1‑131 (MQ=255)
                                                   gccgcGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAaga                                                                                                                     <  1:8728/131‑1 (MQ=255)
                                                        gcgcCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCt                                                                                              <  7:32973/149‑1 (MQ=255)
                                                          gcCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGc                                                                                            <  5:277941/149‑1 (MQ=255)
                                                           ccTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGcc                                                                                           <  1:153752/149‑1 (MQ=255)
                                                                   aCAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGc                                                                                   <  4:270315/149‑1 (MQ=255)
                                                                     aGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGTTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTa                                                                                 <  1:213735/149‑1 (MQ=255)
                                                                     aGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGc                                                                                   >  6:86433/1‑147 (MQ=255)
                                                                     aGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGc                                                                                   <  5:86433/147‑1 (MQ=255)
                                                                       tGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTAt                                                                                >  5:132525/1‑148 (MQ=255)
                                                                           tACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGC                                                                           >  3:91491‑M1/1‑149 (MQ=255)
                                                                                cATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGCTCCGCTATTACGCGCAG                                                                       >  3:185551‑M1/1‑148 (MQ=255)
                                                                                cATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAG                                                                       >  7:253410‑M1/1‑148 (MQ=255)
                                                                                   gcgcCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCGAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTAc                                                                    >  7:258576/1‑148 (MQ=255)
                                                                                        tGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTCGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTCAAGAGCAAGGCTTCTTCCGGCCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCa                                                              >  5:30639/1‑149 (MQ=255)
                                                                                                   aaGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATgg                                                   >  3:76288/1‑149 (MQ=255)
                                                                                                       cTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGcacatcggtcgcagacgcgtaatt                                               >  8:87097‑M1/1‑125 (MQ=255)
                                                                                                       cTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAAc                                               >  2:167260/1‑149 (MQ=255)
                                                                                                        ttCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTCCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACg                                              >  2:141710/1‑149 (MQ=255)
                                                                                                          tttttACCTCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCATGGTTGCAGGTATTGCCGTTTGAAGAGCAAGGCTTCTTCCGTGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGt                                            <  7:111290/145‑1 (MQ=255)
                                                                                                           gCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTg                                           >  3:241799/1‑149 (MQ=255)
                                                                                                              tACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCa                                        >  3:298453/1‑149 (MQ=255)
                                                                                                               aCCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCATCGCTGTATGGCAACTGttata                                        >  8:111290/1‑144 (MQ=255)
                                                                                                                  gCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGc                                    <  1:167260/149‑1 (MQ=255)
                                                                                                                       cGAACAGTACGCAAGGCTATTGGGAGCCCGGTCGCCAGGGCTGCAGGTATTGGCGCTTTAAGAGCATGGCTTCTTGCGGGCCTTCTGCCAGAGCCGCTATTACGCGGAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATGGccacc                                <  6:30639/148‑1 (MQ=255)
                                                                                                                        gAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCGATGccaccgt                              <  2:49612/149‑3 (MQ=255)
                                                                                                                            aGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCG                                                                          >  5:259281‑M1/1‑101 (MQ=255)
                                                                                                                            aGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCG                                                                          <  6:259281‑M1/101‑1 (MQ=255)
                                                                                                                                   aaGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTg                   >  4:210179/1‑149 (MQ=255)
                                                                                                                                     ggCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGcacatcggtcgcagacgcgtaattggtggcgatgccaccgttgccatccagcgc                 >  5:258070‑M1/1‑95 (MQ=255)
                                                                                                                                     ggCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGc                  >  7:11899/1‑148 (MQ=255)
                                                                                                                                       atacctgCAGCCCTGGGGCCCGGGCTGCCGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCt                    >  2:126748/7‑144 (MQ=255)
                                                                                                                                       atacctgCAGCCCTGGGGCCCGGGCTGCCGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCt                    <  1:126748/138‑1 (MQ=255)
                                                                                                                                            ggCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAt           <  4:241799/148‑1 (MQ=255)
                                                                                                                                             gCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAtgt         >  3:248353/1‑149 (MQ=255)
                                                                                                                                             gCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAtgt         >  8:11911/1‑149 (MQ=255)
                                                                                                                                              cAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAtgtg        <  4:248353/149‑1 (MQ=255)
                                                                                                                                               aGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGc       >  4:91491‑M2/1‑148 (MQ=255)
                                                                                                                                                 cccGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGctg     >  4:76288‑M2/1‑146 (MQ=255)
                                                                                                                                                  ccGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGctgc    <  8:78322‑M2/149‑5 (MQ=255)
                                                                                                                                                    gggCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGcacatcggtcgcagacgcgtaattggtggcgatgccaccgttgccatccagcgctgcctcgtactgcgc  <  6:258070‑M1/149‑70 (MQ=255)
                                                                                                                                                    gggCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGcacatcggtcgcagacgcgtaattggtggcgatgccaccgttgccatccagcgctgcctcgtactgcgc  <  7:87097‑M1/149‑70 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTCGGCGACAAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGC  >  NC_002947/4921437‑4921733

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