Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,214,055 T→A 18.4% 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,0550TA18.4% 580.9 / 120.1 234I424I (ATT→ATAycgGputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (107/84);  new base A (22/21);  total (129/105)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

Only 100 of 237 total aligned reads displayed.

CGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATCTGAACCAAAA  >  CP009273/1213909‑1214201
                                                                                                                                                  |                                                                                                                                                  
cGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTg                                                                                                                                                   <  2:217691/148‑1 (MQ=255)
 gTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGa                                                                                                                                                  <  1:406365/148‑1 (MQ=255)
  tGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGAt                                                                                                                                                 >  2:3669387/1‑148 (MQ=255)
     gCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATcac                                                                                                                                              <  1:1692343/148‑1 (MQ=255)
     gCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATcac                                                                                                                                              <  1:2848329/148‑1 (MQ=255)
      ccATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATcaca                                                                                                                                             >  2:1825138/1‑148 (MQ=255)
       cATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAg                                                                                                                                            <  2:1800764/148‑1 (MQ=255)
       cATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAg                                                                                                                                            <  2:1800708/148‑1 (MQ=255)
       cATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAg                                                                                                                                            <  1:1874512/148‑1 (MQ=255)
        aTTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGt                                                                                                                                           <  1:3766105/148‑1 (MQ=255)
         ttGAATGTTGATGAAAGTCTGGTTCAGCTGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGtt                                                                                                                                          >  2:1448301/1‑148 (MQ=255)
         ttGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGtt                                                                                                                                          >  1:560573/1‑148 (MQ=255)
         ttGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGtt                                                                                                                                          >  2:3429741/1‑148 (MQ=255)
         ttGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGtt                                                                                                                                          <  2:1567033/148‑1 (MQ=255)
         ttGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGtt                                                                                                                                          >  1:978388/1‑148 (MQ=255)
          tGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATTATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGttt                                                                                                                                         <  2:1163876/148‑1 (MQ=255)
          tGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGttt                                                                                                                                         >  1:3605195/1‑148 (MQ=255)
          tGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGttt                                                                                                                                         >  2:3151476/1‑148 (MQ=255)
          tGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGttt                                                                                                                                         >  1:531871/1‑148 (MQ=255)
          tGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGttt                                                                                                                                         <  2:31347/148‑1 (MQ=255)
              tGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGtt                                                                                                                                     >  1:3371762/1‑148 (MQ=255)
              tGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGtt                                                                                                                                     >  1:2436386/1‑148 (MQ=255)
                ttGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTgg                                                                                                                                   >  1:2025833/1‑148 (MQ=255)
                 tgatgaAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTgtc                                                                                                                                  >  1:2002116/1‑146 (MQ=255)
                    tgaAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCcgc                                                                                                                               >  1:3520182/1‑148 (MQ=255)
                         gTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTaaa                                                                                                                          >  1:105028/1‑148 (MQ=255)
                         gTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTaaa                                                                                                                          <  2:2093271/148‑1 (MQ=255)
                             ggTTCAGCGGTTGAACATGCTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCa                                                                                                                      >  1:1604806/1‑148 (MQ=255)
                                tCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGAc                                                                                                                   <  1:2115309/148‑1 (MQ=255)
                                 cAGCGGTTGAACCGACTGCCTGAAAATGGTATTGTCATCGCACTGGATGATTGCGGTACTGGCTACTCAGGTCTATCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGAcc                                                                                                                  <  2:1990683/148‑1 (MQ=255)
                                   gCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCa                                                                                                                <  1:249406/148‑1 (MQ=255)
                                    cGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAg                                                                                                               >  1:3720181/1‑148 (MQ=255)
                                      gTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGaa                                                                                                             <  1:1649178/148‑1 (MQ=255)
                                        tGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATc                                                                                                           <  1:759738/148‑1 (MQ=255)
                                        tGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATc                                                                                                           <  2:1221804/148‑1 (MQ=255)
                                        tGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATc                                                                                                           <  2:2318514/148‑1 (MQ=255)
                                         gAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCa                                                                                                          <  2:2984867/148‑1 (MQ=255)
                                          aaCATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCaa                                                                                                         >  1:2982427/1‑148 (MQ=255)
                                          aaCATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCaa                                                                                                         <  1:785049/148‑1 (MQ=255)
                                           aCATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAAc                                                                                                        <  1:1371572/148‑1 (MQ=255)
                                           aCATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAAc                                                                                                        <  1:3274723/148‑1 (MQ=255)
                                             aTACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAAccc                                                                                                      >  2:2701599/1‑148 (MQ=255)
                                              tACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCg                                                                                                     >  1:1433178/1‑148 (MQ=255)
                                               aCTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGa                                                                                                    <  2:810437/148‑1 (MQ=255)
                                               aCTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGa                                                                                                    >  2:2721346/1‑148 (MQ=255)
                                                cTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGaa                                                                                                   <  2:2869705/148‑1 (MQ=255)
                                                   cATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTc                                                                                                >  1:3295812/1‑148 (MQ=255)
                                                   cATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTc                                                                                                >  2:3047730/1‑148 (MQ=255)
                                                   cATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTc                                                                                                <  2:728247/148‑1 (MQ=255)
                                                     tGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTg                                                                                              >  1:647813/1‑148 (MQ=255)
                                                       aaaaTGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGa                                                                                            >  1:2859307/1‑148 (MQ=255)
                                                           tGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATtgt                                                                                        >  2:662480/1‑148 (MQ=255)
                                                           tGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACATCAATATTGATAATATCAAAATTGATCACAGTTTCGATGTCAGCGTAAACGCAGACCCAGAATCAACACGAAATCTGGATagt                                                                                        >  2:3013872/1‑148 (MQ=255)
                                                            ggTTTAGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATtgtg                                                                                       >  2:1616222/1‑148 (MQ=255)
                                                              ttttGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTa                                                                                     <  2:562901/148‑1 (MQ=255)
                                                                 tGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTg                                                                                  >  2:460579/1‑148 (MQ=255)
                                                                   tCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGa                                                                                >  1:1283065/1‑148 (MQ=255)
                                                                   tCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGa                                                                                >  1:1380323/1‑148 (MQ=255)
                                                                    cATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGAt                                                                               <  1:3219933/148‑1 (MQ=255)
                                                                       cGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTg                                                                            >  1:274698/1‑148 (MQ=255)
                                                                           cTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGgcgc                                                                        >  2:2383978/1‑148 (MQ=255)
                                                                           cTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGgcgc                                                                        >  2:3698735/1‑148 (MQ=255)
                                                                            tGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCAGGATTGTGTATTGGATCTGGcgcg                                                                       >  1:3456885/1‑148 (MQ=255)
                                                                               atgatTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTaa                                                                    >  1:1823834/1‑148 (MQ=255)
                                                                                tgatTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTaaa                                                                   <  1:3669387/148‑1 (MQ=255)
                                                                                       ggTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGa                                                            >  1:881591/1‑148 (MQ=255)
                                                                                         tACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATc                                                          <  2:2876348/148‑1 (MQ=255)
                                                                                            tGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGt                                                       <  2:17110/148‑1 (MQ=255)
                                                                                             ggCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTa                                                      <  2:3775574/148‑1 (MQ=255)
                                                                                                    cAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGc                                               >  1:2147995/1‑148 (MQ=255)
                                                                                                      gggCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTg                                             >  2:1328469/1‑148 (MQ=255)
                                                                                                       ggCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGa                                            <  2:161316/148‑1 (MQ=255)
                                                                                                          tttCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAgg                                         >  1:2580917/1‑148 (MQ=255)
                                                                                                          tttCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAgg                                         >  2:3556128/1‑148 (MQ=255)
                                                                                                            tCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGgtg                                       <  1:3770934/148‑1 (MQ=255)
                                                                                                            tCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGgtg                                       >  1:36162/1‑148 (MQ=255)
                                                                                                             cTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGgtgt                                      <  2:1683663/148‑1 (MQ=255)
                                                                                                                 tCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGAGTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTcgaa                                  >  2:3336560/1‑148 (MQ=255)
                                                                                                                    tCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTcgaaacg                               >  2:2900421/1‑148 (MQ=255)
                                                                                                                    tCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTcgaaacg                               >  2:3184497/1‑148 (MQ=255)
                                                                                                                    tCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTcgaaacg                               >  1:3422655/1‑148 (MQ=255)
                                                                                                                    tCATGACCTGCATATTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTcgaaacg                               >  2:3361441/1‑148 (MQ=255)
                                                                                                                     cATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTcgaaacga                              <  2:2025833/148‑1 (MQ=255)
                                                                                                                       tGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTcgaaacgaaa                            <  1:2290989/148‑1 (MQ=255)
                                                                                                                          ccTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGaa                         >  2:3742374/1‑148 (MQ=255)
                                                                                                                            tGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGaaca                       >  2:2583319/1‑148 (MQ=255)
                                                                                                                            tGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGaaca                       <  1:3151476/148‑1 (MQ=255)
                                                                                                                              cATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGaacaac                     >  1:127325/1‑148 (MQ=255)
                                                                                                                              cATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGaacaac                     >  2:614185/1‑148 (MQ=255)
                                                                                                                                 atTGATTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTg                  >  1:1941069/1‑148 (MQ=255)
                                                                                                                                  ttGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGa                 <  2:3385485/148‑1 (MQ=255)
                                                                                                                                     aTTATATCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTa              <  1:2721346/148‑1 (MQ=255)
                                                                                                                                      ttATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTAt             >  2:3596744/1‑148 (MQ=255)
                                                                                                                                       tataTCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATc            >  2:105976/1‑148 (MQ=255)
                                                                                                                                       tataTCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATc            >  1:2089317/1‑148 (MQ=255)
                                                                                                                                       tataTCAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATc            >  2:854588/1‑148 (MQ=255)
                                                                                                                                          atCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATCTGa         >  1:736339/1‑148 (MQ=255)
                                                                                                                                            cAAAATAGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATCTGAAc       >  1:2059244/1‑148 (MQ=255)
                                                                                                                                                aTTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATCTGAACCaaa   <  2:531871/148‑1 (MQ=255)
                                                                                                                                                aTTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATCTGAACCaaa   <  1:3271960/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACATACTGCATGAAAATGGTTTTGTCATCGCACTGGATGATTTCGGTACTGGCTACTCAGGGCTTTCTTATCTTCATGACCTGCATATTGATTATATCAAAATTGATCACAGTTTCGTTGGCCGCGTAAACGCAGACCCAGAATCAACCCGAATTCTGGATTGTGTATTGGATCTGGCGCGTAAACTTTCGATCAGTATCGTCGCTGAAGGTGTCGAAACGAAAGAACAACTTGACTATCTGAACCAAAA  >  CP009273/1213909‑1214201

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: