Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,301,577 G→C 8.1% intergenic (+67/‑1233) PP_3774 → / → PP_5627 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,301,5770GC 8.1% 188.3 / 6.2 62intergenic (+67/‑1233)PP_3774/PP_5627hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base C (2/3);  total (30/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 = 5.02e-01

CCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCT  >  NC_002947/4301429‑4301717
                                                                                                                                                    |                                                                                                                                               
ccAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGcgcg                                                                                                                                                 >  3:46638/1‑149 (MQ=255)
 cAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCgg                                                                                                                                                >  7:233674/1‑149 (MQ=255)
   ggTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCgggg                                                                                                                                              >  7:157273/1‑149 (MQ=255)
       gACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGgggctgg                                                                                                                                          <  3:184978/149‑1 (MQ=255)
       gACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGgggctgg                                                                                                                                          >  7:81183/1‑149 (MQ=255)
             acCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCt                                                                                                                                    <  4:46638/149‑1 (MQ=255)
              ccGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCCCCGCGCTGGGCttttt                                                                                                                                   <  6:271725/149‑4 (MQ=255)
               cGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCCCCGCGCTGGGCtttttg                                                                                                                                  <  5:187890/149‑5 (MQ=255)
                gACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATc                                                                                                                                 <  5:137110/149‑1 (MQ=255)
                 aCGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCCCCGCGCTGGGCtttttgca                                                                                                                                <  4:135314/149‑7 (MQ=255)
                       cccccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTggg                                                                                                                          >  2:111593/1‑149 (MQ=255)
                       cccccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTggg                                                                                                                          >  3:182301/1‑149 (MQ=255)
                        ccccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCtttttgca                                                                                                                                <  3:334310/142‑7 (MQ=255)
                        ccccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCtttttgca                                                                                                                                >  4:334310/1‑136 (MQ=255)
                        ccccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGc                                                                                                                         >  4:229509/1‑149 (MQ=255)
                         cccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCt                                                                                                                        >  3:186785/1‑149 (MQ=255)
                         cccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCt                                                                                                                        <  4:186785/149‑1 (MQ=255)
                         cccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCt                                                                                                                        <  8:233674/149‑1 (MQ=255)
                          ccgccgAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGgggctggg                                                                                                                                         <  3:122338/131‑1 (MQ=255)
                          ccgccgAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGgggctggg                                                                                                                                         >  4:122338/1‑131 (MQ=255)
                          ccgccgAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTg                                                                                                                            >  8:235158/1‑144 (MQ=255)
                          ccgccgAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTg                                                                                                                            <  7:235158/144‑1 (MQ=255)
                                  gACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTAc                                                                                                              >  2:254385/1‑149 (MQ=255)
                                     ggCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTc                                                                                                           >  2:202516/1‑149 (MQ=255)
                                      gCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTcc                                                                                                          <  2:302926/149‑1 (MQ=255)
                                       cTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCg                                                                                                         >  8:4539/1‑149 (MQ=255)
                                          gATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGc                                                                                                        >  5:187174/1‑147 (MQ=255)
                                          gATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGc                                                                                                        <  6:187174/147‑1 (MQ=255)
                                             cTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGcc                                                                                                   >  7:188333/1‑149 (MQ=255)
                                                 atcatcAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTa                                                                                               <  7:4539/149‑1 (MQ=255)
                                                   catcaGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTaca                                                                                             >  5:334137/1‑149 (MQ=255)
                                                      caGCGAGCTCTCGCTCAACCATCGCTGAACGTCGATGCCCAGCACTGCTATCAAAACCACTCAGGGCGTCGCATGAATGCAAAAAGCCAAGCGCGGGGCTCGGCTTCTATCAGCTGGGCTCAACAGTACTTCCGCCAGCCTTTacaaaa                                                                                          >  8:135568/1‑146 (MQ=255)
                                                          gAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTa                                                                                      <  3:227979/149‑1 (MQ=255)
                                                            gCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACt                                                                                    <  8:188333/149‑1 (MQ=255)
                                                                  gCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGt                                                                            >  4:19812/1‑149 (MQ=255)
                                                                       tccatcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTc                                                                       >  8:14515/1‑149 (MQ=255)
                                                                       tccatcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTc                                                                       <  7:100344/149‑1 (MQ=255)
                                                                       tccatcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTc                                                                       <  7:14515/149‑1 (MQ=255)
                                                                       tccatcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTc                                                                       <  8:81183/149‑1 (MQ=255)
                                                                        ccatcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTc                                                                       >  7:224128/1‑148 (MQ=255)
                                                                         catcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCa                                                                     <  3:229509/149‑1 (MQ=255)
                                                                          atcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAg                                                                    <  4:182301/149‑1 (MQ=255)
                                                                           tcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGt                                                                   <  1:202516/149‑1 (MQ=255)
                                                                               tGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGgc                                                               >  1:35568/1‑149 (MQ=255)
                                                                               tGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGgc                                                               >  2:26500/1‑149 (MQ=255)
                                                                                           cccAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATaa                                                    <  8:157273/148‑1 (MQ=255)
                                                                                            ccAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTc                                                                       >  1:291798/1‑128 (MQ=255)
                                                                                            ccAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTc                                                                       <  2:291798/128‑1 (MQ=255)
                                                                                               gCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACtata                                               <  3:19812/149‑1 (MQ=255)
                                                                                                cccTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATaa                                              <  1:26500/149‑1 (MQ=255)
                                                                                                   tGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTtgt                                                         <  4:51471/135‑1 (MQ=255)
                                                                                                   tGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTtgt                                                         >  3:51471/1‑135 (MQ=255)
                                                                                                           aaaCCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTCTGATAACTATAACGCTCGGCAtt                                   >  3:300924/1‑149 (MQ=255)
                                                                                                            aaCCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTgg                                                                 >  4:356831/1‑118 (MQ=255)
                                                                                                            aaCCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTgg                                                                 <  3:356831/118‑1 (MQ=255)
                                                                                                                 cTCAGGTCGTCGCATGAATGGAAAAAGCCCAGGGCGGGGCTTGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCGTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGGTTGTGATAACTATAACGCTAGGCATTGGATGc                             <  4:300924/149‑1 (MQ=255)
                                                                                                                       gCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGa                       >  8:178601/1‑149 (MQ=255)
                                                                                                                             cATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCCCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTc                 <  1:254385/149‑1 (MQ=255)
                                                                                                                                 aaTGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCcggc             <  6:334137/149‑1 (MQ=255)
                                                                                                                                    gatagaaGCCCAGCCCCGCGCTGGGCATCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCggg          >  1:10945/6‑149 (MQ=255)
                                                                                                                                     cAAAAAGCCCAGCGCGGGGCTGCGCTTCTATCAGCTGCGCTGAACAGTACAGACGCCAGCCTTTACACACTGTACT‑‑GGAAGGTACACCACTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGc         >  1:89821/1‑149 (MQ=37)
                                                                                                                                       agaaGCCCAGCCCCGCGCTGGGCTTTTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCt       >  1:204199/3‑149 (MQ=255)
                                                                                                                                        aaaaCCCCGAGCGGGGGTGGGGTTATATAACCTGGGGTGAAAAGATAATACCCCCCCCTCTTACCCACTTTACTAATGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGc         <  1:111593/149‑1 (MQ=255)
                                                                                                                                        aaaGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCt       >  3:200508/1‑148 (MQ=255)
                                                                                                                                            cccAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGAtct  >  3:108970/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCT  >  NC_002947/4301429‑4301717

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: