Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,986,289 T→A 61.8% N4Y (AAT→TAT)  pgsA ← phosphatidylglycerophosphate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,986,2890TA61.8% 228.0 / 296.8 225N4Y (AAT→TAT) pgsAphosphatidylglycerophosphate synthetase
Reads supporting (aligned to +/- strand):  ref base T (46/40);  new base A (48/91);  total (94/131)
Fisher's exact test for biased strand distribution p-value = 5.63e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 233 total aligned reads displayed.

TGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGCGTT  >  CP009273/1986140‑1986438
                                                                                                                                                     |                                                                                                                                                     
 gCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTa                                                                                                                                                      <  2:4763739/150‑1 (MQ=255)
  ccAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGACTGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATaaa                                                                                                                                                     <  2:907352/150‑1 (MQ=255)
  ccAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTaa                                                                                                                                                     >  1:3060376/1‑150 (MQ=255)
   cAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTaaa                                                                                                                                                    <  1:387236/150‑1 (MQ=255)
       aaaCCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGc                                                                                                                                                >  2:2030617/1‑150 (MQ=255)
        aaCCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCa                                                                                                                                               >  1:4792424/1‑150 (MQ=255)
        aaCCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCa                                                                                                                                               <  1:4354921/150‑1 (MQ=255)
          ccATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATa                                                                                                                                             <  2:690115/150‑1 (MQ=255)
          ccATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATa                                                                                                                                             >  1:2683983/1‑150 (MQ=255)
          ccATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATa                                                                                                                                             >  1:4784558/1‑150 (MQ=255)
              cGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGa                                                                                                                                         <  1:15395/150‑1 (MQ=255)
               gAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGAc                                                                                                                                        <  2:3563915/150‑1 (MQ=255)
                 aCCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACgg                                                                                                                                      >  2:2392345/1‑150 (MQ=255)
                     gtcagtcaCCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTaa                                                                                                                                  >  1:3720862/1‑150 (MQ=255)
                      tcagtcaCCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAAc                                                                                                                                 <  2:125326/150‑1 (MQ=255)
                       cagtcaCCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACt                                                                                                                                <  1:1995497/150‑1 (MQ=255)
                       cagtcaCCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGAGAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACt                                                                                                                                >  2:2225507/1‑150 (MQ=255)
                        agtcaCCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTa                                                                                                                               >  1:1267040/1‑150 (MQ=255)
                         gtcaCCGCCGCGACGCAGAACATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCCACGTAGGGAGATTAAAGTGCATAATGAAGGGTAACTAt                                                                                                                              >  1:2134000/1‑150 (MQ=255)
                              cgccgcGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCtgtt                                                                                                                         <  2:4551765/150‑1 (MQ=255)
                              cgccgcGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCtgtt                                                                                                                         <  2:1878078/150‑1 (MQ=255)
                              cgccgcGACGCAGAAAATGAGCGCGGCGGCAAACGGCAACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCtgtt                                                                                                                         <  1:660234/150‑1 (MQ=255)
                               gccgcGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCtgttg                                                                                                                        >  2:2321439/1‑150 (MQ=255)
                                  gcgACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCa                                                                                                                     <  2:3813945/150‑1 (MQ=255)
                                     aCGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTa                                                                                                                  <  2:3246129/150‑1 (MQ=255)
                                      cGCAGAAAATGAGCGCGGCGGCAACCGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAACTTGCATAATGACGGGTAACTATCTGTTGTCAGTaa                                                                                                                 <  2:2975176/150‑1 (MQ=255)
                                      cGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTaa                                                                                                                 >  2:4580373/1‑150 (MQ=255)
                                       gCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAg                                                                                                                <  1:83345/150‑1 (MQ=255)
                                       gCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAg                                                                                                                <  1:3255682/150‑1 (MQ=255)
                                       gCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAg                                                                                                                <  1:513519/150‑1 (MQ=255)
                                         aGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGAt                                                                                                              >  2:2100027/1‑150 (MQ=255)
                                          gAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGAtt                                                                                                             >  1:3544936/1‑150 (MQ=255)
                                          gAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGAtt                                                                                                             <  1:2427719/150‑1 (MQ=255)
                                           aaaaTGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTa                                                                                                            >  1:3219542/1‑150 (MQ=255)
                                               tGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTAcccc                                                                                                        <  2:4314304/150‑1 (MQ=255)
                                               tGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTAcccc                                                                                                        <  2:3526040/150‑1 (MQ=255)
                                                 aGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTa                                                                                                      <  1:3971479/150‑1 (MQ=255)
                                                   cgcgGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGACTGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATg                                                                                                    <  2:814095/150‑1 (MQ=255)
                                                   cgcgGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATg                                                                                                    <  2:3244568/150‑1 (MQ=255)
                                                      ggcggcAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTg                                                                                                 <  1:3654195/150‑1 (MQ=255)
                                                      ggcggcAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTg                                                                                                 >  2:4729129/1‑150 (MQ=255)
                                                        cggcAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCt                                                                                               >  1:696364/1‑150 (MQ=255)
                                                        cggcAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCt                                                                                               <  1:1701305/150‑1 (MQ=255)
                                                          gcAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTAc                                                                                             <  1:2049650/150‑1 (MQ=255)
                                                           cAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACa                                                                                            >  1:2383425/1‑150 (MQ=255)
                                                             aaCGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACaga                                                                                          >  1:2001269/1‑150 (MQ=255)
                                                              aCGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACagag                                                                                         >  2:1057798/1‑150 (MQ=255)
                                                              aCGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACagag                                                                                         >  2:1653009/1‑150 (MQ=255)
                                                              aCGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACagag                                                                                         <  2:2939917/150‑1 (MQ=255)
                                                                  cGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTCTGTTGCTCCAGAGACAt                                                                                     >  2:2113967/1‑150 (MQ=255)
                                                                  cGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACAt                                                                                     <  2:1500771/150‑1 (MQ=255)
                                                                     ccAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCaa                                                                                  <  1:71770/150‑1 (MQ=255)
                                                                       aGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATg                                                                                >  2:2816803/1‑150 (MQ=255)
                                                                       aGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATg                                                                                <  1:1324852/150‑1 (MQ=255)
                                                                            aCAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCa                                                                           <  1:4557572/150‑1 (MQ=255)
                                                                             cAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCaa                                                                          >  1:4830172/1‑150 (MQ=255)
                                                                              aGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAAc                                                                         <  1:4677301/150‑1 (MQ=255)
                                                                                      aaaaGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCagaa                                                                 >  1:1174235/1‑150 (MQ=255)
                                                                                       aaaGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCagaag                                                                <  2:3421569/150‑1 (MQ=255)
                                                                                         aGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGAt                                                              >  1:3394639/1‑150 (MQ=255)
                                                                                         aGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGAt                                                              <  2:931147/150‑1 (MQ=255)
                                                                                           aCCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCt                                                            <  2:4423229/150‑1 (MQ=255)
                                                                                              aaTACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCtttt                                                         >  1:3458011/1‑150 (MQ=255)
                                                                                               aTACAAAGAATGGGATAAGGATGACACGGAACTGTTTAAGCAACGTAGGGATATAAAATTGCATACTGACGGGTAACTTTCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGTGACATCAATGTTTCAACGACCAGAAGATCTTTTc                                                        <  1:1950508/150‑1 (MQ=255)
                                                                                               aTACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTc                                                        >  1:955399/1‑150 (MQ=255)
                                                                                                     aGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAg                                                  >  1:3087884/1‑150 (MQ=255)
                                                                                                      gAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGa                                                 <  1:2247455/150‑1 (MQ=255)
                                                                                                        aTGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGAcc                                               >  1:894723/1‑150 (MQ=255)
                                                                                                          gggATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCtt                                             >  1:2317683/1‑150 (MQ=255)
                                                                                                             aTAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCg                                          >  2:3772820/1‑150 (MQ=255)
                                                                                                             aTAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAAGTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCCGAAGATCTTTTCTGCCAGACCTTGCg                                          <  2:4608207/150‑1 (MQ=255)
                                                                                                             aTAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAAGTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCg                                          <  2:2157309/150‑1 (MQ=255)
                                                                                                               aaGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTGGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGaa                                        <  1:218923/150‑1 (MQ=255)
                                                                                                                aGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGaaa                                       >  2:658255/1‑150 (MQ=255)
                                                                                                                 ggATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAAt                                      <  2:1642321/150‑1 (MQ=255)
                                                                                                                 ggATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAAt                                      <  2:3182491/150‑1 (MQ=255)
                                                                                                                       cacGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCgg                                <  1:138031/150‑1 (MQ=255)
                                                                                                                       cacGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCgg                                <  1:146417/150‑1 (MQ=255)
                                                                                                                          ggAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCAc                             >  1:3500742/1‑150 (MQ=255)
                                                                                                                           gAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACt                            <  2:2540410/150‑1 (MQ=255)
                                                                                                                            aaCAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCAAAATACCCGGCACtt                           <  2:1290825/150‑1 (MQ=255)
                                                                                                                              cAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATTTTGCTACAGAGACATAAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACtttt                         <  2:3948501/150‑1 (MQ=255)
                                                                                                                                gtgtAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCAAAATACCCGGCACTTTTGc                       >  1:1669681/1‑150 (MQ=255)
                                                                                                                                gtgtAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGc                       <  1:168288/150‑1 (MQ=255)
                                                                                                                                 tgtAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCa                      >  2:3865838/1‑150 (MQ=255)
                                                                                                                                    aaGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAAtt                   >  1:4187653/1‑150 (MQ=255)
                                                                                                                                       cAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTcc                >  2:2555517/1‑150 (MQ=255)
                                                                                                                                          ggTAGGGATATAAAATTGGATCATGTCGGGTAACTTTCTGTTGTCAGTTAGATTTCCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCATAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCg             <  1:4632096/149‑1 (MQ=255)
                                                                                                                                          cGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCAAAATACCCGGCACTTTTGCAATTTCCTCg             <  2:118309/150‑1 (MQ=255)
                                                                                                                                          cGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCAAAATACCCGGCACTTTTGCAATTTCCTCg             >  1:1011667/1‑150 (MQ=255)
                                                                                                                                          cGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCg             >  1:1840701/1‑150 (MQ=255)
                                                                                                                                          cGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCg             <  2:2719065/150‑1 (MQ=255)
                                                                                                                                            tAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGAc           <  1:2580726/150‑1 (MQ=255)
                                                                                                                                            tAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGAc           <  1:2581158/150‑1 (MQ=255)
                                                                                                                                            tAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGAc           >  2:851560/1‑150 (MQ=255)
                                                                                                                                             aGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACg          <  2:4536426/150‑1 (MQ=255)
                                                                                                                                               ggATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCAAAATACCCGGCACTTTTGCAATTTCCTCGACGCt        >  1:2987535/1‑150 (MQ=255)
                                                                                                                                               ggATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCAAAATACCCGGCACTTTTGCAATTTCCTCGACGCt        >  1:3013400/1‑150 (MQ=255)
                                                                                                                                                 atatAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTgg      >  1:2054135/1‑150 (MQ=255)
                                                                                                                                                     tAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTGCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGCGtt  <  2:3720861/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGCGTT  >  CP009273/1986140‑1986438

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: