Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,557,815 T→A 17.7% K709N (AAA→AAT glgP ← glycogen phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,557,8150TA17.7% 777.0 / 151.6 300K709N (AAA→AATglgPglycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (111/136);  new base A (31/22);  total (142/158)
Fisher's exact test for biased strand distribution p-value = 9.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 305 total aligned reads displayed.

CCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTT  >  CP009273/3557667‑3557964
                                                                                                                                                    |                                                                                                                                                     
ccGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGtt                                                                                                                                                      >  1:4580112/1‑150 (MQ=255)
 cGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTAGCGTCAGCACCTGATGCAGCTCCTCATCATTCTCGTAGTATTCACGCGGtga                                                                                                                                                     >  1:1304011/1‑148 (MQ=255)
  gCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTg                                                                                                                                                    >  1:1039864/1‑150 (MQ=255)
   ccAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGt                                                                                                                                                   >  1:654308/1‑150 (MQ=255)
    cAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTa                                                                                                                                                  <  1:2725007/150‑1 (MQ=255)
         ccTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCt                                                                                                                                             <  2:2355448/150‑1 (MQ=255)
         ccTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCt                                                                                                                                             <  1:1741564/150‑1 (MQ=255)
         ccTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCt                                                                                                                                             <  2:45000/150‑1 (MQ=255)
          cTGGTAGTGATCGCCGAATTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCATCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTg                                                                                                                                            <  1:4630069/150‑1 (MQ=255)
           tGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTga                                                                                                                                           <  1:559435/150‑1 (MQ=255)
               aGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTgacgac                                                                                                                                       <  2:1609897/150‑1 (MQ=255)
                 tGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGc                                                                                                                                     >  2:412116/1‑150 (MQ=255)
                  gATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCa                                                                                                                                    >  2:1132693/1‑150 (MQ=255)
                   aTCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAg                                                                                                                                   >  2:30662/1‑150 (MQ=255)
                     cGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGtt                                                                                                                                 <  2:2831098/150‑1 (MQ=255)
                        cGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGttctt                                                                                                                              <  2:799380/150‑1 (MQ=255)
                          aaGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTcc                                                                                                                            <  2:237332/150‑1 (MQ=255)
                           aGTTGATCAGCGAATCAACCAGATCGCTATATCGACCCGGATCTTCCGGACTGAAAACACCGAAGCCGATTAGCGTCAGCACCTGATGCAGCTCCACATCATTCTCTTAGAATACACGCGGTATTTAGCAATGACGACGCAGatcttcaa                                                                                                                           >  1:1158724/1‑148 (MQ=255)
                           aGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCa                                                                                                                           <  1:4327969/150‑1 (MQ=255)
                             ttGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCAct                                                                                                                         <  1:2765346/150‑1 (MQ=255)
                             ttGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCAct                                                                                                                         <  1:2771737/150‑1 (MQ=255)
                             ttGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCAct                                                                                                                         >  1:2841578/1‑150 (MQ=255)
                              tGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCActt                                                                                                                        <  1:2372668/150‑1 (MQ=255)
                              tGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCActt                                                                                                                        <  2:1208323/150‑1 (MQ=255)
                              tGATCAGCGAATCAACCAGATCGCGATAGCGAACCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCTTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACTCAGTTCTTCCAatt                                                                                                                        >  2:23401/1‑150 (MQ=255)
                               gATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCActtc                                                                                                                       <  2:4382786/150‑1 (MQ=255)
                                aTCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCActtct                                                                                                                      >  1:879451/1‑150 (MQ=255)
                                 tCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCActtctt                                                                                                                     >  2:1750578/1‑150 (MQ=255)
                                  cAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACttcttc                                                                                                                    >  2:2821083/1‑150 (MQ=255)
                                    gCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCg                                                                                                                  <  2:2517011/150‑1 (MQ=255)
                                    gCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCg                                                                                                                  <  2:4539899/150‑1 (MQ=255)
                                     cGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGc                                                                                                                 >  1:2043039/1‑150 (MQ=255)
                                      gAATCAACCAGATCTCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCt                                                                                                                >  1:3819842/1‑150 (MQ=255)
                                           aaCCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGtt                                                                                                           <  1:1674378/150‑1 (MQ=255)
                                            aCCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTa                                                                                                          <  2:891453/150‑1 (MQ=255)
                                               aGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGCCGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCa                                                                                                       <  2:2050106/150‑1 (MQ=255)
                                               aGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCa                                                                                                       >  1:1067670/1‑150 (MQ=255)
                                               aGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCa                                                                                                       <  1:2382426/150‑1 (MQ=255)
                                               aGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCa                                                                                                       <  1:158122/150‑1 (MQ=255)
                                                gATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCaa                                                                                                      >  2:2296200/1‑150 (MQ=255)
                                                   cgcgATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCaaaaa                                                                                                   <  2:499375/150‑1 (MQ=255)
                                                     cgATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATa                                                                                                 >  2:4268310/1‑150 (MQ=255)
                                                         aGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGa                                                                                             <  1:1845373/150‑1 (MQ=255)
                                                          gCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGat                                                                                            <  1:407904/150‑1 (MQ=255)
                                                           cGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGata                                                                                           >  1:1661013/1‑150 (MQ=255)
                                                            gACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGatat                                                                                          <  1:2062696/150‑1 (MQ=255)
                                                              cccGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTg                                                                                        <  1:2493475/150‑1 (MQ=255)
                                                              cccGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTg                                                                                        <  2:3198206/150‑1 (MQ=255)
                                                               ccGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGt                                                                                       <  2:3790256/150‑1 (MQ=255)
                                                                 ggATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCa                                                                                     >  1:4378494/1‑150 (MQ=255)
                                                                 ggATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCa                                                                                     <  1:4503625/150‑1 (MQ=255)
                                                                  gATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAg                                                                                    <  2:3823841/150‑1 (MQ=255)
                                                                  gATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAg                                                                                    >  2:590722/1‑150 (MQ=255)
                                                                   aTCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGc                                                                                   >  1:3974260/1‑150 (MQ=255)
                                                                    tCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCa                                                                                  <  2:4059435/150‑1 (MQ=255)
                                                                        ccGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGa                                                                              >  1:552957/1‑150 (MQ=255)
                                                                         cGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGAc                                                                             <  2:3952484/150‑1 (MQ=255)
                                                                          ggACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACa                                                                            >  2:2678433/1‑150 (MQ=255)
                                                                          ggACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACa                                                                            >  1:1733199/1‑150 (MQ=255)
                                                                             cTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGa                                                                         <  1:1988189/150‑1 (MQ=255)
                                                                              tGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGAt                                                                        <  1:756742/150‑1 (MQ=255)
                                                                              tGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGAt                                                                        >  2:3236867/1‑150 (MQ=255)
                                                                               gAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATc                                                                       <  1:3255886/150‑1 (MQ=255)
                                                                                  tACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAg                                                                    <  1:1492427/150‑1 (MQ=255)
                                                                                       cGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCAtct                                                               >  2:1398872/1‑150 (MQ=255)
                                                                                       cGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCAGTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCAtct                                                               >  1:2806768/1‑150 (MQ=255)
                                                                                       cGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCAtct                                                               >  2:1478901/1‑150 (MQ=255)
                                                                                        gCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATCCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCAtctc                                                              <  2:4088963/150‑1 (MQ=255)
                                                                                         cTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCg                                                             >  2:1056308/1‑150 (MQ=255)
                                                                                         cTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCg                                                             <  1:2999545/150‑1 (MQ=255)
                                                                                         cTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCg                                                             <  2:1330785/150‑1 (MQ=255)
                                                                                           gCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGAc                                                           >  2:3923012/1‑150 (MQ=255)
                                                                                                 ttGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGc                                                     >  2:3126264/1‑150 (MQ=255)
                                                                                                     gTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCg                                                 >  1:3487184/1‑150 (MQ=255)
                                                                                                     gTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCg                                                 >  1:1875877/1‑150 (MQ=255)
                                                                                                       cAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTc                                               >  2:245215/1‑150 (MQ=255)
                                                                                                          cACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCaa                                            <  1:793957/150‑1 (MQ=255)
                                                                                                           aCCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAAc                                           >  2:754099/1‑150 (MQ=255)
                                                                                                                 tGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCg                                     <  2:3610695/150‑1 (MQ=255)
                                                                                                                    aGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATa                                  >  1:2838166/1‑150 (MQ=255)
                                                                                                                    aGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATa                                  <  2:2427447/150‑1 (MQ=255)
                                                                                                                     gCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAg                                 <  1:3540156/150‑1 (MQ=255)
                                                                                                                         ctcATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAg                             >  1:3094817/1‑150 (MQ=255)
                                                                                                                          tcATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAgc                            <  1:3877314/150‑1 (MQ=255)
                                                                                                                           cATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAgcg                           <  1:4544173/150‑1 (MQ=255)
                                                                                                                            aTCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAgcgc                          >  1:4136462/1‑150 (MQ=255)
                                                                                                                              cTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCAc                        >  1:2737410/1‑150 (MQ=255)
                                                                                                                               tttCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCAcc                       >  1:3929612/1‑150 (MQ=255)
                                                                                                                                    cGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTaa                  <  1:2545753/150‑1 (MQ=255)
                                                                                                                                    cGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCGTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTaa                  <  2:1531481/150‑1 (MQ=255)
                                                                                                                                      tagtaTTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAgc                <  1:4419109/150‑1 (MQ=255)
                                                                                                                                      tagtaTTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAgc                >  2:2255185/1‑150 (MQ=255)
                                                                                                                                          aTTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCaa            >  2:2172721/1‑150 (MQ=255)
                                                                                                                                             cACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACt         >  1:3379774/1‑150 (MQ=255)
                                                                                                                                             cACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACt         >  2:3561732/1‑150 (MQ=255)
                                                                                                                                                gcgGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCa      >  1:3906523/1‑150 (MQ=255)
                                                                                                                                                 cgGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCAt     <  2:3410873/150‑1 (MQ=255)
                                                                                                                                                   gTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGt   <  1:833806/150‑1 (MQ=255)
                                                                                                                                                   gTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGt   >  1:263490/1‑150 (MQ=255)
                                                                                                                                                    tttGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGtt  <  2:393482/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTT  >  CP009273/3557667‑3557964

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: