Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 537,464 T→A 7.2% noncoding (57/66 nt) PP_mr07 → P26 ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947537,4640TA 7.2% 311.4 / 11.7 97noncoding (57/66 nt)PP_mr07P26 ncRNA
Reads supporting (aligned to +/- strand):  ref base T (45/45);  new base A (3/4);  total (48/49)
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.78e-01

GGTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAG  >  NC_002947/537318‑537612
                                                                                                                                                  |                                                                                                                                                    
ggTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGcc                                                                                                                                                   >  2:2109919/1‑150 (MQ=255)
     tAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCgg                                                                                                                                              >  1:343333/1‑150 (MQ=255)
       gCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGACACCGGcc                                                                                                                                            >  1:88592/1‑150 (MQ=255)
             ggCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGTGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTc                                                                                                                                      <  1:659697/150‑1 (MQ=255)
             ggCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTc                                                                                                                                      <  1:1701912/150‑1 (MQ=255)
                gTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCAc                                                                                                                                                 >  2:1227437/1‑136 (MQ=255)
                gTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCAc                                                                                                                                                 <  1:1227437/136‑1 (MQ=255)
                gTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGt                                                                                                                                   <  1:1794862/150‑1 (MQ=255)
                gTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGt                                                                                                                                   >  1:464955/1‑150 (MQ=255)
                  ttCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTa                                                                                                                                 <  1:1809218/150‑1 (MQ=255)
                     aaaaGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGCCGCTAAAGTATAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtg                                                                                                                              <  2:1174343/149‑1 (MQ=255)
                      aaaGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtgg                                                                                                                             >  1:1597556/1‑150 (MQ=255)
                      aaaGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtgg                                                                                                                             >  1:1880106/1‑150 (MQ=255)
                       aagaagCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtggt                                                                                                                            <  2:343333/150‑1 (MQ=255)
                         gaagCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtggtgg                                                                                                                          <  1:2107483/150‑1 (MQ=255)
                            gCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCg                                                                                                                       >  1:2181536/1‑150 (MQ=255)
                            gCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCg                                                                                                                       <  1:565519/150‑1 (MQ=255)
                            gCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCg                                                                                                                       >  1:2375654/1‑150 (MQ=255)
                               gCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGc                                                                                                                                                    <  1:2374742/118‑1 (MQ=255)
                               gCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGc                                                                                                                                                    >  2:2374742/1‑118 (MQ=255)
                                     gACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCg                                                                                                              >  1:2021998/1‑150 (MQ=255)
                                     gACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTACTGCGGTCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCg                                                                                                              <  1:1233842/150‑1 (MQ=255)
                                          tAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTccc                                                                                                         <  2:956943/150‑1 (MQ=255)
                                          tAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCACGACCGTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTccc                                                                                                         >  1:1869891/1‑150 (MQ=255)
                                                  agCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGa                                                                                                 >  2:1793115/1‑150 (MQ=255)
                                                   gCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAg                                                                                                >  1:1197664/1‑150 (MQ=255)
                                                   gCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAg                                                                                                <  2:1197664/150‑1 (MQ=255)
                                                    cTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGc                                                                                               >  1:2303550/1‑150 (MQ=255)
                                                     tGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtggt                                                                                                                            <  2:123442/120‑1 (MQ=255)
                                                     tGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtggt                                                                                                                            >  1:123442/1‑120 (MQ=255)
                                                      ggAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTg                                                                                             >  1:831419/1‑150 (MQ=255)
                                                       gaagaaGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCCGGTGGCATAAGCCACCAGcc                                                                                                                                            >  1:883185/1‑102 (MQ=255)
                                                       gaagaaGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCCGGTGGCATAAGCCACCAGcc                                                                                                                                            <  2:883185/102‑1 (MQ=255)
                                                           aagCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGa                                                                                        <  2:2021998/150‑1 (MQ=255)
                                                             gCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGaca                                                                                      >  2:2089790/1‑150 (MQ=255)
                                                             gCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCCGGTGGCATAAGCCAccag                                                                                                                                              >  1:1845946/1‑92 (MQ=255)
                                                             gCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCCGGTGGCATAAGCCAccag                                                                                                                                              <  2:1845946/94‑3 (MQ=255)
                                                                ggCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACAccc                                                                                   <  2:2375654/150‑1 (MQ=255)
                                                                     tAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGt                                                                                                                     <  2:958281/111‑1 (MQ=255)
                                                                     tAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGt                                                                                                                     >  1:958281/1‑111 (MQ=255)
                                                                     tAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCAtgtg                                                                              >  2:751663/1‑150 (MQ=255)
                                                                      aaaGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATgtgg                                                                             <  1:626253/150‑1 (MQ=255)
                                                                         gTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATACCGAGCTGCAAGACACCCATgtgggt                                                                           <  1:667712/149‑1 (MQ=255)
                                                                         gTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATgtgggt                                                                           >  2:667712/1‑149 (MQ=255)
                                                                          tCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATgtgggtgg                                                                         <  2:88592/150‑1 (MQ=255)
                                                                           cGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGgc                                                                        <  1:1724155/150‑1 (MQ=255)
                                                                              gCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCACatgggtgtcttg                                                                     <  1:1113135/150‑6 (MQ=255)
                                                                                tGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATaa                                                                                                    <  2:225223/117‑1 (MQ=255)
                                                                                tGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATaa                                                                                                    >  1:225223/1‑117 (MQ=255)
                                                                                  aaGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTcc                                                                                                                                     >  1:681912/1‑82 (MQ=255)
                                                                                  aaGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTcc                                                                                                                                     <  2:681912/82‑1 (MQ=255)
                                                                                    gTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCaa                                                               <  2:381494/150‑1 (MQ=255)
                                                                                     tAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAg                                                              <  2:847659/150‑1 (MQ=255)
                                                                                     tAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAg                                                              <  1:2109919/150‑1 (MQ=255)
                                                                                        tctcGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTa                                                                                                                                 <  2:352825/80‑1 (MQ=255)
                                                                                        tctcGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTa                                                                                                                                 >  1:352825/1‑80 (MQ=255)
                                                                                        tctcGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCAt                                                                                 >  1:1843709/1‑128 (MQ=255)
                                                                                        tctcGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCAt                                                                                 <  2:1843709/128‑1 (MQ=255)
                                                                                        tctcGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATgtgg                                                                             <  1:1015653/132‑1 (MQ=255)
                                                                                        tctcGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATgtgg                                                                             >  2:1015653/1‑132 (MQ=255)
                                                                                        tctcGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAgggg                                                           >  2:1057051/1‑150 (MQ=255)
                                                                                        tctcGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAgggg                                                           >  2:1883768/1‑150 (MQ=255)
                                                                                         ctcGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGt                                                          <  2:464955/150‑1 (MQ=255)
                                                                                          tcGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGtt                                                         >  1:2183455/1‑150 (MQ=255)
                                                                                            gACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCAtgt                                                                               >  1:909133/1‑126 (MQ=255)
                                                                                            gACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCAtgt                                                                               <  2:909133/126‑1 (MQ=255)
                                                                                                         cAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtg                                                                                                                              >  1:953969/1‑66 (MQ=255)
                                                                                                         cAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtg                                                                                                                              <  2:953969/66‑1 (MQ=255)
                                                                                                         cAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTgg                                          >  2:2368721/1‑150 (MQ=255)
                                                                                                               gAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTc                                    >  1:322655/1‑150 (MQ=255)
                                                                                                                     tGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCccgccg                              <  2:1597556/150‑1 (MQ=255)
                                                                                                                          gAAAGGCTGATGGCTGGTGGCTTATGACACCGGCCTTTTTCCGTTATTGTATTCCTGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGgaga                         >  1:1750432/1‑150 (MQ=255)
                                                                                                                           aaaGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa                        <  1:751663/150‑1 (MQ=255)
                                                                                                                           aaaGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa                        >  1:2174585/1‑150 (MQ=255)
                                                                                                                           aaaGGCTGATGGCTGGTGGCATATGCCACCGGCCTTTTGCCGTAATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCCGCCGCCGGGAGaa                        >  2:44642/1‑150 (MQ=255)
                                                                                                                             aGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTTCTCCCGGCGGGagcagc                      <  2:1814503/150‑1 (MQ=255)
                                                                                                                              ggCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAg                                                              >  2:2322878/1‑109 (MQ=255)
                                                                                                                              ggCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAg                                                              <  1:2322878/109‑1 (MQ=255)
                                                                                                                              ggCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTTCTCCCGGCGGGAGCAGcc                     <  2:2181536/150‑1 (MQ=255)
                                                                                                                                 tGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCaaa                  >  2:1314026/1‑150 (MQ=255)
                                                                                                                                  gATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAAc                 >  1:101923/1‑150 (MQ=255)
                                                                                                                                  gATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAAc                 <  1:1314026/150‑1 (MQ=255)
                                                                                                                                   aTGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATaa                                                                                                    >  1:912198/1‑66 (MQ=255)
                                                                                                                                   aTGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATaa                                                                                                    <  2:912198/66‑1 (MQ=255)
                                                                                                                                    tGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTTCTCCCGGCGGGAGCAGCCagaaaa               <  2:2183455/150‑6 (MQ=255)
                                                                                                                                    aggccggtggCATAAGCCACCAGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGc                      >  2:1858770/6‑143 (MQ=255)
                                                                                                                                    aggccggtggCATAAGCCACCAGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGc                      <  1:1858770/138‑1 (MQ=255)
                                                                                                                                     ggCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAg              >  1:309553/1‑150 (MQ=255)
                                                                                                                                      gCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGc             <  2:831419/150‑1 (MQ=255)
                                                                                                                                       cTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCa            <  1:10209/150‑1 (MQ=255)
                                                                                                                                        tggtggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGTGAGAAGCCAAACAAGCAg           >  1:2272096/1‑150 (MQ=255)
                                                                                                                                        tggtggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGc             <  2:652101/148‑1 (MQ=255)
                                                                                                                                        tggtggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGc             >  1:652101/1‑148 (MQ=255)
                                                                                                                                        cggtggCATAAGCCACCAGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAg           <  1:2346485/149‑1 (MQ=255)
                                                                                                                                            ggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGa       >  2:973230/1‑150 (MQ=255)
                                                                                                                                               ttATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCa    >  2:58768/1‑150 (MQ=255)
                                                                                                                                                 aTGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAg  <  1:1793115/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAG  >  NC_002947/537318‑537612

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