Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,791,509 C→A 6.4% P231P (CCG→CCT pheT ← phenylalanine tRNA synthetase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,791,5090CA 6.4% 522.0 / 10.5 140P231P (CCG→CCTpheTphenylalanine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (61/70);  new base A (5/4);  total (66/74)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 141 total aligned reads displayed.

CCGCACCACAATGCCGCCTTCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTC  >  CP009273/1791365‑1791656
                                                                                                                                                |                                                                                                                                                   
ccGCACCACAATGCCGCCTTCAATTCGATCTTTATCGAAAGCGTGCATCGGCTGGCCAAGTTCGAGCAGCACATAGTTGGAGACGTCAACAACTGCATAGAGAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGa                                                                                                                                                  >  1:1668515/1‑148 (MQ=255)
ccGCACCACAATGCCGCCTTCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGa                                                                                                                                                  <  2:1510368/148‑1 (MQ=255)
ccGCACCACAATGCCGCCTTCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGa                                                                                                                                                  <  2:1510390/148‑1 (MQ=255)
   caccacAATGCCGCCTTCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGtt                                                                                                                                               <  1:2626718/148‑1 (MQ=255)
       acaATGCCGCCTTCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGgcg                                                                                                                                           <  1:2385299/148‑1 (MQ=255)
        cactGCCGCCTTCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGgcgc                                                                                                                                          <  2:2509559/145‑1 (MQ=255)
        caATGCCGCCTTCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGgcgc                                                                                                                                          >  2:1137741/1‑148 (MQ=255)
         aaTGCCGCCTTCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCt                                                                                                                                         <  2:114674/148‑1 (MQ=255)
         aaTGCCGCCTTCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCt                                                                                                                                         <  1:574088/148‑1 (MQ=255)
         aaTGCCGCCTTCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCt                                                                                                                                         <  2:2697429/148‑1 (MQ=255)
          aTGCCGCCTTCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCtt                                                                                                                                        <  2:2073297/148‑1 (MQ=255)
             ccgccTTCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTaa                                                                                                                                     <  1:3115685/148‑1 (MQ=255)
               gccTTCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACg                                                                                                                                   <  1:1403565/148‑1 (MQ=255)
                ccTTCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGt                                                                                                                                  <  2:2410610/148‑1 (MQ=255)
                ccTTCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGt                                                                                                                                  >  1:3191247/1‑148 (MQ=255)
                  ttCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTa                                                                                                                                <  1:694665/148‑1 (MQ=255)
                  ttCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTa                                                                                                                                <  2:1189691/148‑1 (MQ=255)
                   tCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTaa                                                                                                                               >  2:175084/1‑148 (MQ=255)
                   tCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTaa                                                                                                                               <  2:2771844/148‑1 (MQ=255)
                    cAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAAt                                                                                                                              >  2:2148778/1‑148 (MQ=255)
                      aTGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGc                                                                                                                            <  2:949469/148‑1 (MQ=255)
                       tGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGcc                                                                                                                           >  1:1964251/1‑148 (MQ=255)
                        gCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCt                                                                                                                          <  2:1784213/148‑1 (MQ=255)
                        gCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCt                                                                                                                          >  1:44522/1‑148 (MQ=255)
                        gCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCt                                                                                                                          >  2:1060317/1‑148 (MQ=255)
                         cGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCtt                                                                                                                         <  2:2363836/148‑1 (MQ=255)
                         cGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCtt                                                                                                                         <  1:1107691/148‑1 (MQ=255)
                         cGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCtt                                                                                                                         <  1:1302125/148‑1 (MQ=255)
                         agATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCtt                                                                                                                         <  1:1303578/147‑1 (MQ=255)
                           aTCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCtttt                                                                                                                       <  1:207693/148‑1 (MQ=255)
                           aTCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCtttt                                                                                                                       <  2:1264498/148‑1 (MQ=255)
                            tCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTa                                                                                                                      >  2:728638/1‑148 (MQ=255)
                            tCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTa                                                                                                                      <  2:2561491/148‑1 (MQ=255)
                            tCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTa                                                                                                                      >  2:809171/1‑148 (MQ=255)
                             cTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTAc                                                                                                                     >  1:2208463/1‑148 (MQ=255)
                               ttATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACca                                                                                                                   >  1:1186271/1‑148 (MQ=255)
                                     aaaGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGc                                                                                                             <  1:1764136/148‑1 (MQ=255)
                                      aaGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGcc                                                                                                            >  1:1769311/1‑148 (MQ=255)
                                       aGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCa                                                                                                           >  1:1782119/1‑148 (MQ=255)
                                       aGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCa                                                                                                           >  1:2295689/1‑148 (MQ=255)
                                         cGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGAGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAg                                                                                                         >  2:807131/1‑148 (MQ=255)
                                          gTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGa                                                                                                        <  1:269978/148‑1 (MQ=255)
                                           tGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGAt                                                                                                       >  2:633111/1‑148 (MQ=255)
                                           tGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGAt                                                                                                       >  2:1231440/1‑148 (MQ=255)
                                            gCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATa                                                                                                      <  1:424652/148‑1 (MQ=255)
                                            gCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATa                                                                                                      <  1:685297/148‑1 (MQ=255)
                                            gCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATa                                                                                                      >  2:950614/1‑148 (MQ=255)
                                             cATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATaa                                                                                                     <  1:1253392/148‑1 (MQ=255)
                                                cGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGAGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAAcgc                                                                                                  >  2:3110090/1‑148 (MQ=255)
                                                 ggCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAAcgcg                                                                                                 <  2:817603/148‑1 (MQ=255)
                                                  gCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCgg                                                                                                <  2:2886116/148‑1 (MQ=255)
                                                  gCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCgg                                                                                                <  2:2886174/148‑1 (MQ=255)
                                                  gCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGAGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCgg                                                                                                >  2:1666383/1‑148 (MQ=255)
                                                         cAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCt                                                                                         <  1:2487380/148‑1 (MQ=255)
                                                          aGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCtt                                                                                        <  1:2840103/148‑1 (MQ=255)
                                                            ttCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTcc                                                                                      >  1:1080209/1‑148 (MQ=255)
                                                             tCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCg                                                                                     >  2:3211449/1‑148 (MQ=255)
                                                             tCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCAGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCg                                                                                     >  2:3211909/1‑148 (MQ=255)
                                                                 gcagcaCATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCNACAGAGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGgcgc                                                                                 <  1:1680951/148‑1 (MQ=255)
                                                                 gcagcaCATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGgcgc                                                                                 <  2:2208463/148‑1 (MQ=255)
                                                                   agcaCATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCtt                                                                               <  1:2288518/148‑1 (MQ=255)
                                                                    gcaCATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTc                                                                              >  1:591136/1‑148 (MQ=255)
                                                                     cacaTAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCg                                                                             >  1:82526/1‑148 (MQ=255)
                                                                      acaTAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGa                                                                            <  1:2432148/148‑1 (MQ=255)
                                                                      acaTAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGa                                                                            <  1:2919644/148‑1 (MQ=255)
                                                                           gTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTa                                                                       >  1:189566/1‑148 (MQ=255)
                                                                           gTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCGTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTa                                                                       >  1:2458624/1‑148 (MQ=255)
                                                                            ttGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTaa                                                                      >  1:2091665/1‑148 (MQ=255)
                                                                            ttGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTaa                                                                      >  1:594799/1‑148 (MQ=255)
                                                                                tGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCgg                                                                  <  2:1494284/148‑1 (MQ=255)
                                                                                 gACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGc                                                                 <  1:3252383/148‑1 (MQ=255)
                                                                                        acaacTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCgtgt                                                          >  2:777771/1‑148 (MQ=255)
                                                                                             tGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGtcgtcg                                                     >  2:2530683/1‑148 (MQ=255)
                                                                                             tGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGtcgtcg                                                     <  1:728638/148‑1 (MQ=255)
                                                                                               cATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGAt                                                   >  1:613327/1‑148 (MQ=255)
                                                                                               cATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGAt                                                   >  1:1696958/1‑148 (MQ=255)
                                                                                               cATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGAt                                                   >  2:699059/1‑148 (MQ=255)
                                                                                                     tAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGc                                             >  2:2139970/1‑148 (MQ=255)
                                                                                                      aGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCa                                            >  1:2146754/1‑148 (MQ=255)
                                                                                                      aGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGAGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCa                                            <  2:3090105/148‑1 (MQ=255)
                                                                                                        aaCGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCAcc                                          >  1:1363459/1‑148 (MQ=255)
                                                                                                         aCGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCa                                         <  1:1437476/148‑1 (MQ=255)
                                                                                                                ccGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGaa                                  <  2:1769311/148‑1 (MQ=255)
                                                                                                                      cgacgCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGAttt                            >  1:3261247/1‑148 (MQ=255)
                                                                                                                        acgCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTcc                          <  1:1044074/148‑1 (MQ=255)
                                                                                                                            aGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGtt                      >  2:3179268/1‑148 (MQ=255)
                                                                                                                              tttttCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGa                    <  1:962416/148‑1 (MQ=255)
                                                                                                                              tttttCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGa                    <  1:1705926/148‑1 (MQ=255)
                                                                                                                               ttttCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGaa                   <  1:1425016/148‑1 (MQ=255)
                                                                                                                                tttctttcATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAAc                  <  1:1262974/148‑1 (MQ=255)
                                                                                                                                     ttcATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCg             >  2:1237494/1‑148 (MQ=255)
                                                                                                                                     ttcATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCg             >  2:2223686/1‑148 (MQ=255)
                                                                                                                                     ttcATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCg             >  2:1440184/1‑148 (MQ=255)
                                                                                                                                     ttcATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCg             >  2:1606471/1‑148 (MQ=255)
                                                                                                                                      tcATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCgg            <  2:2687542/148‑1 (MQ=255)
                                                                                                                                         tCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAg         >  2:1358788/1‑148 (MQ=255)
                                                                                                                                         tCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAg         <  1:1060317/148‑1 (MQ=255)
                                                                                                                                             caGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTgg     <  1:2417169/148‑1 (MQ=255)
                                                                                                                                               gCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGtt   >  2:2446096/1‑148 (MQ=255)
                                                                                                                                                cGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTc  >  2:2446593/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CCGCACCACAATGCCGCCTTCAATGCGATCTTTATCGAAAGCGTGCATCGGCTGGCCCAGTTCGAGCAGCACATAGTTGGTGACGTCAACAACTGCATCGATAGAACGGATCCCGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTC  >  CP009273/1791365‑1791656

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: