Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,203,359 C→T 5.8% A37V (GCG→GTG)  PP_2810 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,203,3590CT 5.8% 208.5 / 2.3 69A37V (GCG→GTG) PP_2810hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/31);  new base T (2/2);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

GAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCAC  >  NC_002947/3203212‑3203504
                                                                                                                                                   |                                                                                                                                                  
gAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTgcg                                                                                                                                                   >  6:20682/1‑149 (MQ=255)
      ccGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGATCGGTACAGGCTGCGCCCCAg                                                                                                                                             <  5:224542/149‑1 (MQ=255)
        ggTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAggg                                                                                                                                           <  2:134741/149‑1 (MQ=255)
             aaGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCgg                                                                                                                                      <  1:277270/149‑1 (MQ=255)
              aGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGTGCCCCAGGGGGCg                                                                                                                                       >  6:338032/1‑147 (MQ=255)
               gCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAg                                                                                                                                    <  2:61102/149‑1 (MQ=255)
                cccTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTCCGCCCCCTGGGGCGCAgc                                                                                                                                   <  7:127883/149‑2 (MQ=255)
                   ttCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGgcgc                                                                                                                                >  8:14467/1‑149 (MQ=255)
                      gTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGTTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTgg                                                                                                                             <  6:51845/149‑1 (MQ=255)
                      gTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTgg                                                                                                                             >  3:234187/1‑149 (MQ=255)
                           aGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTc                                                                                                                         <  2:222292/148‑1 (MQ=255)
                                 ctccAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGtt                                                                                                                   <  1:24211/148‑1 (MQ=255)
                                 ctccAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGtt                                                                                                                   <  8:218686/148‑1 (MQ=255)
                                  tccAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGa                                                                                                                 >  2:349893/1‑149 (MQ=255)
                                        gAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGAcgccgc                                                                                                           >  8:50158/1‑149 (MQ=255)
                                             aaGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTgggg                                                                                                      >  4:288310/1‑149 (MQ=255)
                                               gTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGa                                                                                                                                     <  5:187912/116‑1 (MQ=255)
                                               gTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGa                                                                                                                                     >  6:187912/1‑116 (MQ=255)
                                                  aCTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGc                                                                                                   >  3:128705/1‑147 (MQ=255)
                                                   cTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCgag                                                                                                <  4:234187/149‑1 (MQ=255)
                                                      cACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGTGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCgaga                                                                                               <  5:338032/147‑1 (MQ=255)
                                                       aCGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAgcg                                                                                            >  4:359817/1‑149 (MQ=255)
                                                         ggTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTgg                                                                                                                             <  5:264430/114‑1 (MQ=255)
                                                         ggTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTgg                                                                                                                             >  6:264430/1‑114 (MQ=255)
                                                             ggAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGcgccgc                                                                                        >  5:297091/1‑147 (MQ=255)
                                                               aaaGTGTCGTATCCGCCCGCCATATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGt                                                                                    >  3:270775/1‑149 (MQ=255)
                                                                aaGTGTCGTATCCGCCCGCCATATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGt                                                                                    <  4:270775/148‑1 (MQ=255)
                                                                     tCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAgc                                                                              >  2:66686/1‑149 (MQ=255)
                                                                     tCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAgc                                                                              <  6:39848/149‑1 (MQ=255)
                                                                     tCGTATCCGCCAGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGcgccgc                                                                                        <  6:200901/139‑1 (MQ=255)
                                                                     tCGTATCAGCCCGCCGTATTTTCAGGCAGCTGCTGGTTAACCTTGCGGCCTCGACAGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGcgccgc                                                                                        >  5:200901/1‑139 (MQ=255)
                                                                         aTCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAgcggc                                                                           >  8:90633/1‑148 (MQ=255)
                                                                          tCCGCCCCCCGCAATTTCACGCAGCTGCTGGGTAACCATGCGGCCTCGCCCGCCCTCTGTACCGTACACCCTGCGCCCCACGGGGCGGAGGCGCTGGCTTCGAAGCTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGa                                                                         <  2:48133/149‑1 (MQ=255)
                                                                             gcccgccGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGatat                                                                      >  5:299249/1‑149 (MQ=255)
                                                                                     tATTTTCAGGCAGCTGCTGGTTATCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCt                                                              >  1:258895/1‑149 (MQ=255)
                                                                                       ttttCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGGTATTCCAGCCTgg                                                            >  3:114634/1‑149 (MQ=255)
                                                                                       ttttCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTg                                                             <  3:359817/148‑1 (MQ=255)
                                                                                       ttttCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTg                                                             <  6:297091/148‑1 (MQ=255)
                                                                                          tCAGGCAGCTGCTGGTTATCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGc                                                          <  2:258895/148‑1 (MQ=255)
                                                                                            aGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGcaca                                                       <  5:20682/149‑1 (MQ=255)
                                                                                              gCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGc                                                     >  8:120132/1‑149 (MQ=255)
                                                                                              gCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGc                                                     >  5:15744/1‑149 (MQ=255)
                                                                                              gCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGc                                                     >  3:216741/1‑149 (MQ=255)
                                                                                                aGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGc                                                     >  8:127839/1‑147 (MQ=255)
                                                                                                     ctgGTTACCCTTTCGTCCTCGACCGCTATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGGTCGAAGTTGACGCCGCTGGGGGGCTAGAGCGCCGCCAGTTCTAGCGGCTATATTACAGCCTGGTCACAGCCAGGgc                                               <  8:270114/148‑1 (MQ=255)
                                                                                                     ctgGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAggg                                                <  7:14467/147‑1 (MQ=255)
                                                                                                     ctgGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAggg                                                <  1:66686/147‑1 (MQ=255)
                                                                                                             ccTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTACTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑ggggg                                     >  1:124881/1‑149 (MQ=255)
                                                                                                               ttGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑gggggg                                    >  2:287409/1‑147 (MQ=255)
                                                                                                               ttGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑ggggg                                     >  3:290768/1‑147 (MQ=255)
                                                                                                                tGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAg                                  >  8:229014/1‑149 (MQ=255)
                                                                                                                 gCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAg                                  >  4:189503/1‑148 (MQ=255)
                                                                                                                  cGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGa                                <  7:127839/149‑1 (MQ=255)
                                                                                                                     ccTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGATg                              >  6:199837/1‑148 (MQ=255)
                                                                                                                     ccTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGATg                              >  2:71282/1‑148 (MQ=255)
                                                                                                                       tCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGGCAGGATGGa                            <  1:287409/149‑1 (MQ=255)
                                                                                                                       tCGACCGCGATCTGTTC‑GTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGATGGAGt                          >  7:279112/1‑149 (MQ=255)
                                                                                                                         gACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGATGGAGt                          <  7:90633/148‑1 (MQ=255)
                                                                                                                            cgcgATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGATGGAGTTACg                      <  7:229014/149‑1 (MQ=255)
                                                                                                                            cgcgATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGATGGAGTTACg                      <  2:354550/149‑1 (MQ=255)
                                                                                                                                    gTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGATGGa                            <  2:20670/135‑1 (MQ=255)
                                                                                                                                    gTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGATGGa                            >  1:20670/1‑135 (MQ=255)
                                                                                                                                        ggTACAGGCTGCGCCCGAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGATGGAGTTACGGCCAGGTACGAg          <  4:216741/149‑1 (MQ=255)
                                                                                                                                        ggTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGATGGAGTTACGGCCAGGTACGAg          <  3:288310/149‑1 (MQ=255)
                                                                                                                                         gTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGATGGAGTTACGGCCAGGTACGAgg         >  3:239718/1‑149 (MQ=255)
                                                                                                                                              ggCTGTGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTgg                                                            <  3:218555/94‑1 (MQ=255)
                                                                                                                                              ggCTGTGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTgg                                                            >  4:218555/1‑94 (MQ=255)
                                                                                                                                              gcctccgcCCCCTGGGGCGGAGGCGCTGGGTTCGAAGTTGCCGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAg     >  8:262720/6‑148 (MQ=255)
                                                                                                                                               gCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCa   >  6:13710/1‑149 (MQ=255)
                                                                                                                                               gCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCa   >  2:35057/1‑149 (MQ=255)
                                                                                                                                               gCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCa   >  1:161370/1‑149 (MQ=255)
                                                                                                                                                ctccgcCCCCTGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCAc  >  1:52581/4‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCA‑GGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCAC  >  NC_002947/3203212‑3203504

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