Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,615,370 C→A 53.0% S272Y (TCT→TAT)  purM → phosphoribosylaminoimidazole synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,615,3700CA53.0% 67.8 / inf 296S272Y (TCT→TAT) purMphosphoribosylaminoimidazole synthetase
Reads supporting (aligned to +/- strand):  ref base C (91/48);  new base A (92/65);  total (183/113)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

Only 100 of 297 total aligned reads displayed.

CTGGCACCGACCCGCATTTACGTGAAGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATG  >  CP009273/2615222‑2615519
                                                                                                                                                    |                                                                                                                                                     
 tgGCACCGACCCGCATTTACGTGAAGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCtt                                                                                                                                                     >  2:808424/1‑150 (MQ=255)
 tgGCACCGACCCGCATTTACGTGAAGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCtt                                                                                                                                                     >  2:4064045/1‑150 (MQ=255)
   gCACCGACCCGCATTTACGTGAAGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTcc                                                                                                                                                   >  1:1976942/1‑150 (MQ=255)
    cACCGACCCGCATTTACGTGAAGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCt                                                                                                                                                  >  2:4124259/1‑150 (MQ=255)
    cACCGACCCGCATTTACGTGAAGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCt                                                                                                                                                  <  1:1016609/150‑1 (MQ=255)
       cGACCCGCATTTACGTGAAGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGc                                                                                                                                               >  1:210090/1‑150 (MQ=255)
        gACCCGCATTTACGTGAAGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCa                                                                                                                                              >  2:4245449/1‑150 (MQ=255)
         aCCCGCATTTACGTGAAGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAg                                                                                                                                             <  2:3382391/150‑1 (MQ=255)
              cATTTACGTGAAGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGcc                                                                                                                                        >  1:4065578/1‑150 (MQ=255)
               aTTTACGTGAAGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCg                                                                                                                                       >  2:1835781/1‑150 (MQ=255)
                 ttACGTGAAGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGa                                                                                                                                     >  2:4523387/1‑150 (MQ=255)
                      tGAAGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAgtgt                                                                                                                                >  2:4496342/1‑150 (MQ=255)
                        aaGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTc                                                                                                                              >  1:1349196/1‑150 (MQ=255)
                         aGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCa                                                                                                                             >  1:4060756/1‑150 (MQ=255)
                         aGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATAGGCCTGGCAGTGGCCGGAAGTGTTCa                                                                                                                             <  2:1190442/150‑1 (MQ=255)
                         aGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACATGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCCGATAATACTCAGGCAGTGAGTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCa                                                                                                                             >  1:21553/1‑150 (MQ=255)
                                  tGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGc                                                                                                                    >  1:59310/1‑150 (MQ=255)
                                  tGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGc                                                                                                                    >  2:2320759/1‑150 (MQ=255)
                                  tGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGc                                                                                                                    >  2:1791751/1‑150 (MQ=255)
                                    gAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCaa                                                                                                                  >  2:636370/1‑150 (MQ=255)
                                        tgattgaAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACgg                                                                                                              <  2:3303443/150‑1 (MQ=255)
                                          attgaAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCa                                                                                                            >  2:572877/1‑150 (MQ=255)
                                          attgaAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCa                                                                                                            >  1:619724/1‑150 (MQ=255)
                                           ttgaAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAg                                                                                                           <  2:2013321/150‑1 (MQ=255)
                                           ttgaAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAg                                                                                                           <  1:3495540/150‑1 (MQ=255)
                                            tgaAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAgg                                                                                                          >  1:146646/1‑150 (MQ=255)
                                               aaaGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTATGGGAAAACATTCCGCGCGTATTGCCAGATAATAATCAGTCAGTGATTGATGAATCTTCCTGGCAGTGGCAGGAAGTGTTCAACTGGCTGCAAACGGCAGGTaa                                                                                                       >  1:1544291/1‑150 (MQ=255)
                                                aaGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAAc                                                                                                      >  2:3735324/1‑150 (MQ=255)
                                                     cGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTTTTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGa                                                                                                 <  1:3553777/150‑1 (MQ=255)
                                                     cGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGa                                                                                                 >  1:2668521/1‑150 (MQ=255)
                                                     cGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGa                                                                                                 <  1:3824236/150‑1 (MQ=255)
                                                      gATTTGCATTCCTTTGCGCTCCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTTATTGTTGTATATTCCTGGCAGTGGCCGGATGTTTTCAAGTGGCTGCAAACGGCAGGTTACGTTgtt                                                                                                <  2:487792/150‑3 (MQ=255)
                                                        tgtgCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCa                                                                                              >  1:3569100/1‑150 (MQ=255)
                                                        tgtgCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATAGTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCa                                                                                              <  2:345381/150‑1 (MQ=255)
                                                         gtgCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGACTAGTCCGGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCAc                                                                                             <  2:1722478/150‑1 (MQ=255)
                                                             aTGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATg                                                                                         >  1:3968284/1‑150 (MQ=255)
                                                             aTGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATg                                                                                         >  2:240541/1‑150 (MQ=255)
                                                              tGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGa                                                                                        >  1:4112867/1‑150 (MQ=255)
                                                                 cATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAAt                                                                                     >  1:17457/1‑150 (MQ=255)
                                                                 cATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAAt                                                                                     >  2:2750897/1‑150 (MQ=255)
                                                                 cATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGGATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAAt                                                                                     <  1:4595820/150‑1 (MQ=255)
                                                                  aTTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATg                                                                                    >  2:2846649/1‑150 (MQ=255)
                                                                  aTTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATg                                                                                    >  2:1286519/1‑150 (MQ=255)
                                                                      cgcACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATc                                                                                <  1:1655556/150‑1 (MQ=255)
                                                                         aCCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCa                                                                             >  1:4585665/1‑150 (MQ=255)
                                                                           cTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCAcc                                                                           <  1:1176893/150‑1 (MQ=255)
                                                                           cTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCAcc                                                                           <  1:3399410/150‑1 (MQ=255)
                                                                            tGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCt                                                                          >  2:2453189/1‑150 (MQ=255)
                                                                             gACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCtt                                                                         >  1:451788/1‑150 (MQ=255)
                                                                             gACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCtt                                                                         >  1:1639535/1‑150 (MQ=255)
                                                                             gACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCtt                                                                         <  2:1175454/150‑1 (MQ=255)
                                                                             gACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATGGCCAGATAATACTCAGGCAGTGATTGATGAATAGTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCtt                                                                         <  2:223324/150‑1 (MQ=255)
                                                                              aCCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTc                                                                        <  1:1315934/150‑1 (MQ=255)
                                                                              aCCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTc                                                                        <  1:2374617/150‑1 (MQ=255)
                                                                              aCCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTc                                                                        <  1:2975627/150‑1 (MQ=255)
                                                                                 ggcggcggcTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAAc                                                                     >  2:429537/1‑150 (MQ=255)
                                                                                  gcggcggcTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACt                                                                    <  1:2930557/150‑1 (MQ=255)
                                                                                      cggcTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTgcgg                                                                >  1:2429787/1‑150 (MQ=255)
                                                                                      cggcTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTgcgg                                                                >  1:505152/1‑150 (MQ=255)
                                                                                       ggcTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTgcggc                                                               >  1:326316/1‑150 (MQ=255)
                                                                                       ggcTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTgcggc                                                               >  2:3351570/1‑150 (MQ=255)
                                                                                            cTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATAGTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCgg                                                          <  2:1580608/150‑1 (MQ=255)
                                                                                             tGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCggg                                                         >  1:1829595/1‑150 (MQ=255)
                                                                                              gggAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGga                                                        >  1:3618745/1‑150 (MQ=255)
                                                                                              gggAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGga                                                        >  1:1870939/1‑150 (MQ=255)
                                                                                               ggAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGgat                                                       >  1:4364464/1‑150 (MQ=255)
                                                                                               ggAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGgat                                                       >  2:1555328/1‑150 (MQ=255)
                                                                                                    aCATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGatta                                                  >  1:2048692/1‑150 (MQ=255)
                                                                                                      aTTCCTCGCGTATTGCCAGATAATACTCATGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGATGTGTTCAACTGGCTGCAAACGGCATGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAAATTCGGCGTCGGGATGattatt                                                <  2:1405014/150‑1 (MQ=255)
                                                                                                      aTTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGattatt                                                >  1:4271266/1‑150 (MQ=255)
                                                                                                      aTTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGattatt                                                >  2:1881787/1‑150 (MQ=255)
                                                                                                           gcgcgTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCcct                                           >  1:2268534/1‑150 (MQ=255)
                                                                                                           gcgcgTATTGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCcct                                           <  2:4408162/150‑1 (MQ=255)
                                                                                                            cgcgTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCcctg                                          >  2:458228/1‑150 (MQ=255)
                                                                                                                   tGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTc                                   <  1:3446722/150‑1 (MQ=255)
                                                                                                                   tGCCAGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTc                                   <  2:1218700/150‑1 (MQ=255)
                                                                                                                   tGCCAGATAATACTCAGGCAGTGATTGATGAATAGTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTc                                   <  2:212421/150‑1 (MQ=255)
                                                                                                                       aGATAATACTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGa                               >  1:3202629/1‑150 (MQ=255)
                                                                                                                         ataataCTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAg                             >  1:2835654/1‑150 (MQ=255)
                                                                                                                             taCTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGa                         >  2:3107859/1‑150 (MQ=255)
                                                                                                                             taCTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGa                         <  2:3582154/150‑1 (MQ=255)
                                                                                                                               cTCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACa                       >  1:2181282/1‑150 (MQ=255)
                                                                                                                                tCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACaa                      >  1:660949/1‑150 (MQ=255)
                                                                                                                                tCAGGCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACaa                      >  1:1900975/1‑150 (MQ=255)
                                                                                                                                 cAGGCAGTGATTGATGACGCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACaaa                     <  2:2600961/150‑1 (MQ=255)
                                                                                                                                 cAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACaaa                     >  1:4202581/1‑150 (MQ=255)
                                                                                                                                   ggCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGc                   >  2:3889204/1‑150 (MQ=255)
                                                                                                                                   ggCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGc                   >  2:2425827/1‑150 (MQ=255)
                                                                                                                                   ggCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGc                   >  1:4326752/1‑150 (MQ=255)
                                                                                                                                   ggCAGTGATTGATGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGc                   >  1:2743857/1‑150 (MQ=255)
                                                                                                                                     cAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGccc                 >  2:2933661/1‑150 (MQ=255)
                                                                                                                                     cAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGccc                 >  1:2892686/1‑150 (MQ=255)
                                                                                                                                       gTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTc               >  1:4624465/1‑150 (MQ=255)
                                                                                                                                       gTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTc               >  2:934968/1‑150 (MQ=255)
                                                                                                                                          attgatGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGcc            >  2:3517228/1‑150 (MQ=255)
                                                                                                                                           ttgatGAATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGccc           <  2:3991211/150‑1 (MQ=255)
                                                                                                                                            tgatgaATATTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCt          <  1:2507495/150‑1 (MQ=255)
                                                                                                                                              atgaATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGc        >  1:372429/1‑150 (MQ=255)
                                                                                                                                                gaATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCACGCCCTGCTc      <  2:122858/150‑1 (MQ=255)
                                                                                                                                                   tCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAAt   >  2:2164698/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGGCACCGACCCGCATTTACGTGAAGTCAGTGCTGGAGTTGATTGAAAAGGTCGATGTGCATGCCATTGCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATG  >  CP009273/2615222‑2615519

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: