Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,214,055 T→A 27.5% I424I (ATT→ATA ycgG → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,214,0550TA27.5% 537.5 / 251.1 316I424I (ATT→ATAycgGputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (134/95);  new base A (50/37);  total (184/132)
Fisher's exact test for biased strand distribution p-value = 8.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 326 total aligned reads displayed.

GAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATCTGAACCAAAATTA  >  CP009273/1213906‑1214204
                                                                                                                                                     |                                                                                                                                                     
gAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAAtt                                                                                                                                                       <  2:4869636/150‑1 (MQ=255)
gAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAAtt                                                                                                                                                       <  1:5257727/150‑1 (MQ=255)
 aGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTg                                                                                                                                                      >  2:4076823/1‑150 (MQ=255)
   cGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGAt                                                                                                                                                    <  2:2831417/150‑1 (MQ=255)
   cGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGAt                                                                                                                                                    <  1:825447/150‑1 (MQ=255)
   cGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGAt                                                                                                                                                    <  1:1078584/150‑1 (MQ=255)
    gTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATc                                                                                                                                                   <  2:4971095/150‑1 (MQ=255)
     tGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATca                                                                                                                                                  >  2:4607951/1‑150 (MQ=255)
     tGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATca                                                                                                                                                  <  1:205302/150‑1 (MQ=255)
       aGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATcaca                                                                                                                                                <  2:3211718/150‑1 (MQ=255)
       aGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATcaca                                                                                                                                                <  1:551950/150‑1 (MQ=255)
        gCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAg                                                                                                                                               <  2:2971173/150‑1 (MQ=255)
         ccATTGAATGTTTATGAAAGTCTGGTTCAGCGGTTTAAAATACAGCAAGAAAAATGTATTGAAATAGAAATTGATGATTTATGTACTGGCTACTCAGGGCTTTATTAACTTCATTACCTTCATATTTATTATATCAAAATTGATCACAGt                                                                                                                                              >  2:2458315/1‑150 (MQ=255)
         ccATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGt                                                                                                                                              >  2:4349079/1‑150 (MQ=255)
          cATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGtt                                                                                                                                             <  2:1431420/150‑1 (MQ=255)
          cATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGtt                                                                                                                                             >  1:958994/1‑150 (MQ=255)
           aTTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGttt                                                                                                                                            <  1:3180388/150‑1 (MQ=255)
              gAATGTTGATGAAAGTCGGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGt                                                                                                                                         <  1:1988636/150‑1 (MQ=255)
                aTGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTg                                                                                                                                       >  1:1355896/1‑150 (MQ=255)
                aTGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTg                                                                                                                                       <  1:1447936/150‑1 (MQ=255)
                       tgaAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGt                                                                                                                                >  2:1128422/1‑150 (MQ=255)
                        gaAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTa                                                                                                                               >  2:582333/1‑150 (MQ=255)
                         aaaGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTaa                                                                                                                              >  2:5194036/1‑150 (MQ=255)
                           aGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAAc                                                                                                                            <  2:1820896/150‑1 (MQ=255)
                              cTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCa                                                                                                                         >  2:4682185/1‑150 (MQ=255)
                               tGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAg                                                                                                                        <  1:489424/150‑1 (MQ=255)
                                   tCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGAccc                                                                                                                    >  1:1571559/1‑150 (MQ=255)
                                     aGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAg                                                                                                                  >  1:3363526/1‑150 (MQ=255)
                                     aGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAg                                                                                                                  >  1:4469722/1‑150 (MQ=255)
                                        ggTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAAt                                                                                                               <  1:4102200/150‑1 (MQ=255)
                                        ggTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAAt                                                                                                               >  1:4303596/1‑150 (MQ=255)
                                         gTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATc                                                                                                              >  2:3768640/1‑150 (MQ=255)
                                         gTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATc                                                                                                              >  2:74614/1‑150 (MQ=255)
                                                  aCTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAAt                                                                                                     >  1:2612572/1‑150 (MQ=255)
                                                   cTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAAtt                                                                                                    >  2:2775439/1‑150 (MQ=255)
                                                    tGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTc                                                                                                   <  1:1655338/150‑1 (MQ=255)
                                                     gCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCt                                                                                                  <  2:4854654/150‑1 (MQ=255)
                                                       aTGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTgg                                                                                                <  1:5166059/150‑1 (MQ=255)
                                                       aTGAAAATGGTTTTGTCATCGAACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTgg                                                                                                <  2:1324131/150‑1 (MQ=255)
                                                           aaaTGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATtg                                                                                            >  2:3928783/1‑150 (MQ=255)
                                                            aaTGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATtgt                                                                                           <  1:2527297/150‑1 (MQ=255)
                                                                  tttGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTg                                                                                     >  2:2309345/1‑150 (MQ=255)
                                                                       cATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTTCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGTTTGTGTATTGGATCt                                                                                <  1:5119426/150‑1 (MQ=255)
                                                                        aTCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTg                                                                               >  2:1709082/1‑150 (MQ=255)
                                                                        aTCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTg                                                                               >  2:4261842/1‑150 (MQ=255)
                                                                        aTCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTg                                                                               >  2:4021453/1‑150 (MQ=255)
                                                                         tCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTgg                                                                              <  1:3213163/150‑1 (MQ=255)
                                                                          cGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGgc                                                                             <  1:3985255/150‑1 (MQ=255)
                                                                           gCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGgcg                                                                            >  2:2027342/1‑150 (MQ=255)
                                                                               tGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTa                                                                        >  1:1795057/1‑150 (MQ=255)
                                                                                 gatgatTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTaaa                                                                      >  1:4692096/1‑150 (MQ=255)
                                                                                 gatgatTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTaaa                                                                      >  1:587351/1‑150 (MQ=255)
                                                                                 gatgagGTCGGTACGGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTaaa                                                                      <  1:1004178/150‑1 (MQ=255)
                                                                                  atgatTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAAc                                                                     >  1:2190253/1‑150 (MQ=255)
                                                                                      tttCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTc                                                                 <  1:1301701/150‑1 (MQ=255)
                                                                                          ggTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATc                                                             >  2:1170296/1‑150 (MQ=255)
                                                                                           gcaCTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCa                                                            <  2:3812581/148‑1 (MQ=255)
                                                                                                ggCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTAtc                                                       >  1:1139188/1‑150 (MQ=255)
                                                                                                ggCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTAtc                                                       >  2:3181189/1‑150 (MQ=255)
                                                                                                 gCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTAtcg                                                      <  1:3871881/150‑1 (MQ=255)
                                                                                                   tACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTAtcgtc                                                    <  1:4253883/150‑1 (MQ=255)
                                                                                                    aCTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTAGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTAtcgtcg                                                   >  1:4601114/1‑150 (MQ=255)
                                                                                                     cTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGc                                                  >  2:2743191/1‑150 (MQ=255)
                                                                                                     cTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGc                                                  >  1:3512323/1‑150 (MQ=255)
                                                                                                      tCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCt                                                 <  1:2323219/150‑1 (MQ=255)
                                                                                                         gggCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGaa                                              <  2:4003572/150‑1 (MQ=255)
                                                                                                          ggCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAg                                             <  1:3768640/150‑1 (MQ=255)
                                                                                                          ggCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAg                                             >  2:5290243/1‑150 (MQ=255)
                                                                                                           gCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAgg                                            <  1:4932867/150‑1 (MQ=255)
                                                                                                             tttCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGgtg                                          >  1:3817254/1‑150 (MQ=255)
                                                                                                               tCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTc                                        >  2:2912562/1‑150 (MQ=255)
                                                                                                               tCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTc                                        >  2:3555358/1‑150 (MQ=255)
                                                                                                                  tATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTcgaa                                     <  2:1368515/150‑1 (MQ=255)
                                                                                                                    tCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTcgaaac                                   <  2:3452135/150‑1 (MQ=255)
                                                                                                                           gACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGaa                            >  1:2685765/1‑150 (MQ=255)
                                                                                                                            aCCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGaac                           <  1:285866/150‑1 (MQ=255)
                                                                                                                             ccTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCAGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGaaca                          >  2:2145171/1‑150 (MQ=255)
                                                                                                                              cTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGaacaa                         >  2:4542998/1‑150 (MQ=255)
                                                                                                                               tGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGaacaac                        >  1:793144/1‑150 (MQ=255)
                                                                                                                                gCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACt                       >  2:3114408/1‑150 (MQ=255)
                                                                                                                                gCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACt                       <  1:666153/150‑1 (MQ=255)
                                                                                                                                 cATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACtt                      >  1:2993176/1‑150 (MQ=255)
                                                                                                                                  atatTGTTTATATCAAAATTGATCACAGTTTAGATGTCAGCGTAAACGAAGAACCAGAATCAAACCGAATTCTGGATTTTGTATTGTATCTGTCTCGTAAAATTTCTATAAGTATATTAGCTGAAGTTGTAGAAAATAAATAAAAAATTg                     >  2:690353/1‑150 (MQ=255)
                                                                                                                                  atatTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTg                     <  2:5128951/150‑1 (MQ=255)
                                                                                                                                   tatTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGa                    >  1:2042010/1‑150 (MQ=255)
                                                                                                                                   tatTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGa                    >  1:2219220/1‑150 (MQ=255)
                                                                                                                                    atTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGAc                   >  1:3628736/1‑150 (MQ=255)
                                                                                                                                    atTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGAc                   >  2:1091681/1‑150 (MQ=255)
                                                                                                                                      tGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGGGTCGAAACGAAAGAACAACTTGACTa                 >  1:4518858/1‑150 (MQ=255)
                                                                                                                                        aTTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATc               >  1:4420917/1‑150 (MQ=255)
                                                                                                                                        aTTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATc               >  1:212627/1‑150 (MQ=255)
                                                                                                                                         ttATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATCt              >  2:4166056/1‑150 (MQ=255)
                                                                                                                                         ttATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATCt              >  2:2577033/1‑150 (MQ=255)
                                                                                                                                           atatCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATCTGa            >  1:73921/1‑150 (MQ=255)
                                                                                                                                               cAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATCTGAACCa        >  2:838535/1‑150 (MQ=255)
                                                                                                                                               aaaaaTTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACACAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCTCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGATAGATCAACTTGACTATCTGAACCa        <  2:271711/149‑1 (MQ=255)
                                                                                                                                                 aaaTAGATCTCAGTTTCGTTGGCCGAGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGATGGTGTCGAAACGAAAGAACAACTTGACTATCTGAACCaaa      <  2:4809843/150‑1 (MQ=255)
                                                                                                                                                  aaTTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATCTGAACCaaaa     >  1:3939248/1‑150 (MQ=255)
                                                                                                                                                    ttGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATCTGAACCAAAAtt   <  2:1536523/150‑1 (MQ=255)
                                                                                                                                                     tGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATCTGAACCAAAATta  >  2:1619616/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATCTGAACCAAAATTA  >  CP009273/1213906‑1214204

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: