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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,986,2890TA92.1% 944.9 / 53.0 278N4Y (AAT→TAT) pgsAphosphatidylglycerophosphate synthetase
Reads supporting (aligned to +/- strand):  ref base T (10/12);  new base A (111/145);  total (121/157)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 283 total aligned reads displayed.

CCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGCGTT  >  CP009273/1986142‑1986438
                                                                                                                                                   |                                                                                                                                                     
ccAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGCATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTaa                                                                                                                                                     <  2:1033995/150‑1 (MQ=255)
  aGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAAt                                                                                                                                                   >  2:295824/1‑150 (MQ=255)
   gAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAAtt                                                                                                                                                  <  1:2966091/150‑1 (MQ=255)
    aaaaCCATCGAACCAGTCAGTCACCGCCGCGACGCTGTAAATGAGCGCGGCGGCAAACGTCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTg                                                                                                                                                 <  1:4533896/150‑1 (MQ=255)
    aaaaCCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTg                                                                                                                                                 >  1:4145515/1‑150 (MQ=255)
     aaaCCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGACTGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGc                                                                                                                                                >  2:4218963/1‑150 (MQ=255)
       aCCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCAt                                                                                                                                              <  2:3449411/150‑1 (MQ=255)
        ccATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATa                                                                                                                                             <  2:4722069/150‑1 (MQ=255)
          aTCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAAt                                                                                                                                           <  1:2731706/150‑1 (MQ=255)
              aaCCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACg                                                                                                                                       <  1:4900312/150‑1 (MQ=255)
              aaCCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACg                                                                                                                                       <  2:1713398/150‑1 (MQ=255)
                  agtcagtcACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATTAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTa                                                                                                                                   <  2:1111421/150‑1 (MQ=255)
                  agtcagtcACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTa                                                                                                                                   <  2:4805731/150‑1 (MQ=255)
                   gtcagtcaCCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTaa                                                                                                                                  <  2:3685330/150‑1 (MQ=255)
                   gtcagtcaCCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTaa                                                                                                                                  >  1:2132884/1‑150 (MQ=255)
                     cagtcaCCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACt                                                                                                                                <  2:346259/150‑1 (MQ=255)
                            cgccgcGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCtgtt                                                                                                                         <  2:3837683/150‑1 (MQ=255)
                             gccgcGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCtgttg                                                                                                                        <  1:1942177/150‑1 (MQ=255)
                             gccgcGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCtgttg                                                                                                                        <  1:3498746/150‑1 (MQ=255)
                               cgcgACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTc                                                                                                                      <  2:3420513/150‑1 (MQ=255)
                                  gACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGt                                                                                                                   >  1:3236446/1‑150 (MQ=255)
                                  gACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGt                                                                                                                   <  1:3632769/150‑1 (MQ=255)
                                    cGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTTACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTaa                                                                                                                 <  2:4825530/150‑1 (MQ=255)
                                     gCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAg                                                                                                                <  1:3572058/150‑1 (MQ=255)
                                        gAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGAtt                                                                                                             <  1:4388613/150‑1 (MQ=255)
                                         aaaaTGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTa                                                                                                            <  1:2550479/150‑1 (MQ=255)
                                              gAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCt                                                                                                       <  1:2518825/150‑1 (MQ=255)
                                              gAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCt                                                                                                       <  2:1287363/150‑1 (MQ=255)
                                              gAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCt                                                                                                       >  1:3985008/1‑150 (MQ=255)
                                              gAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACAAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCt                                                                                                       <  1:874208/150‑1 (MQ=255)
                                                 cgcgGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATg                                                                                                    >  1:1680661/1‑150 (MQ=255)
                                                 cgcgGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATg                                                                                                    <  1:1465802/150‑1 (MQ=255)
                                                 cgcgGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATg                                                                                                    <  2:1523913/150‑1 (MQ=255)
                                                       ggcAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTa                                                                                              >  1:2227026/1‑150 (MQ=255)
                                                        gcAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTAc                                                                                             <  1:1480018/150‑1 (MQ=255)
                                                            aCGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACagag                                                                                         <  2:2717235/150‑1 (MQ=255)
                                                            aCGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACagag                                                                                         <  2:2529252/150‑1 (MQ=255)
                                                            aCGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACagag                                                                                         <  2:2029826/150‑1 (MQ=255)
                                                             cGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAgaga                                                                                        <  2:3306539/150‑1 (MQ=255)
                                                                 gACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATc                                                                                    <  2:4906910/150‑1 (MQ=255)
                                                                  aCCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCa                                                                                   >  1:947723/1‑150 (MQ=255)
                                                                    cAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAAt                                                                                 <  2:1865357/150‑1 (MQ=255)
                                                                     aGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATg                                                                                >  2:701627/1‑150 (MQ=255)
                                                                     aGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATg                                                                                >  2:701731/1‑150 (MQ=255)
                                                                      ggTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGt                                                                               <  1:3560302/150‑1 (MQ=255)
                                                                        tGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGttt                                                                             >  2:658461/1‑150 (MQ=255)
                                                                           cAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCaa                                                                          >  2:4997996/1‑150 (MQ=255)
                                                                           cAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCaa                                                                          >  1:1311183/1‑150 (MQ=255)
                                                                            aGGCATATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGTACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGTCGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAAc                                                                         <  2:628939/150‑1 (MQ=255)
                                                                            aGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAAc                                                                         <  1:2663283/150‑1 (MQ=255)
                                                                              gCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGa                                                                       >  2:4359508/1‑150 (MQ=255)
                                                                               cAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGAc                                                                      >  1:1349047/1‑150 (MQ=255)
                                                                                aGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGAcc                                                                     <  2:838452/150‑1 (MQ=255)
                                                                                aGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGAcc                                                                     >  2:3561950/1‑150 (MQ=255)
                                                                                  aTAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCag                                                                   <  1:2464301/150‑1 (MQ=255)
                                                                                  aTAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCag                                                                   <  2:2776252/150‑1 (MQ=255)
                                                                                    aaaaGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCagaa                                                                 >  1:4591465/1‑150 (MQ=255)
                                                                                    aaaaGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCagaa                                                                 >  1:3814203/1‑150 (MQ=255)
                                                                                      aaGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCagaaga                                                               >  1:4293206/1‑150 (MQ=255)
                                                                                        gACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATc                                                             <  2:4550118/150‑1 (MQ=255)
                                                                                         aCCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCt                                                            >  2:5057337/1‑150 (MQ=255)
                                                                                           cAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCttt                                                          >  1:1496927/1‑150 (MQ=255)
                                                                                           cAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCttt                                                          <  1:4567894/150‑1 (MQ=255)
                                                                                               aCAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTg                                                      <  2:2240470/150‑1 (MQ=255)
                                                                                                cAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGc                                                     >  2:753284/1‑150 (MQ=255)
                                                                                                 aaaGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGcc                                                    >  2:2005516/1‑150 (MQ=255)
                                                                                                 aaaGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGcc                                                    >  2:396921/1‑150 (MQ=255)
                                                                                                 aaaGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGcc                                                    >  2:4603258/1‑150 (MQ=255)
                                                                                                   aGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAg                                                  <  1:473752/150‑1 (MQ=255)
                                                                                                   aGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAg                                                  <  1:186803/150‑1 (MQ=255)
                                                                                                      aTGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGAcc                                               <  2:4137040/150‑1 (MQ=255)
                                                                                                        gggATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATACTGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCtt                                             <  2:3113005/150‑1 (MQ=255)
                                                                                                         ggATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTg                                            <  1:4224230/150‑1 (MQ=255)
                                                                                                         ggATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTg                                            <  1:5125096/150‑1 (MQ=255)
                                                                                                           aTAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCg                                          <  2:3508166/150‑1 (MQ=255)
                                                                                                            tAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGa                                         <  2:1528521/150‑1 (MQ=255)
                                                                                                               ggATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATACTGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAAt                                      <  2:1419957/150‑1 (MQ=255)
                                                                                                               ggATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAAt                                      >  1:1609773/1‑150 (MQ=255)
                                                                                                               ggATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAAt                                      <  1:2293256/150‑1 (MQ=255)
                                                                                                                gATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATa                                     <  2:975217/150‑1 (MQ=255)
                                                                                                                gATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATa                                     >  1:750427/1‑150 (MQ=255)
                                                                                                                 aTGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATAc                                    >  1:1100082/1‑150 (MQ=255)
                                                                                                                 aTGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATAc                                    >  1:4050099/1‑150 (MQ=255)
                                                                                                                   gACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATAccc                                  <  1:3298132/150‑1 (MQ=255)
                                                                                                                        ggAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCAc                             <  2:2112000/150‑1 (MQ=255)
                                                                                                                         gAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACt                            <  2:4672377/150‑1 (MQ=255)
                                                                                                                         gAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACt                            <  2:1896553/150‑1 (MQ=255)
                                                                                                                            cAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACtttt                         <  2:4303438/150‑1 (MQ=255)
                                                                                                                            cAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACtttt                         >  1:2216156/1‑150 (MQ=255)
                                                                                                                                   aGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAAttt                  <  2:2155990/150‑1 (MQ=255)
                                                                                                                                    gCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTc                 <  1:1152842/150‑1 (MQ=255)
                                                                                                                                     cAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTcc                >  2:2035354/1‑150 (MQ=255)
                                                                                                                                        cGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCg             <  2:4563138/150‑1 (MQ=255)
                                                                                                                                         gTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGa            <  2:2363796/150‑1 (MQ=255)
                                                                                                                                         gTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGa            >  1:3042397/1‑150 (MQ=255)
                                                                                                                                          tAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGAc           >  2:815421/1‑150 (MQ=255)
                                                                                                                                             ggATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCt        >  2:3944816/1‑150 (MQ=255)
                                                                                                                                             ggATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCt        <  1:642574/150‑1 (MQ=255)
                                                                                                                                               atatAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTgg      >  2:2533659/1‑150 (MQ=255)
                                                                                                                                                   tAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGCGtt  >  1:2368963/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGCGTT  >  CP009273/1986142‑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: