Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,508,588 C→A 6.4% D229Y (GAT→TAT)  dam ← DNA adenine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,508,5880CA 6.4% 932.5 / 27.0 248D229Y (GAT→TAT) damDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (98/134);  new base A (3/13);  total (101/147)
Fisher's exact test for biased strand distribution p-value = 7.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 264 total aligned reads displayed.

TTTTTTCGCGGGTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGCAAT  >  CP009273/3508439‑3508737
                                                                                                                                                     |                                                                                                                                                     
 tttttCGCGGGTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACatcgtatcg                                                                                                                                                      <  2:2470145/150‑1 (MQ=255)
 tttttCGCGGGTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACatcgtatcg                                                                                                                                                      <  2:2470143/150‑1 (MQ=255)
  ttttCGCGGGTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACatcgtatagt                                                                                                                                                     <  1:2648755/150‑1 (MQ=255)
     tCGCGGGTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGAt                                                                                                                                                  <  2:1365325/150‑1 (MQ=255)
     tCGCGGGTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGAt                                                                                                                                                  >  2:3501758/1‑150 (MQ=255)
     tCGCGGGTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGAt                                                                                                                                                  <  1:2060789/150‑1 (MQ=255)
     tCGCGGGTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCATGAt                                                                                                                                                  <  2:3584141/150‑1 (MQ=255)
        cgGGTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTgg                                                                                                                                               >  1:1888999/1‑150 (MQ=255)
        cgGGTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTgg                                                                                                                                               <  2:4313702/150‑1 (MQ=255)
          ggTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGgag                                                                                                                                             <  1:301032/150‑1 (MQ=255)
           gTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGGATCGTGATTGgaga                                                                                                                                            <  2:2123538/150‑1 (MQ=255)
             gAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATc                                                                                                                                          >  2:2709120/1‑150 (MQ=255)
                 cGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCa                                                                                                                                      <  1:215263/150‑1 (MQ=255)
                  gACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCAc                                                                                                                                     >  1:3778489/1‑150 (MQ=255)
                     tCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTgg                                                                                                                                  <  1:4224218/150‑1 (MQ=255)
                      ccTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGa                                                                                                                                 <  1:4342524/150‑1 (MQ=255)
                         ggTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATc                                                                                                                                          <  2:2199819/138‑1 (MQ=255)
                         ggTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAata                                                                                                                              >  1:3932522/1‑150 (MQ=255)
                          gTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAatat                                                                                                                             >  2:1837836/1‑150 (MQ=255)
                             tGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGgc                                                                                                                          >  1:2263311/1‑150 (MQ=255)
                                aCAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCt                                                                                                                       >  2:977318/1‑150 (MQ=255)
                                   aaGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTTCGTCGAACTTTGACGACATGCAATTTTTCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCaa                                                                                                                    >  1:1799522/1‑150 (MQ=255)
                                      ccAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCa                                                                                                                 <  2:2798823/150‑1 (MQ=255)
                                      ccAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCa                                                                                                                 >  1:2830394/1‑150 (MQ=255)
                                       cagcagTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAg                                                                                                                >  2:3282494/1‑150 (MQ=255)
                                              tCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCg                                                                                                         >  2:1302954/1‑150 (MQ=255)
                                              tCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCg                                                                                                         <  2:1691757/150‑1 (MQ=255)
                                              tCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCg                                                                                                         <  2:363739/150‑1 (MQ=255)
                                              tCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGAGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCg                                                                                                         <  1:3270398/150‑1 (MQ=255)
                                               cGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCgg                                                                                                        >  1:429916/1‑150 (MQ=255)
                                                gTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATAGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGc                                                                                                       >  2:511061/1‑150 (MQ=255)
                                                   cACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGAt                                                                                                    >  1:3961413/1‑150 (MQ=255)
                                                    aCCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGAtc                                                                                                   >  1:2692849/1‑150 (MQ=255)
                                                    aCCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGAtc                                                                                                   >  1:4494125/1‑150 (MQ=255)
                                                     ccTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGAtct                                                                                                  <  1:3563187/150‑1 (MQ=255)
                                                     ccTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGAtct                                                                                                  <  1:926403/150‑1 (MQ=255)
                                                         tttACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGc                                                                                              <  2:2378533/150‑1 (MQ=255)
                                                          ttACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGcc                                                                                             <  1:4080300/150‑1 (MQ=255)
                                                           tACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCa                                                                                            >  2:3833385/1‑150 (MQ=255)
                                                           tACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCa                                                                                            <  1:4311234/150‑1 (MQ=255)
                                                           tACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCa                                                                                            <  2:1314839/150‑1 (MQ=255)
                                                            aCGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAg                                                                                           <  2:3056645/150‑1 (MQ=255)
                                                             cGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGa                                                                                          <  1:678260/150‑1 (MQ=255)
                                                               tgtgCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATg                                                                                        >  1:1444204/1‑150 (MQ=255)
                                                               tgtgCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCACCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATg                                                                                        >  2:3668861/1‑150 (MQ=255)
                                                                 tgCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATgcg                                                                                      <  1:4484848/150‑1 (MQ=255)
                                                                 tgCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATgcg                                                                                      >  1:342433/1‑150 (MQ=255)
                                                                 tgCCGCCGTTGCTGCTGATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACCTCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATgcg                                                                                      <  2:750796/150‑1 (MQ=255)
                                                                     gccgTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCttg                                                                                  >  1:1223088/1‑150 (MQ=255)
                                                                       cgTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATAGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCttgtt                                                                                <  2:2754898/150‑1 (MQ=255)
                                                                        gTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCttgttg                                                                               >  2:2414779/1‑150 (MQ=255)
                                                                          tgctgcTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTgttgtt                                                                             <  1:3490108/150‑1 (MQ=255)
                                                                           gctgctTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTc                                                                            >  1:1367402/1‑150 (MQ=255)
                                                                            ctgctTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATAGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCa                                                                           <  2:2505030/150‑1 (MQ=255)
                                                                            ctgctTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATAGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCa                                                                           <  2:205082/150‑1 (MQ=255)
                                                                               ctTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGc                                                                        >  2:1858963/1‑150 (MQ=255)
                                                                               ctTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATAGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGc                                                                        <  1:34653/150‑1 (MQ=255)
                                                                                    aCTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATAGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTaaa                                                                   <  1:511061/150‑1 (MQ=255)
                                                                                       gCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACt                                                                <  1:2093671/150‑1 (MQ=255)
                                                                                       gCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACt                                                                <  1:4049375/150‑1 (MQ=255)
                                                                                          tCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGtt                                                             >  1:1255880/1‑150 (MQ=255)
                                                                                           cGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGttt                                                            >  2:315282/1‑150 (MQ=255)
                                                                                            gAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTtg                                                           >  1:2703051/1‑150 (MQ=255)
                                                                                              aCTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTtgtg                                                         <  2:4486365/150‑1 (MQ=255)
                                                                                               cTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTtgtgt                                                        <  2:832665/150‑1 (MQ=255)
                                                                                               cTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTtgtgt                                                        <  2:2668293/150‑1 (MQ=255)
                                                                                                tttGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTtgtgtg                                                       >  2:4415167/1‑150 (MQ=255)
                                                                                                tttGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTtgtgtg                                                       <  2:2568689/150‑1 (MQ=255)
                                                                                                tttGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTtgtgtg                                                       <  2:266486/150‑1 (MQ=255)
                                                                                                tttGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTtgtgtg                                                       <  2:1205346/150‑1 (MQ=255)
                                                                                                 ttGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATAGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGa                                                      <  1:1043149/150‑1 (MQ=255)
                                                                                                    atgacATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATAc                                                   <  2:458182/148‑1 (MQ=255)
                                                                                                     cgacATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATAGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACg                                                  >  1:612321/1‑150 (MQ=255)
                                                                                                      gacATGCAATTGTGCGCGATGATACCACTCACGAGTTAACATCGTTTCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGc                                                 <  1:2128229/150‑1 (MQ=255)
                                                                                                       acATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGcc                                                <  2:3778110/150‑1 (MQ=255)
                                                                                                            cAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTaaa                                           <  1:3310080/150‑1 (MQ=255)
                                                                                                               ttttGCGCGCTGATACCACTCTCGCGTTAACATCGTATCGTTATTGGAGATCAGCTCTGGAATATGGCGCTCAACCTGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGtt                                        <  2:2516497/150‑1 (MQ=255)
                                                                                                               ttttGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGtt                                        <  2:2901741/150‑1 (MQ=255)
                                                                                                                tttGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTg                                       <  1:746817/150‑1 (MQ=255)
                                                                                                                 ttGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTgg                                      <  1:1303015/150‑1 (MQ=255)
                                                                                                                  tGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGc                                     >  1:1537760/1‑150 (MQ=255)
                                                                                                                   gcgcgcTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCg                                    <  2:4234339/150‑1 (MQ=255)
                                                                                                                      cgcTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTc                                 <  2:1261695/150‑1 (MQ=255)
                                                                                                                      cgcTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTc                                 <  2:1261473/150‑1 (MQ=255)
                                                                                                                       gcTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCg                                <  1:1995862/150‑1 (MQ=255)
                                                                                                                             aCCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACa                          <  2:4051572/150‑1 (MQ=255)
                                                                                                                                aCTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACAGCg                       >  2:3416693/1‑150 (MQ=255)
                                                                                                                                  tCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTCAAGTTGGCGGTCGCAGACAGCGGt                     >  2:3240972/1‑150 (MQ=255)
                                                                                                                                         ttAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAg              >  2:1583352/1‑150 (MQ=255)
                                                                                                                                          tAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAgg             <  2:2541329/150‑1 (MQ=255)
                                                                                                                                          tAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAgg             <  2:552364/150‑1 (MQ=255)
                                                                                                                                          tAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAgg             <  2:693828/150‑1 (MQ=255)
                                                                                                                                           aaCATCGTATCGTGATTGGATATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGc            >  1:2982069/1‑150 (MQ=255)
                                                                                                                                           aaCATCGTATCGTGATTGGATATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGc            >  1:2984559/1‑150 (MQ=255)
                                                                                                                                           aaCATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCACACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGc            >  2:4390609/1‑150 (MQ=255)
                                                                                                                                            aCATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCg           >  1:244509/1‑150 (MQ=255)
                                                                                                                                              atcgtatcgtGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGa         >  1:3619122/1‑150 (MQ=255)
                                                                                                                                                 gtatcgtGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGAGGATCg      >  1:1839093/1‑150 (MQ=255)
                                                                                                                                                  tatcgtGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGc     >  1:1491759/1‑150 (MQ=255)
                                                                                                                                                     cgtGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGCAAt  >  2:3555208/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTTTTCGCGGGTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGCAAT  >  CP009273/3508439‑3508737

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: