Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,127,712 T→C 6.2% intergenic (+29/+118) PP_0983 → / ← PP_0984 membrane protein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,127,7120TC 6.2% 878.6 / 28.7 275intergenic (+29/+118)PP_0983/PP_0984membrane protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (135/123);  new base C (7/10);  total (142/133)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 276 total aligned reads displayed.

GGCCAGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCACCCGCACCAGCT  >  NC_002947/1127564‑1127861
                                                                                                                                                    |                                                                                                                                                     
    aGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAggg                                                                                                                                                  <  1:6329821/150‑1 (MQ=255)
      gACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTc                                                                                                                                                <  2:6489728/150‑1 (MQ=255)
         cTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGc                                                                                                                                             <  1:3215756/150‑1 (MQ=255)
             ttGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGTGCTGCACAAGAAAGCCGACCTCAGGGTCGGCtgtt                                                                                                                                         >  2:561494/1‑150 (MQ=255)
             ttGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCtttc                                                                                                                                         <  2:594765/150‑2 (MQ=255)
                 gCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGt                                                                                                                                     >  1:6305366/1‑150 (MQ=255)
                 gCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGt                                                                                                                                     <  1:4649485/150‑1 (MQ=255)
                 gCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCCTGAGGTCGGCtttcttgt                                                                                                                                     <  2:4248907/150‑6 (MQ=255)
                  ccGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACGAAAAAGCCGACCCTGAGGTCGGCTTTCTtgtg                                                                                                                                    <  1:6656074/150‑1 (MQ=38)
                  ccGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTg                                                                                                                                    <  1:738778/150‑1 (MQ=255)
                  ccGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTg                                                                                                                                    <  1:6618875/150‑1 (MQ=255)
                  ccGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTg                                                                                                                                    <  2:1688494/150‑1 (MQ=255)
                   cGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTGTTCGTgg                                                                                                                                   >  2:109846/1‑150 (MQ=255)
                   cGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCCTGAGGTCGGCTTTCTtgtgc                                                                                                                                   <  1:1018338/150‑2 (MQ=255)
                    gTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTggg                                                                                                                                  >  2:1448422/1‑150 (MQ=255)
                         gCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTg                                                                                                                             >  2:3998609/1‑150 (MQ=255)
                             ggTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCt                                                                                                                         <  2:6181937/150‑1 (MQ=255)
                              gTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTg                                                                                                                        >  2:2077600/1‑150 (MQ=255)
                               tATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTgg                                                                                                                       >  1:4963101/1‑150 (MQ=255)
                               tATTCGGCGTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTgg                                                                                                                       <  2:6606300/150‑1 (MQ=255)
                                 ttCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGcc                                                                                                                     >  1:640801/1‑150 (MQ=255)
                                   cGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCg                                                                                                                   <  1:5424943/150‑1 (MQ=255)
                                   cGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGCGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCg                                                                                                                   <  1:1192931/150‑1 (MQ=255)
                                    ggCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCgg                                                                                                                  >  2:5350099/1‑150 (MQ=255)
                                     gCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGc                                                                                                                 <  2:221625/150‑1 (MQ=255)
                                       ttCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCct                                                                                                               >  2:5181931/1‑150 (MQ=255)
                                       ttCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCct                                                                                                               >  2:1182181/1‑150 (MQ=255)
                                          ccGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCtt                                                                                                            >  1:5327587/1‑150 (MQ=255)
                                          ccGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTGTTCGTGGGTGTTGATCTGGCCCGGCCTCtt                                                                                                            >  1:2233503/1‑150 (MQ=255)
                                            gCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTcg                                                                                                          >  1:5251043/1‑150 (MQ=255)
                                            gCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTcg                                                                                                          >  1:1249537/1‑150 (MQ=255)
                                            gCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTcg                                                                                                          >  1:469311/1‑150 (MQ=255)
                                            gCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTcg                                                                                                          >  1:5654386/1‑150 (MQ=255)
                                              cACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTcgcg                                                                                                        >  2:2583026/1‑150 (MQ=255)
                                              cACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTcgcg                                                                                                        >  1:34547/1‑150 (MQ=255)
                                                 tgctgGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGt                                                                                                     >  1:2890750/1‑150 (MQ=255)
                                                 tgctgGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGt                                                                                                     >  1:6135087/1‑150 (MQ=255)
                                                 tgctgGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTGTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCgggg                                                                                                     >  2:1032845/1‑149 (MQ=255)
                                                   ctgGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTaa                                                                                                   <  2:6798414/150‑1 (MQ=255)
                                                      gCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGcc                                                                                                >  1:2384302/1‑150 (MQ=255)
                                                         gtggtgATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCcgc                                                                                             <  2:5046396/150‑1 (MQ=255)
                                                         gtggtgATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCcgc                                                                                             >  1:1465620/1‑150 (MQ=255)
                                                            gtgATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGcc                                                                                          >  2:6365796/1‑150 (MQ=255)
                                                            gtgATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAAGAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCcgcgaa                                                                                          <  2:1373278/150‑3 (MQ=255)
                                                              gATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTa                                                                                        >  1:1955508/1‑150 (MQ=255)
                                                              gATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTa                                                                                        >  1:5186803/1‑150 (MQ=255)
                                                                   cAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGt                                                                                   >  1:2729751/1‑150 (MQ=255)
                                                                    aGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTa                                                                                  >  2:2509328/1‑150 (MQ=255)
                                                                          cATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCgcgc                                                                            <  1:4205553/150‑1 (MQ=255)
                                                                           aTCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCa                                                                           >  2:2627379/1‑150 (MQ=255)
                                                                                  ccTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCa                                                                    <  1:6242800/150‑1 (MQ=255)
                                                                                  ccTTGGCCGGTGTGTGGTTGTGGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCa                                                                    <  1:248486/150‑1 (MQ=255)
                                                                                   cTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTaa                                                                                                   >  1:3681456/1‑118 (MQ=255)
                                                                                   cTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAg                                                                   >  1:5841686/1‑150 (MQ=255)
                                                                                    ttGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAgg                                                                  >  2:3567053/1‑150 (MQ=255)
                                                                                    ttGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAgg                                                                  >  1:3308947/1‑150 (MQ=255)
                                                                                    ttGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAgg                                                                  >  1:3006769/1‑150 (MQ=255)
                                                                                     tGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGa                                                                 <  1:6666139/150‑1 (MQ=255)
                                                                                          ggTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAggg                                                            <  1:2120738/150‑1 (MQ=255)
                                                                                              tgtgGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGTCCg                                                        >  1:904796/1‑150 (MQ=255)
                                                                                              tgtgGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCg                                                        >  1:6775149/1‑150 (MQ=255)
                                                                                               gtgGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGa                                                       >  2:2892936/1‑150 (MQ=255)
                                                                                               gtgGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGa                                                       <  1:3903433/150‑1 (MQ=255)
                                                                                                tgGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGCAAGGGCCCGaa                                                      >  2:4438921/1‑150 (MQ=255)
                                                                                                tgGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGaa                                                      >  1:5019689/1‑150 (MQ=255)
                                                                                                 ggTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAg                                                     <  2:4692735/150‑1 (MQ=255)
                                                                                                 ggTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAg                                                     >  1:1112113/1‑150 (MQ=255)
                                                                                                  gttgttGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAgg                                                    >  1:5377303/1‑150 (MQ=255)
                                                                                                     gttgCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCaa                                                 >  1:6234078/1‑150 (MQ=255)
                                                                                                         cgcCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTcc                                             >  2:828592/1‑150 (MQ=255)
                                                                                                          gcCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCa                                            >  1:1751472/1‑150 (MQ=255)
                                                                                                          gcCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCa                                            >  2:6267654/1‑150 (MQ=255)
                                                                                                          gcCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCa                                            >  1:5942801/1‑150 (MQ=255)
                                                                                                          gcCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGAATCAGGAAAGGGCCCGAAGGCAACTCCa                                            <  1:974014/150‑1 (MQ=255)
                                                                                                           ccGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCaa                                           <  2:2460202/150‑1 (MQ=255)
                                                                                                            cGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAAt                                          >  1:4812996/1‑150 (MQ=255)
                                                                                                             gggCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATc                                         <  1:5828/150‑1 (MQ=255)
                                                                                                                  gggTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTg                                    <  1:3203782/150‑1 (MQ=255)
                                                                                                                   ggTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTgg                                   >  2:3316455/1‑150 (MQ=255)
                                                                                                                    gTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGc                                  <  2:5016429/150‑1 (MQ=255)
                                                                                                                    gTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGc                                  <  1:4056718/150‑1 (MQ=255)
                                                                                                                          gCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTAttt                            >  2:5942561/1‑150 (MQ=255)
                                                                                                                             tGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttctt                         >  2:5508863/1‑150 (MQ=255)
                                                                                                                             tGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttctt                         >  1:4014291/1‑150 (MQ=255)
                                                                                                                             tGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttctt                         >  1:815269/1‑150 (MQ=255)
                                                                                                                              gAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttcttc                        <  1:4365679/150‑1 (MQ=255)
                                                                                                                              gAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttcttc                        <  2:2769015/150‑1 (MQ=255)
                                                                                                                                 cTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGc                     <  2:5093827/150‑1 (MQ=255)
                                                                                                                                 cTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTACTGc                     <  2:4573090/150‑1 (MQ=255)
                                                                                                                                     acgaaaaaGCCGACCCTGAGGTCGGCTTTCTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCg                 >  2:1757368/6‑150 (MQ=255)
                                                                                                                                     acaAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCg                 <  1:2027693/150‑1 (MQ=255)
                                                                                                                                      caAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCgg                >  1:3598852/1‑150 (MQ=255)
                                                                                                                                       gaaaaaGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTCCCGCGCAGGCATCAGGACCGGGCCCGACGGCACCTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGc               >  2:3703940/4‑150 (MQ=255)
                                                                                                                                       gaaaaaGCCGACCCTGAGGTCGGCTTTCTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGc               >  1:156977/4‑150 (MQ=255)
                                                                                                                                       aaGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGc               >  2:5766192/1‑150 (MQ=255)
                                                                                                                                       aaGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGc               >  1:6828103/1‑150 (MQ=255)
                                                                                                                                        aGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGGGTTGATCTGGCCCGGCCTCGTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCCTCCGGCAAGGGCCCGAAGGCAACTCCACTCCTCTGGCTTATTTCTTCTGCGTCGGCa              >  2:1533466/1‑150 (MQ=255)
                                                                                                                                          aaaGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCAcc            >  1:4439099/1‑150 (MQ=255)
                                                                                                                                            aGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCACCCg          >  1:5496215/1‑150 (MQ=255)
                                                                                                                                                   cTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCACCCGCACCAGc   >  1:3912442/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCCAGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCACCCGCACCAGCT  >  NC_002947/1127564‑1127861

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