Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,859,400 G→C 7.8% H3Q (CAC→CAG)  PP_1662 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,859,4000GC 7.8% 267.6 / 9.7 89H3Q (CAC→CAG) PP_1662hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (44/37);  new base C (3/4);  total (47/42)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

CCCACCAGTTGACCTTGTTGGTACACATCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGG  >  NC_002947/1859260‑1859544
                                                                                                                                            |                                                                                                                                                 
cccACCAGTTGACCTTGTTGGTACACATCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACACACTGTCGCag                                                                                                                                           <  1:162179/149‑2 (MQ=255)
 ccaccaGTTGACCTTGTTGGTACACATCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACg                                                                                                                                          >  1:24412/1‑149 (MQ=255)
  caccaGTTGACCTTGTTGGTACACATCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACACACTGTCGCagc                                                                                                                                          >  2:246690/1‑146 (MQ=255)
  caccaGTTGACCTTGTTGGTACACATCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTACAAGATGGCTGCGACACACTGTCGCagc                                                                                                                                          <  1:246690/148‑3 (MQ=255)
  caccaGTTGACCTTGTTGGTACACATCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAATCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGt                                                                                                                                         <  4:160795/149‑1 (MQ=255)
       gTTGACCTTGTTGGTACACATCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGtcttct                                                                                                                                    <  3:249449/149‑1 (MQ=255)
        ttGACCTTGTTGGTACACATCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGTACATTCTCGATTCGTCGCGTTCGGTTACTTTGTTAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTcttctt                                                                                                                                   <  1:341963/149‑1 (MQ=255)
        ttGACCTTGTTGGTACACATCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTcttctt                                                                                                                                   <  3:145734/149‑1 (MQ=255)
           aCCTTGTTGGTACACATCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTta                                                                                                                                 >  2:267052/1‑148 (MQ=255)
               tgttgGTACACATCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATagag                                                                                                                            >  5:70844/1‑149 (MQ=255)
                gttgGTACACATCGCAGATCCGCAGGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTATACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATaga                                                                                                                             <  3:306404/147‑1 (MQ=255)
                    gTACACATCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTg                                                                                                                        <  8:188486/148‑1 (MQ=255)
                    gTACACATCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGt                                                                                                                       <  7:39372/149‑1 (MQ=255)
                    gTACACATCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGt                                                                                                                       <  8:317808/149‑1 (MQ=255)
                     tACACATCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGtt                                                                                                                      >  3:1637/1‑149 (MQ=255)
                      acacATCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCAGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGt                                                                                                                       >  4:198237/1‑147 (MQ=255)
                        acaTGCCAGATCCGCACGACAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGGGTTCGGTTACTTGGACAGTCATATGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGc                                                                                                                   <  5:410136/149‑1 (MQ=255)
                        acaTCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAg                                                                                                                    >  8:103648/1‑148 (MQ=255)
                        acaTCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAg                                                                                                                    >  3:206867/1‑148 (MQ=255)
                        acaTCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAg                                                                                                                    >  2:396940/1‑148 (MQ=255)
                        acaTCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGc                                                                                                                   <  4:206867/149‑1 (MQ=255)
                            cccAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGAc                                                                                                               >  7:91032/1‑149 (MQ=255)
                             ccAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGAc                                                                                                               >  7:248623/1‑148 (MQ=255)
                             ccAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGAc                                                                                                               >  8:313576/1‑148 (MQ=255)
                              cAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACt                                                                                                              >  5:86806/1‑148 (MQ=255)
                               aGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACtaa                                                                                                            <  5:209655/149‑1 (MQ=255)
                                    cGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACtaagtaa                                                                                                        <  8:91032/148‑1 (MQ=255)
                                       aCGCCAGAGGCCAGGACTTCGTGGGAAATTCTCGATTCGTCGCGTTCGGTTAATTGGACAGTCCTCGGTAAGCTCAACTTCTTCAAGATGGATGCGACAGTGTGTCGCTCGTCTTCTTTATAGAGTTTGTGAGCAGACTAactaagta                                                                                                     >  6:254133/1‑148 (MQ=255)
                                          ccAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTGCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAAc                                                                                                 <  6:86806/149‑1 (MQ=255)
                                           cAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTcttctt                                                                                                                                   >  5:413459/1‑114 (MQ=255)
                                           cAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTcttctt                                                                                                                                   <  6:413459/114‑1 (MQ=255)
                                           cAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTaaaa                                                                                                  >  6:192443/1‑147 (MQ=255)
                                              cGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCAGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGc                                                                                             <  3:198237/149‑1 (MQ=255)
                                                gCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCAt                                                                                           >  1:203080/1‑149 (MQ=255)
                                                     gCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTGGGACAGGCATCGGTAAACTCCACTTCTTCAAGCTGGCTGCGAGAGTGTGGCGCACGTCTTCTTTATAGAGTTTGTTAGCAGa                                                                                                                <  7:126945/123‑1 (MQ=255)
                                                        ttCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAg                                                                                   >  7:300345/1‑149 (MQ=255)
                                                              ggACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCAt                                                                             >  2:403246/1‑149 (MQ=255)
                                                                aCATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATg                                                                            >  4:161597/1‑148 (MQ=255)
                                                                    tctcGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCt                                                                       >  3:59146/1‑149 (MQ=255)
                                                                     ctcGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTc                                                                      >  3:332133/1‑149 (MQ=255)
                                                                      tcGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTc                                                                      >  8:37941/1‑148 (MQ=255)
                                                                      tcGATTCGTCGCGTTCGGTTACTTGGACACTCAACGGTAAACTCCACTTCTTCAAGATCGCTGCGACAGTGTGACGCACGTCATCTTTATAGAGTATCTTACCAGACTAAGTAAGTAAAAAAGGCATGAAAAAAGAATCATCAAGcac                                                                      >  2:92629/1‑146 (MQ=255)
                                                                       cGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCg                                                                     >  6:373069/1‑148 (MQ=255)
                                                                        gATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGaa                                                                   <  6:88978/149‑1 (MQ=255)
                                                                        gATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGaa                                                                   <  6:316032/149‑1 (MQ=255)
                                                                         aTTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAg                                                                  >  5:100801/1‑149 (MQ=255)
                                                                         aTTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAg                                                                  >  8:63198/1‑149 (MQ=255)
                                                                             gtcgCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GAt                                                             >  6:311887/1‑149 (MQ=255)
                                                                                      ggTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCtgg                                                    <  8:248623/149‑1 (MQ=255)
                                                                                      ggTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCtgg                                                    <  6:70844/149‑1 (MQ=255)
                                                                                       gTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCtggg                                                   >  2:251981/1‑149 (MQ=255)
                                                                                       gTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCtggg                                                   >  7:335377/1‑149 (MQ=255)
                                                                                        ttACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGCTGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGGAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCtgggt                                                  <  4:1637/149‑1 (MQ=255)
                                                                                        ttACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCtgggt                                                  >  6:275925/1‑149 (MQ=255)
                                                                                           cTTGGAGAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCtgggtgg                                                <  4:59146/148‑1 (MQ=255)
                                                                                           cTTGGACAGTCATCGGTAATCTCCACTTCTTCAAGATTGCTTCGACAGTTTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCtgggtgg                                                <  5:144474/148‑1 (MQ=255)
                                                                                           cTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCtgggtggg                                               >  5:217244/1‑149 (MQ=255)
                                                                                           cTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCtgggtgg                                                <  5:319747/148‑1 (MQ=255)
                                                                                            ttGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGc                                              >  1:203529/1‑149 (MQ=255)
                                                                                              ggACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCtt                                            >  2:223108/1‑149 (MQ=255)
                                                                                                  aGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTc                                        <  2:24412/149‑1 (MQ=255)
                                                                                                  aGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTc                                        >  3:43987/1‑149 (MQ=255)
                                                                                                     cATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGATCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCg                                       >  6:40166/1‑147 (MQ=255)
                                                                                                      aTCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTTGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCgca                                     <  1:358479/148‑1 (MQ=255)
                                                                                                      aTCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCttt                                           <  7:292108/142‑1 (MQ=255)
                                                                                                      aTCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCttt                                           >  8:292108/1‑142 (MQ=255)
                                                                                                      aTCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCgca                                     <  3:161597/148‑1 (MQ=255)
                                                                                                      aTCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCgca                                     <  1:251981/148‑1 (MQ=255)
                                                                                                          ttAAACTCAACTTCTTCAAGATGGCTGCTACATTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGAGATGCTGCCTGGGTGGGCTTTTTCGCAGCaa                                 <  1:267052/148‑1 (MQ=255)
                                                                                                          gTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCgcagca                                  >  4:186803/1‑147 (MQ=255)
                                                                                                          gTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCgcagca                                  >  8:193864/1‑147 (MQ=255)
                                                                                                              aCTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGaa                             >  1:258199/1‑148 (MQ=255)
                                                                                                                tCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTc                                        <  3:138632/135‑1 (MQ=255)
                                                                                                                tCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTc                                        >  4:138632/1‑135 (MQ=255)
                                                                                                                  cACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAg                          <  4:43987/147‑1 (MQ=255)
                                                                                                                          cAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCa                   <  8:300345/146‑1 (MQ=255)
                                                                                                                               tGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCgg           <  7:313576/149‑1 (MQ=255)
                                                                                                                                 gCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGTTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCgggg         <  7:103648/149‑1 (MQ=255)
                                                                                                                                 gCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCgggg         <  2:203080/149‑1 (MQ=255)
                                                                                                                                 cgtGCGACACACTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCgggg         >  4:181842/3‑149 (MQ=255)
                                                                                                                                  gtGCGACACACTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGc        <  3:181842/148‑1 (MQ=255)
                                                                                                                                  cTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGc        >  2:249840/1‑149 (MQ=255)
                                                                                                                                  cTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGc        >  5:369274/1‑149 (MQ=255)
                                                                                                                                  cTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCACGCACtgggctgggc        >  2:256437/1‑149 (MQ=255)
                                                                                                                                   tGCGACACACTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCg            <  6:95713/144‑1 (MQ=255)
                                                                                                                                   tGCGACACACTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCg            >  5:95713/1‑144 (MQ=255)
                                                                                                                                       aCAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCCCAGTGCGGGGCttttc   <  1:396940/149‑4 (MQ=255)
                                                                                                                                        cAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAgg  >  1:201472/1‑149 (MQ=255)
                                                                                                                                        cAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAgg  <  6:100801/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
CCCACCAGTTGACCTTGTTGGTACACATCCCAGATCCGCACGCCAGCGGCCAGGCCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGG‑GATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGG  >  NC_002947/1859260‑1859544

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: ‑