Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,716,114 C→A 5.3% A548A (GCG→GCT glyS ← glycine tRNA synthetase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,716,1140CA 5.3% 1110.1 / 21.7 298A548A (GCG→GCTglySglycine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (153/129);  new base A (6/10);  total (159/139)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 306 total aligned reads displayed.

CTCAGCACTTCGTCAGATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCACCGCTTCTTCGGTCAGCGTTTGCAGATC  >  CP009273/3715965‑3716263
                                                                                                                                                     |                                                                                                                                                     
  cAGCACTTCGTCAGATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGc                                                                                                                                                     >  2:4991072/1‑150 (MQ=255)
  cAGCACTTCGTCAGATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGc                                                                                                                                                     >  2:3327444/1‑150 (MQ=255)
   aGCACTTCGTCAGATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGcc                                                                                                                                                    >  1:2078868/1‑150 (MQ=255)
   aGCACTTCGTCAGATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACAGcc                                                                                                                                                    >  1:4553908/1‑150 (MQ=255)
    gCACTTCGTCAGATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCt                                                                                                                                                   >  2:4538807/1‑150 (MQ=255)
    gCACTTCGTCAGATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCt                                                                                                                                                   <  2:2265952/150‑1 (MQ=255)
    gCACTTCGTCAGATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCt                                                                                                                                                   <  1:1539655/150‑1 (MQ=255)
    gCACTTCGTCAGATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCt                                                                                                                                                   <  1:5068489/150‑1 (MQ=255)
    gCACTTCGTCAGATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACAGCCt                                                                                                                                                   <  2:1641758/150‑1 (MQ=255)
     cACTTCGTCAGATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTg                                                                                                                                                  >  2:3450389/1‑150 (MQ=255)
          cGTCAGATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATgg                                                                                                                                             >  2:3308133/1‑150 (MQ=255)
           gTCAGATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACAGCCTGGATGgt                                                                                                                                            <  1:1617690/150‑1 (MQ=255)
            tCAGATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGgtg                                                                                                                                           >  1:2192716/1‑150 (MQ=255)
             cAGATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGATGCGCCAGTACCGCCTGGATGgtgt                                                                                                                                          >  2:3546209/1‑150 (MQ=255)
               gATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCa                                                                                                                                        >  1:4204837/1‑150 (MQ=255)
               gATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCa                                                                                                                                        >  2:2040293/1‑150 (MQ=255)
               gATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCa                                                                                                                                        <  2:4375415/150‑1 (MQ=255)
                   tCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACgg                                                                                                                                    <  1:1025672/150‑1 (MQ=255)
                    cGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGt                                                                                                                                   >  2:4863792/1‑150 (MQ=255)
                      ccAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGtat                                                                                                                                 >  1:3187656/1‑150 (MQ=255)
                       cAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACAGCCTGGATGGTGTCAACGGtata                                                                                                                                >  2:3678526/1‑150 (MQ=255)
                        aGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATaa                                                                                                                               >  1:4585315/1‑150 (MQ=255)
                          aaTGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAAcc                                                                                                                             >  2:4552965/1‑150 (MQ=255)
                                  aTACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCt                                                                                                                     >  2:693712/1‑150 (MQ=255)
                                  aTACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACAGCCTGGATGGTGTCAACGGTATAACCTTCGTCCt                                                                                                                     <  2:3851489/150‑1 (MQ=255)
                                   tACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTg                                                                                                                    <  1:4679158/150‑1 (MQ=255)
                                   tACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTg                                                                                                                    <  1:3661316/150‑1 (MQ=255)
                                       cGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATAc                                                                                                                <  2:359480/150‑1 (MQ=255)
                                       cGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATAc                                                                                                                >  2:1148435/1‑150 (MQ=255)
                                       cGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATAc                                                                                                                >  1:3788474/1‑150 (MQ=255)
                                        gTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATAcc                                                                                                               >  2:2553806/1‑150 (MQ=255)
                                         tttGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCa                                                                                                              >  1:1322770/1‑150 (MQ=255)
                                              tCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGgcgc                                                                                                         >  1:35524/1‑150 (MQ=255)
                                                gccgccgccAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCgg                                                                                                       <  1:930629/150‑1 (MQ=255)
                                                  cgccgccAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGaa                                                                                                     <  2:59992/150‑1 (MQ=255)
                                                  cgccgccAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACAGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGaa                                                                                                     >  2:1085164/1‑150 (MQ=255)
                                                   gccgccAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAg                                                                                                    >  1:2765168/1‑150 (MQ=255)
                                                   gccgccAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAg                                                                                                    <  1:1289559/150‑1 (MQ=255)
                                                     cgccAGTGCCGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCg                                                                                                  >  1:1013988/1‑150 (MQ=255)
                                                     cgccAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCg                                                                                                  >  1:1386028/1‑150 (MQ=255)
                                                      gccAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGa                                                                                                 >  2:4741645/1‑150 (MQ=255)
                                                      gccAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGa                                                                                                 <  2:2363891/150‑1 (MQ=255)
                                                        cAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGAcc                                                                                               <  2:1283666/150‑1 (MQ=255)
                                                        cAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGAcc                                                                                               >  1:2959475/1‑150 (MQ=255)
                                                         aGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCg                                                                                              >  1:1881595/1‑150 (MQ=255)
                                                          gTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGa                                                                                             >  1:497124/1‑150 (MQ=255)
                                                             cagcagCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCa                                                                                          >  2:1522778/1‑150 (MQ=255)
                                                               gcagcTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACAGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATa                                                                                        >  2:4593506/1‑150 (MQ=255)
                                                                cagcTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATaa                                                                                       >  1:289077/1‑150 (MQ=255)
                                                                 agcTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATaaa                                                                                      >  1:2145104/1‑150 (MQ=255)
                                                                          cAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAAc                                                                             <  1:4143678/150‑1 (MQ=255)
                                                                            gggTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAAcat                                                                           >  1:2530288/1‑150 (MQ=255)
                                                                             ggTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAAcatc                                                                          <  2:1130080/150‑1 (MQ=255)
                                                                              gTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAAcatca                                                                         <  1:4891684/150‑1 (MQ=255)
                                                                                aCGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACatcatc                                                                       >  2:4778555/1‑150 (MQ=255)
                                                                                 cGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACAtcatca                                                                      >  1:2039575/1‑150 (MQ=255)
                                                                                     aaTGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTa                                                                  <  1:1428687/150‑1 (MQ=255)
                                                                                       tGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACg                                                                <  2:3255359/150‑1 (MQ=255)
                                                                                       tGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACg                                                                <  1:3013084/150‑1 (MQ=255)
                                                                                        gCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGt                                                               <  1:443281/150‑1 (MQ=255)
                                                                                           aTACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTgg                                                            <  1:4391266/150‑1 (MQ=255)
                                                                                               cGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCAtt                                                        >  2:1731456/1‑150 (MQ=255)
                                                                                                gCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTa                                                       <  1:3389990/150‑1 (MQ=255)
                                                                                                 cTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAg                                                      >  2:703585/1‑150 (MQ=255)
                                                                                                 cTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAg                                                      >  1:3670661/1‑150 (MQ=255)
                                                                                                   ttCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTc                                                    <  1:5115631/150‑1 (MQ=255)
                                                                                                    tCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCa                                                   <  1:1585200/150‑1 (MQ=255)
                                                                                                     cATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAg                                                  >  1:883929/1‑150 (MQ=255)
                                                                                                     cATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAg                                                  <  1:1809674/150‑1 (MQ=255)
                                                                                                       tACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCt                                                <  1:3450389/150‑1 (MQ=255)
                                                                                                         cGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTa                                              <  2:4514718/150‑1 (MQ=255)
                                                                                                          gggCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTAt                                             >  2:3702187/1‑150 (MQ=255)
                                                                                                             cATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGc                                          <  1:4362885/150‑1 (MQ=255)
                                                                                                             cATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCAGTAGTCAGCTTATCGc                                          >  2:1372364/1‑150 (MQ=255)
                                                                                                                 gAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAACATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATa                                      >  1:477726/1‑150 (MQ=255)
                                                                                                                  aaaTCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATAc                                     <  1:201189/150‑1 (MQ=255)
                                                                                                                  aaaTCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATAc                                     <  2:2907263/150‑1 (MQ=255)
                                                                                                                  aaaTCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATAc                                     >  1:4043035/1‑150 (MQ=255)
                                                                                                                    aTCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAg                                   >  1:3572537/1‑150 (MQ=255)
                                                                                                                       aGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAGGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACg                                >  1:542715/1‑150 (MQ=255)
                                                                                                                         ccGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCa                              <  1:4824085/150‑1 (MQ=255)
                                                                                                                         ccGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCa                              <  2:3140929/150‑1 (MQ=255)
                                                                                                                          cGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCAc                             >  2:3883048/1‑150 (MQ=255)
                                                                                                                           ggACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCAcc                            <  1:795890/150‑1 (MQ=255)
                                                                                                                            gACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCACCg                           <  2:4204837/150‑1 (MQ=255)
                                                                                                                             aCGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCACCGc                          >  1:4253101/1‑150 (MQ=255)
                                                                                                                              cGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCACCGct                         <  1:961226/150‑1 (MQ=255)
                                                                                                                               gAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCACCGctt                        >  2:4647486/1‑150 (MQ=255)
                                                                                                                                   cGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCACCGCttcttc                    >  2:405122/1‑150 (MQ=255)
                                                                                                                                   cGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCACCGCttcttc                    <  1:3346700/150‑1 (MQ=255)
                                                                                                                                   cGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCACCGCttcttc                    <  2:497124/150‑1 (MQ=255)
                                                                                                                                    ggACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCACCGCTTCTTCg                   >  2:3039515/1‑150 (MQ=255)
                                                                                                                                      acgacgCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCACCGCTTCTTCGGt                 <  1:65865/150‑1 (MQ=255)
                                                                                                                                      acgacgCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCACCGCTTCTTCGGt                 <  2:5149582/150‑1 (MQ=255)
                                                                                                                                       cgacgCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCACCGCTTCTTCGGTc                >  1:3866833/1‑150 (MQ=255)
                                                                                                                                           gcgcCAGTACCGCCTGGAGGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCACCGCTTCTTCGGTCAGCg            <  2:4915831/150‑1 (MQ=255)
                                                                                                                                               cAGTACAGCCTGGTTGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCACCGCTTCTTCGGTCAGCGTTTg        <  1:3687350/150‑1 (MQ=255)
                                                                                                                                                aGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCACCGCTTCTTCGGTCAGCGTTTGc       >  1:4181367/1‑150 (MQ=255)
                                                                                                                                                 gTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCACCGCTTCTTCGGTCAGCGTTTGCa      >  2:2735125/1‑150 (MQ=255)
                                                                                                                                                     cGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCACCGCTTCTTCGGTCAGCGTTTGCAGATc  <  1:1712270/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTCAGCACTTCGTCAGATTTCGCCAGAATGTTAGATACACGTTTGTTCGCCGCCGCCAGTGCAGCAGCTGCATCCAGGGTACGGAAATGCGATACCGCTTTCATACGGGCATCGAAATCAGCCGGACGAGTCGGACGACGCGCCAGTACCGCCTGGATGGTGTCAACGGTATAACCTTCGTCCTGATACCAGGCGCGGAAGCGACCGAGCATAAAGTCGATAACATCATCAACTACGTTGGCATTAGTCAGCTTATCGCCATACAGACGCACCGCTTCTTCGGTCAGCGTTTGCAGATC  >  CP009273/3715965‑3716263

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: