Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,388,575 G→A 5.4% P471S (CCG→TCG)  proS ← proline‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,388,5750GA 5.4% 825.6 / 22.7 258P471S (CCG→TCG) proSproline‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (133/111);  new base A (5/9);  total (138/120)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 261 total aligned reads displayed.

GGCAAACTTGATGCCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCC  >  NC_002947/1388426‑1388723
                                                                                                                                                       |                                                                                                                                                    
 gCAAACTTGATGCCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGa                                                                                                                                                     >  2:1059684/1‑150 (MQ=255)
 gCAAACTTGATGCCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGa                                                                                                                                                     >  1:5817625/1‑150 (MQ=255)
 gCAAACTTGATGCCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGa                                                                                                                                                     <  1:5020314/150‑1 (MQ=255)
 gCAAACTTGATGCCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGa                                                                                                                                                     <  1:2698957/150‑1 (MQ=255)
 gCAAACTTGATGCCTGGGCAGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGAA‑‑ATCTGGAACGGa                                                                                                                                                     >  2:6044683/1‑150 (MQ=255)
   aaaCTTGATGCCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGc                                                                                                                                                   <  1:2502601/150‑1 (MQ=255)
     aCTTGATGCCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCa                                                                                                                                                 <  1:3152621/150‑1 (MQ=255)
       ttGATGCCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAgg                                                                                                                                               <  2:4308139/150‑1 (MQ=255)
       ttGATGCCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAgg                                                                                                                                               <  2:4899116/150‑1 (MQ=255)
         gATGCCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGc                                                                                                                                             <  1:3839464/150‑1 (MQ=255)
         gATGCCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGc                                                                                                                                             <  1:2449531/150‑1 (MQ=255)
                gggCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTc                                                                                                                                      <  1:1083687/150‑1 (MQ=255)
                gggCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTc                                                                                                                                      <  1:6364847/150‑1 (MQ=255)
                  gCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCa                                                                                                                                    <  2:6866032/150‑1 (MQ=255)
                  gCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCa                                                                                                                                    <  2:4644561/150‑1 (MQ=255)
                  gCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCa                                                                                                                                    <  2:466286/150‑1 (MQ=255)
                  gCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCa                                                                                                                                    <  2:5146335/150‑1 (MQ=255)
                   cTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAg                                                                                                                                   <  1:6433395/150‑1 (MQ=255)
                    tGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACAAGGGCATCATCTGGAACGACGCCCTGGCTCCGTTCCa                                                                                                                                    <  1:6806514/150‑1 (MQ=37)
                      gTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAgatg                                                                                                                                >  1:3284147/1‑150 (MQ=255)
                       ttttCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAgatga                                                                                                                               >  1:3101193/1‑150 (MQ=255)
                       ttttCTTGACGCTGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGTTTTCGTAGCGCAGTGGGAACAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAgatga                                                                                                                               >  2:3180668/1‑150 (MQ=255)
                         ttCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGatgatg                                                                                                                             >  1:4557670/1‑150 (MQ=255)
                           cTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCATCATCTGGAACGACGCCCTGGCTCCGTTCCAgattgcc                                                                                                                             <  2:5588990/150‑5 (MQ=37)
                            ttGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGccc                                                                                                                          >  2:3983701/1‑150 (MQ=255)
                            ttGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGccc                                                                                                                          <  2:6183565/150‑1 (MQ=255)
                            ttGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGccc                                                                                                                          <  1:2263748/150‑1 (MQ=255)
                            ttGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGccc                                                                                                                          >  1:3972822/1‑150 (MQ=255)
                            ttGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGACAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGccc                                                                                                                          <  1:3442236/150‑1 (MQ=255)
                              gTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTGTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCtt                                                                                                                        >  2:5961943/1‑150 (MQ=255)
                                 gcgGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTgtc                                                                                                                     >  1:5447468/1‑150 (MQ=255)
                                  cgGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTgtcg                                                                                                                    >  1:5172807/1‑150 (MQ=255)
                                     tcgtcCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCa                                                                                                                 >  2:640977/1‑150 (MQ=255)
                                       gtcCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATa                                                                                                               >  1:5243179/1‑150 (MQ=255)
                                        tcCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAg                                                                                                              <  1:3822563/150‑1 (MQ=255)
                                        tcCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAg                                                                                                              <  2:1449709/150‑1 (MQ=255)
                                          cagcagCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGct                                                                                                            <  1:4317565/150‑1 (MQ=255)
                                           agcagcACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGctc                                                                                                           >  1:4464768/1‑150 (MQ=255)
                                              agcaCTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGc                                                                                                        >  2:592575/1‑150 (MQ=255)
                                                caCTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTc                                                                                                      >  1:2686314/1‑150 (MQ=255)
                                                caCTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTc                                                                                                      <  2:3900514/150‑1 (MQ=255)
                                                   ttCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGag                                                                                                   >  2:4166661/1‑149 (MQ=255)
                                                     cATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATgg                                                                                                 >  1:4669864/1‑150 (MQ=255)
                                                       tAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCa                                                                                               >  1:5594330/1‑150 (MQ=255)
                                                        aGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAg                                                                                              >  1:5096130/1‑150 (MQ=255)
                                                         gCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGc                                                                                             <  1:904457/150‑1 (MQ=255)
                                                          ccGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCt                                                                                            >  2:1260936/1‑150 (MQ=255)
                                                             gcggcAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCa                                                                                         <  1:3552219/150‑1 (MQ=255)
                                                                 cagtcagtTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACca                                                                                     <  1:208398/150‑1 (MQ=255)
                                                                  agtcagtTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACcac                                                                                    >  1:1496755/1‑150 (MQ=255)
                                                                    tcagtTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCacac                                                                                  <  2:744065/150‑1 (MQ=255)
                                                                    tcagtTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCacac                                                                                  >  2:871020/1‑150 (MQ=255)
                                                                      agtTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACgg                                                                                >  2:718547/1‑150 (MQ=255)
                                                                       gtTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACggg                                                                               >  1:2606292/1‑150 (MQ=255)
                                                                        ttCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGa                                                                              <  1:2056785/150‑1 (MQ=255)
                                                                        ttCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGa                                                                              >  2:3420886/1‑150 (MQ=255)
                                                                             cATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGc                                                                         <  2:4495000/150‑1 (MQ=255)
                                                                                  aGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATg                                                                    >  2:6310576/1‑150 (MQ=255)
                                                                                  aGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATg                                                                    >  1:4988361/1‑150 (MQ=255)
                                                                                       gTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACCCGGGACACGCCAATGCCGTa                                                               <  2:3654623/150‑1 (MQ=255)
                                                                                       gTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                               <  2:2172501/150‑1 (MQ=255)
                                                                                       gTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                               <  2:5417433/150‑1 (MQ=255)
                                                                                        tCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAg                                                              <  1:5120916/150‑1 (MQ=255)
                                                                                           gTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCaa                                                           <  2:1900140/150‑1 (MQ=255)
                                                                                           gTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCaa                                                           <  2:5206080/150‑1 (MQ=255)
                                                                                           gTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCaa                                                           >  2:2862948/1‑150 (MQ=255)
                                                                                            ttGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAAc                                                          >  2:2090968/1‑150 (MQ=255)
                                                                                                ttCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCAt                                                      >  1:4264255/1‑150 (MQ=255)
                                                                                                     ggACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACa                                                 <  1:572125/150‑1 (MQ=255)
                                                                                                        ccACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTa                                              >  2:2617430/1‑150 (MQ=255)
                                                                                                        ccACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTa                                              <  2:1637562/150‑1 (MQ=255)
                                                                                                         cACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTAc                                             <  2:3530097/150‑1 (MQ=255)
                                                                                                         cACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTAc                                             >  1:6623094/1‑150 (MQ=255)
                                                                                                         cACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTAc                                             >  2:1836165/1‑150 (MQ=255)
                                                                                                          aCGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACg                                            >  1:4837213/1‑150 (MQ=255)
                                                                                                          aCGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACg                                            >  2:4832385/1‑150 (MQ=255)
                                                                                                          aCGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACg                                            >  1:882438/1‑150 (MQ=255)
                                                                                                          aCGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACg                                            >  2:1386985/1‑150 (MQ=255)
                                                                                                                 tttCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCt                                     >  1:1577731/1‑150 (MQ=255)
                                                                                                                 tttCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCt                                     >  2:1160177/1‑150 (MQ=255)
                                                                                                                 tttCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCt                                     <  2:5817625/150‑1 (MQ=255)
                                                                                                                 tttCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCt                                     >  2:2347710/1‑150 (MQ=255)
                                                                                                                     gTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGcc                                 <  2:5087343/150‑1 (MQ=255)
                                                                                                                       aGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGt                               >  1:6557352/1‑150 (MQ=255)
                                                                                                                           cAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTtctc                           <  2:3681974/150‑1 (MQ=255)
                                                                                                                              tGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGcc                        >  2:1707310/1‑150 (MQ=255)
                                                                                                                                ggACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCa                      <  1:4872597/150‑1 (MQ=255)
                                                                                                                                 gACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAg                     <  1:6125479/150‑1 (MQ=255)
                                                                                                                                 gACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAg                     <  1:2018436/150‑1 (MQ=255)
                                                                                                                                  aCCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGc                    >  1:3862152/1‑150 (MQ=255)
                                                                                                                                   ccAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCa                   >  2:4809422/1‑150 (MQ=255)
                                                                                                                                   ccAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCa                   >  2:2502542/1‑150 (MQ=255)
                                                                                                                                    cAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCAc                  >  1:1160223/1‑150 (MQ=255)
                                                                                                                                     aGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCAcc                 >  2:5336354/1‑150 (MQ=255)
                                                                                                                                       ggCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTg               <  1:1839241/150‑1 (MQ=255)
                                                                                                                                         cA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGc             >  2:4707760/1‑150 (MQ=255)
                                                                                                                                                 ggAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCa       >  1:5980426/1‑150 (MQ=255)
                                                                                                                                                  gAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACa                                                 <  1:763330/107‑1 (MQ=255)
                                                                                                                                                   aaCGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAgc     <  2:6221262/150‑1 (MQ=255)
                                                                                                                                                    aCGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAgcg    >  1:6023385/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
GGCAAACTTGATGCCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCC  >  NC_002947/1388426‑1388723

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