Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,402,520 T→A 5.3% intergenic (+44/‑108) PP_3010 → / → PP_3011 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,402,5200TA 5.3% 228.6 / 2.6 75intergenic (+44/‑108)PP_3010/PP_3011hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/35);  new base A (2/2);  total (38/37)
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.81e-01

TTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTGGAAGCAAGAA  >  NC_002947/3402372‑3402664
                                                                                                                                                    |                                                                                                                                                      
ttGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCtt                                                                                                                                                        >  7:65584/1‑149 (MQ=255)
   cGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCttttt                                                                                                                                                     <  4:303803/149‑1 (MQ=255)
     gggCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTtgt                                                                                                                                                   >  2:225134/1‑149 (MQ=255)
        cACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACAAAAAGCCCGCTCTCGACGGCGGGCTATTggtgtg                                                                                                                                                <  6:173200/149‑1 (MQ=38)
           gggAAGCCAGCGCCACTTCACGGCCTATCATCGTGATGCCTGGTACGACGATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACAAAAAGCCCGCTCTCGACGGCGGGCTATTggtgtgggg                                                                                                                                             <  8:141164/149‑4 (MQ=37)
                     cgcCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGcc                                                                                                                                   <  7:174883/149‑1 (MQ=255)
                      gcCACTTCACGGCCTATCAACGTGATGGCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGccc                                                                                                                                  >  1:197671/1‑149 (MQ=255)
                        cACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCg‑‑g                                                                                                                              <  1:225134/149‑1 (MQ=255)
                              aCGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCa                                                                                                                            >  8:321353/1‑145 (MQ=255)
                              aCGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCa                                                                                                                            <  7:321353/145‑1 (MQ=255)
                               cGGCCTATCAACGTGATGCCTGGTACGACCATGGCATGTTGCTGGATGCTGCACTGATGTGTTAAAAAGCGTAATCATCCAACCCCCTCACCTATAGCCCGGCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGtt                                                                                                                      <  7:175630/149‑1 (MQ=255)
                                 gCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTc                                                                                                                     >  6:183489/1‑148 (MQ=255)
                                  ccTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCaa                                                                                                                   <  8:65584/149‑1 (MQ=255)
                                      tCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCttttt                                                                                                                                                     >  5:17108/1‑114 (MQ=255)
                                      tCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCttttt                                                                                                                                                     <  6:17108/114‑1 (MQ=255)
                                      tCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCaaaaaa                                                                                                               <  5:183489/149‑1 (MQ=255)
                                              aTGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCggg                                                                                                       <  2:207484/149‑1 (MQ=255)
                                                  cTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGc                                                                                                                                               <  6:325839/108‑1 (MQ=255)
                                                  cTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGc                                                                                                                                               >  5:325839/1‑108 (MQ=255)
                                                    ggTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGc                                                                                                 >  7:115487/1‑149 (MQ=255)
                                                        cGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGt                                                                                             <  2:197671/149‑1 (MQ=255)
                                                         gACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTa                                                                                           <  3:186651/149‑1 (MQ=255)
                                                          aCCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCgg                                                                                                        >  5:220060/1‑136 (MQ=255)
                                                          aCCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCgg                                                                                                        <  6:220060/136‑1 (MQ=255)
                                                                    gTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGagca                                                                                <  3:299517/149‑1 (MQ=255)
                                                                       gCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGagcagc                                                                              <  2:3549/148‑1 (MQ=255)
                                                                        cTGGATGCTGCACGGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGt                                                                            >  4:321284/1‑149 (MQ=255)
                                                                                tgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTg                                                                    >  2:271475/1‑148 (MQ=255)
                                                                                tgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCAGCGCCTGCCAG‑‑GCAACGTTCAAAAAAATGGCGGGTGGAGCCGGTCTGGCCTGAGCATCTTTGAtttt                                                                    >  4:198047/1‑147 (MQ=255)
                                                                                     tGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTgg                                                                                                    <  2:288977/113‑1 (MQ=255)
                                                                                     tGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTgg                                                                                                    >  1:288977/1‑113 (MQ=255)
                                                                                     tGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTg                                                             >  7:49833/1‑149 (MQ=255)
                                                                                       aTGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGcc                                                           >  5:134226/1‑149 (MQ=255)
                                                                                         gtgtTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCgg                                                         >  3:225343/1‑149 (MQ=255)
                                                                                            ttCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGccc                                                      >  7:229451/1‑149 (MQ=255)
                                                                                             tCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGcccc                                                     >  1:188210/1‑149 (MQ=255)
                                                                                              cAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGCGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGTCGCGTCGAGCCGGTATGGCCTGAGCAGCGTTGCTTTTCATTTGCCGGAcccg                                                    >  2:358344/1‑149 (MQ=255)
                                                                                                 aaGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCcccgccc                                                 >  7:187921/1‑149 (MQ=255)
                                                                                                 aaGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCcccgccc                                                 >  3:283530/1‑149 (MQ=255)
                                                                                                    cGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTg                                                                                                     >  1:299415/1‑97 (MQ=255)
                                                                                                    cGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTg                                                                                                     <  2:299415/97‑1 (MQ=255)
                                                                                                       aaTCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGAGGAGGCGTTCAAAAAAAGGGCCGGTTGAACCCGTAATGGCATAACAGCCTTTCATTTGAATTTCCCGGCcccgccgt                                                >  1:111309/1‑146 (MQ=255)
                                                                                                           atcCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCg                                       <  7:200672/149‑1 (MQ=255)
                                                                                                           atcCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCg                                       >  8:200672/1‑149 (MQ=255)
                                                                                                           atcCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCg                                       >  5:70493/1‑149 (MQ=255)
                                                                                                              cAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGtt                                     <  4:283530/148‑1 (MQ=255)
                                                                                                              cAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTg                                    <  6:134226/149‑1 (MQ=255)
                                                                                                                 tgggcacacaAAAAGCCCGCTCTCGACGGCGGGCTATTGGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTAt                                 >  7:304810/5‑149 (MQ=37)
                                                                                                                   cccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTAttt                               >  1:76482/1‑149 (MQ=255)
                                                                                                                    gcacacaAAAAGCCCGCTCTCGACGGCGGGCTATTGGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtg                              >  2:231826/2‑149 (MQ=38)
                                                                                                                    ccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTAt                                                                                          >  4:107283/1‑91 (MQ=255)
                                                                                                                    ccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTAt                                                                                          <  3:107283/91‑1 (MQ=255)
                                                                                                                    ccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtg                              >  8:240809/1‑149 (MQ=255)
                                                                                                                    ccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTCGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtg                              >  2:285965/1‑149 (MQ=255)
                                                                                                                      acacCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGc                                                                                       >  5:331044/1‑92 (MQ=255)
                                                                                                                      acacCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGc                                                                                       <  6:331044/92‑1 (MQ=255)
                                                                                                                      acacCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtgga                            >  6:34259/1‑149 (MQ=255)
                                                                                                                       cacCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggat                           <  8:49833/149‑1 (MQ=255)
                                                                                                                       cacCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggat                           <  4:225343/149‑1 (MQ=255)
                                                                                                                        acaaaaaGCCCGCTCTCGACGGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggatg                          >  1:292371/7‑149 (MQ=255)
                                                                                                                        acCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggatg                          <  5:172121/149‑1 (MQ=255)
                                                                                                                        acCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggatg                          >  2:29401/1‑149 (MQ=255)
                                                                                                                         ccAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggatgg                         <  8:229451/149‑1 (MQ=255)
                                                                                                                          cAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTa                                              <  3:101556/127‑1 (MQ=255)
                                                                                                                          cAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTa                                              >  4:101556/1‑127 (MQ=255)
                                                                                                                          cAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggatgga                        <  5:34259/149‑1 (MQ=255)
                                                                                                                          cAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggatgga                        <  3:321284/149‑1 (MQ=255)
                                                                                                                          aaaaaGCCCGCTCTCGACGGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggatg                          >  8:271567/5‑147 (MQ=255)
                                                                                                                               gCCCTCCGTCGAGAGCTGGCTTTTTGTGTGCCAATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGTGTGGAGCCGGTATGGACTGAGCAGCGTTGCTTTGCATTTGCCGTCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGt                    <  6:151606/148‑1 (MQ=255)
                                                                                                                                     cgTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGt             <  1:271475/149‑1 (MQ=255)
                                                                                                                                     agTCGAGAGCGGGCTTGTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAACAAGGGCGTGTGGAGCCGGTATGGCCGGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGt             <  3:198047/148‑1 (MQ=255)
                                                                                                                                      gTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATAATCGTTGTATTTGGATGGAGAGTCGCAGGTg            <  1:29401/149‑1 (MQ=255)
                                                                                                                                       tCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTg            >  3:48250/1‑148 (MQ=255)
                                                                                                                                         gagaGCGGGCTTTTTGTGTGCGCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggatgga                        <  6:139903/134‑1 (MQ=255)
                                                                                                                                         gagaGCGGGCTTTTTGTGTGCGCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggatgga                        >  5:139903/1‑134 (MQ=255)
                                                                                                                                             gCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTGGAAGCaa     <  2:292371/149‑1 (MQ=255)
                                                                                                                                             gCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGGGGAAGCaa     <  6:70493/149‑1 (MQ=255)
                                                                                                                                                ggCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTAttt                               >  2:205853/1‑120 (MQ=255)
                                                                                                                                                ggCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTAttt                               <  1:205853/120‑1 (MQ=255)
                                                                                                                                                 ggtttttGTGTGCCCATCGCCTGCCGG‑‑GCAAAGTTAAAAAAAAGGGCGGGTGGAGCCGGTTTGGCCTGATGAGCGGTGCTTTGCCCTTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTGGAAGCaagaa  <  1:285965/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCG‑‑GCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTGGAAGCAAGAA  >  NC_002947/3402372‑3402664

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: