Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,902,302 C→A 73.5% D242Y (GAT→TAT)  pyrG ← CTP synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,902,3020CA73.5% 379.7 / 242.5 287D242Y (GAT→TAT) pyrGCTP synthetase
Reads supporting (aligned to +/- strand):  ref base C (42/34);  new base A (86/125);  total (128/159)
Fisher's exact test for biased strand distribution p-value = 3.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 287 total aligned reads displayed.

GACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGC  >  CP009273/2902156‑2902450
                                                                                                                                                  |                                                                                                                                                    
  cTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGt                                                                                                                                                 >  2:870973/1‑150 (MQ=255)
  cTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGt                                                                                                                                                 <  1:3719970/150‑1 (MQ=255)
    tCACTTACCGTGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACATATTCGCTTCGGGTCATTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCt                                                                                                                                               <  1:1945725/150‑1 (MQ=255)
    tCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCt                                                                                                                                               >  2:915128/1‑150 (MQ=255)
        ttACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCa                                                                                                                                           >  1:2682665/1‑150 (MQ=255)
         tACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCag                                                                                                                                          >  2:3136141/1‑150 (MQ=255)
           ccGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCagag                                                                                                                                        <  1:4250433/150‑1 (MQ=255)
            cGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAgaga                                                                                                                                       <  1:4177507/150‑1 (MQ=255)
             gggTTCGCGACTTCGAAGATAACCTGTTCCCATTCGGACCGATACGCTTCCGGGCAGTTTAAGCTGAATCGTGTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGACCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGaa                                                                                                                                      <  2:848923/150‑1 (MQ=255)
              ggTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGaaa                                                                                                                                     <  1:3343038/150‑1 (MQ=255)
               gTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAAt                                                                                                                                    >  1:4712608/1‑150 (MQ=255)
                 tCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATAGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATaa                                                                                                                                  >  1:1501123/1‑150 (MQ=255)
                  cGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAAc                                                                                                                                 >  2:4193864/1‑150 (MQ=255)
                   gCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAAcc                                                                                                                                >  1:3510569/1‑150 (MQ=255)
                   gCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAAcc                                                                                                                                <  1:2623197/150‑1 (MQ=255)
                   gCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCGGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAAcc                                                                                                                                >  1:1014070/1‑150 (MQ=255)
                      tcttcGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCt                                                                                                                             <  2:4213852/150‑1 (MQ=255)
                       cttcGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCGGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCtt                                                                                                                            >  1:4554103/1‑150 (MQ=255)
                         tcGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCtttt                                                                                                                          <  1:3932623/150‑1 (MQ=255)
                          cGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCttttt                                                                                                                         <  1:2196615/150‑1 (MQ=255)
                            aaGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTcc                                                                                                                       >  2:5257451/1‑150 (MQ=255)
                             aGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCg                                                                                                                      <  2:3793054/150‑1 (MQ=255)
                              gATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCggg                                                                                                                     >  1:4634620/1‑150 (MQ=255)
                              gATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCgg                                                                                                                     >  2:2193516/1‑150 (MQ=255)
                               aTAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGa                                                                                                                    <  2:4867412/150‑1 (MQ=255)
                               aTAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGa                                                                                                                    <  2:3224660/150‑1 (MQ=255)
                                 aaCCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAAc                                                                                                                  >  2:2750114/1‑150 (MQ=255)
                                   ccTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACAt                                                                                                                <  2:5153494/150‑1 (MQ=255)
                                   ccTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACAt                                                                                                                <  2:2844718/150‑1 (MQ=255)
                                    cTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACAtt                                                                                                               >  2:1856324/1‑150 (MQ=255)
                                    cTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACAtt                                                                                                               <  2:4164897/150‑1 (MQ=255)
                                        tCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAg                                                                                                           >  1:2855786/1‑150 (MQ=255)
                                         cacATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTAAAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGa                                                                                                          <  1:4448812/148‑1 (MQ=255)
                                          ccATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGaa                                                                                                         >  2:2914238/1‑150 (MQ=255)
                                           cATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAAc                                                                                                        >  2:5068658/1‑150 (MQ=255)
                                           cATTCGGACAGAGTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACACGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAAc                                                                                                        <  2:512051/150‑1 (MQ=255)
                                            aTTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACa                                                                                                       >  1:2581485/1‑150 (MQ=255)
                                               cGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATg                                                                                                    <  1:5084093/150‑1 (MQ=255)
                                                ggACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAAttc                                                                                                   >  1:4592956/1‑148 (MQ=255)
                                                      aTTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCtt                                                                                             <  2:261663/150‑1 (MQ=255)
                                                        tCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCg                                                                                           >  1:3104825/1‑150 (MQ=255)
                                                        tCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCg                                                                                           <  1:298781/150‑1 (MQ=255)
                                                         cGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGc                                                                                          >  1:688183/1‑150 (MQ=255)
                                                         cGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGc                                                                                          >  2:2024867/1‑150 (MQ=255)
                                                              ccGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCAcgtt                                                                                     <  2:4584500/150‑1 (MQ=255)
                                                              ccGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCAcgtt                                                                                     <  1:4288972/150‑1 (MQ=255)
                                                                gggCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCAcgttcg                                                                                   <  2:5283051/150‑1 (MQ=255)
                                                                 ggCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCAcgttcgt                                                                                  >  2:1240331/1‑150 (MQ=255)
                                                                 ggCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCAcgttcgt                                                                                  >  2:3582095/1‑150 (MQ=255)
                                                                  gCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCAcgttcgtt                                                                                 <  1:1240331/150‑1 (MQ=255)
                                                                    aGTTTAAGCTGAATCGTTTACAAATATACTCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGttcgttcg                                                                               <  2:1577837/150‑1 (MQ=255)
                                                                     gTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATc                                                                                               <  2:3060177/133‑1 (MQ=255)
                                                                     gTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGc                                                                              <  1:95149/150‑1 (MQ=255)
                                                                      tttAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGcc                                                                             <  1:2024867/150‑1 (MQ=255)
                                                                      tttAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGcc                                                                             <  2:582299/150‑1 (MQ=255)
                                                                        tAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCgg                                                                           >  2:1117691/1‑150 (MQ=255)
                                                                               aaTCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGTAATCTTCGCACGTTCGTTCGCCGGAACAgcg                                                                    <  1:1082501/150‑1 (MQ=255)
                                                                                 tCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCAGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGctcg                                                                  <  2:1035819/150‑1 (MQ=255)
                                                                                   gTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGAt                                                                <  2:2280452/150‑1 (MQ=255)
                                                                                       aCAAATATAATCGTCCAGCCCCGGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGa                                                            >  1:3184888/1‑150 (MQ=255)
                                                                                        cAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGaa                                                           <  1:906843/150‑1 (MQ=255)
                                                                                         aaaTATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAAc                                                          >  1:773528/1‑150 (MQ=255)
                                                                                          aaTATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACg                                                         >  2:3462180/1‑150 (MQ=255)
                                                                                          aaTATAATCGACCAGCCCCTGAGATTTCAACCGGCCCGGAATTTAAGAGATGGAATAGACGTCGATCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTGCGCCCGTTCGTACGCCGGACCAGCGCGAGCTgcccg                                                         >  1:1343083/1‑146 (MQ=255)
                                                                                             ataATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACa                                                      <  2:748186/150‑1 (MQ=255)
                                                                                             ataATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACa                                                      >  2:2532519/1‑150 (MQ=255)
                                                                                             ataATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACa                                                      >  1:4202865/1‑150 (MQ=255)
                                                                                                  cGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCa                                                 >  2:388867/1‑150 (MQ=255)
                                                                                                   gTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAg                                                <  2:2302858/150‑1 (MQ=255)
                                                                                                   gTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAg                                                >  2:3686845/1‑150 (MQ=255)
                                                                                                   gTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAg                                                <  1:4521112/150‑1 (MQ=255)
                                                                                                     ccAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGa                                              <  1:3304109/150‑1 (MQ=255)
                                                                                                      cAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGAt                                             <  2:1687057/150‑1 (MQ=255)
                                                                                                      cAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGAt                                             <  1:4769654/150‑1 (MQ=255)
                                                                                                      cAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGAt                                             <  1:2108343/150‑1 (MQ=255)
                                                                                                       aGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATg                                            <  1:2414504/150‑1 (MQ=255)
                                                                                                        gCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGt                                           <  2:3283778/150‑1 (MQ=255)
                                                                                                          cccTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCa                                         <  2:2104200/150‑1 (MQ=255)
                                                                                                                 aTTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGa                                  <  1:1667446/150‑1 (MQ=255)
                                                                                                                  tttCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGAt                                 <  1:5068658/150‑1 (MQ=255)
                                                                                                                    tCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATAc                               >  2:3758355/1‑150 (MQ=255)
                                                                                                                    tCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATAc                               >  1:1583505/1‑150 (MQ=255)
                                                                                                                     cAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATAcc                              >  1:4341885/1‑150 (MQ=255)
                                                                                                                       aCAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGa                            >  2:1018807/1‑150 (MQ=255)
                                                                                                                        cAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGAt                           <  1:1763073/150‑1 (MQ=255)
                                                                                                                        cAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGAt                           <  1:2275576/150‑1 (MQ=255)
                                                                                                                        cAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGAt                           <  1:5237427/150‑1 (MQ=255)
                                                                                                                              cGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGagag                     <  2:4004033/150‑1 (MQ=255)
                                                                                                                              cGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGagag                     <  2:1338834/150‑1 (MQ=255)
                                                                                                                                     ttATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGctct              <  1:1706341/150‑1 (MQ=255)
                                                                                                                                      tataGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCtt             >  1:3395340/1‑150 (MQ=255)
                                                                                                                                       ataGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCttt            >  2:1756575/1‑150 (MQ=255)
                                                                                                                                       ataGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCttt            >  1:1275748/1‑150 (MQ=255)
                                                                                                                                       ataGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCttt            >  1:2339942/1‑150 (MQ=255)
                                                                                                                                       ataGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCttt            >  2:4829183/1‑150 (MQ=255)
                                                                                                                                             ggAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACaga      >  1:3848200/1‑150 (MQ=255)
                                                                                                                                              gAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACagag     >  2:4473037/1‑150 (MQ=255)
                                                                                                                                              gAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACagag     >  1:142136/1‑150 (MQ=255)
                                                                                                                                              gAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACagag     >  1:4977556/1‑150 (MQ=255)
                                                                                                                                                 tCGACGTCTTTCAGAGAAATTACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGc  <  1:4059011/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGC  >  CP009273/2902156‑2902450

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: