Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 279,334 C→A 90.2% I141I (ATC→ATA)  yagF → CP4‑6 prophage; dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273279,3340CA90.2% 773.7 / 79.0 294I141I (ATC→ATA) yagFCP4‑6 prophage; dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base C (14/14);  new base A (123/142);  total (137/157)
Fisher's exact test for biased strand distribution p-value = 8.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

Only 100 of 294 total aligned reads displayed.

GCATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGACGC  >  CP009273/279189‑279479
                                                                                                                                                 |                                                                                                                                                 
 cATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAgtgt                                                                                                                                               <  1:3057885/149‑1 (MQ=255)
    gcagatgcagGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTc                                                                                                                                            <  2:3168311/149‑1 (MQ=255)
          gcagGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCgccgcc                                                                                                                                      <  2:79817/149‑1 (MQ=255)
            agGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTg                                                                                                                                    >  2:1177431/1‑149 (MQ=255)
            agGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTg                                                                                                                                    >  2:3805823/1‑149 (MQ=255)
             ggcggcggcGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGa                                                                                                                                   >  2:4683211/1‑149 (MQ=255)
              gcggcggcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGAt                                                                                                                                  <  1:4161413/149‑1 (MQ=255)
               cggcggcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATc                                                                                                                                 <  1:586451/149‑1 (MQ=255)
               cggcggcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATc                                                                                                                                 <  2:3499394/149‑1 (MQ=255)
                  cggcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGc                                                                                                                              <  1:4036014/149‑1 (MQ=255)
                   ggcgAAGGAGATCACCCGCACTGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCt                                                                                                                             <  2:3221346/149‑1 (MQ=255)
                   ggcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCACCGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCt                                                                                                                             >  2:2477255/1‑149 (MQ=255)
                   ggcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCt                                                                                                                             <  2:1960765/149‑1 (MQ=255)
                      gAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTccc                                                                                                                          >  1:1951078/1‑149 (MQ=255)
                       aaGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCt                                                                                                                         >  2:3399578/1‑149 (MQ=255)
                       aaGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCt                                                                                                                         >  2:2030743/1‑149 (MQ=255)
                       aaGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCt                                                                                                                         <  2:2104245/149‑1 (MQ=255)
                            gaTCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGa                                                                                                                    >  2:1270042/1‑149 (MQ=255)
                               cACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGAcgc                                                                                                                 <  1:4461835/149‑1 (MQ=255)
                                  ccGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACgcggc                                                                                                              >  1:1637725/1‑149 (MQ=255)
                                  ccGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACgcggc                                                                                                              >  1:4025613/1‑149 (MQ=255)
                                    gCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGcggcgg                                                                                                            >  2:2672426/1‑149 (MQ=255)
                                    gCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGcggcgg                                                                                                            >  1:395107/1‑149 (MQ=255)
                                    gCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGcggcgg                                                                                                            <  1:959988/149‑1 (MQ=255)
                                          gCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTg                                                                                                      <  2:514879/149‑1 (MQ=255)
                                            gggATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGAt                                                                                                    >  1:4304580/1‑149 (MQ=255)
                                                tCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTCTGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGc                                                                                                >  1:1337669/1‑149 (MQ=255)
                                                  ccGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGt                                                                                              <  2:4291141/149‑1 (MQ=255)
                                                  ccGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGt                                                                                              <  2:2577429/149‑1 (MQ=255)
                                                  ccGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGt                                                                                              <  1:255715/149‑1 (MQ=255)
                                                   cGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTa                                                                                             <  1:3637406/149‑1 (MQ=255)
                                                     ttCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGc                                                                                           <  1:4666206/149‑1 (MQ=255)
                                                        gcgGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGAc                                                                                        <  1:914325/149‑1 (MQ=255)
                                                        gcgGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGAc                                                                                        <  2:2287731/149‑1 (MQ=255)
                                                        gcgGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGAc                                                                                        >  1:3766905/1‑149 (MQ=255)
                                                         cgGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGAcc                                                                                       <  1:254762/149‑1 (MQ=255)
                                                            ccTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGc                                                                                    >  2:3494456/1‑149 (MQ=255)
                                                             cTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCg                                                                                   >  2:1323235/1‑149 (MQ=255)
                                                              ttCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGa                                                                                  <  1:525802/149‑1 (MQ=255)
                                                              ttCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGa                                                                                  <  2:2000915/149‑1 (MQ=255)
                                                              ttCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGa                                                                                  >  1:4516091/1‑149 (MQ=255)
                                                                 gTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATaa                                                                               <  1:1928725/149‑1 (MQ=255)
                                                                   cAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAg                                                                             <  1:3217902/149‑1 (MQ=255)
                                                                    aGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAgg                                                                            <  1:720725/149‑1 (MQ=255)
                                                                      cGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGc                                                                          <  1:395864/149‑1 (MQ=255)
                                                                      cGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGc                                                                          <  2:694655/149‑1 (MQ=255)
                                                                        aTCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTg                                                                        >  1:3697406/1‑149 (MQ=255)
                                                                         tCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGc                                                                       <  1:1405880/149‑1 (MQ=255)
                                                                          ccGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGcc                                                                      <  1:2185482/149‑1 (MQ=255)
                                                                          ccGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGcc                                                                      <  1:2419718/149‑1 (MQ=255)
                                                                           cGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGccc                                                                     >  1:109222/1‑149 (MQ=255)
                                                                           cGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGccc                                                                     >  1:622605/1‑149 (MQ=255)
                                                                              gCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGcc                                                                  <  1:3283517/149‑1 (MQ=255)
                                                                                gACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGccac                                                                >  2:2038225/1‑149 (MQ=255)
                                                                                    ggCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATg                                                            >  2:1828352/1‑149 (MQ=255)
                                                                                     gcgcTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGa                                                           >  2:3383163/1‑149 (MQ=255)
                                                                                      cgcTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGAt                                                          <  1:3080162/149‑1 (MQ=255)
                                                                                       gcTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGAtt                                                         <  1:2827609/149‑1 (MQ=255)
                                                                                       gcTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGAtt                                                         >  1:2936683/1‑149 (MQ=255)
                                                                                       gcTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGAtt                                                         <  1:434217/149‑1 (MQ=255)
                                                                                          cGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTgcg                                                      >  2:3099646/1‑149 (MQ=255)
                                                                                           gCAGGGCACGCACGGTATGTTCGAGTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGGTCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTgcgc                                                     <  1:842767/149‑1 (MQ=255)
                                                                                            cAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCt                                                    <  1:3195243/149‑1 (MQ=255)
                                                                                             aGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTg                                                   >  1:2310485/1‑149 (MQ=255)
                                                                                             aGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTg                                                   <  2:3932889/149‑1 (MQ=255)
                                                                                             aGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGAGCAAGGGCTGCCCGCCACCATGATTGCGCTg                                                   >  2:1810063/1‑149 (MQ=255)
                                                                                              gggCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTCGCGACCTGCGAGAAAGGGCTGCCCGCCACCATGATTGCGCTgg                                                  >  1:2068822/1‑149 (MQ=255)
                                                                                              gggCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTgg                                                  >  1:2709036/1‑149 (MQ=255)
                                                                                                  acgcacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCcgc                                              <  2:2421531/149‑1 (MQ=255)
                                                                                                  acgcacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCcgc                                              >  2:2605140/1‑149 (MQ=255)
                                                                                                  acgcacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCcgc                                              >  1:2915314/1‑149 (MQ=255)
                                                                                                  acgcacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCcgc                                              <  2:2964429/149‑1 (MQ=255)
                                                                                                  acgcacgGTATGTTCGAGTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCcgc                                              <  2:2872155/149‑1 (MQ=255)
                                                                                                    gcacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGa                                            >  1:170871/1‑149 (MQ=255)
                                                                                                     cacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGAt                                           <  2:2745383/149‑1 (MQ=255)
                                                                                                      acgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATg                                          >  1:2857973/1‑149 (MQ=255)
                                                                                                        ggTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCa                                        <  1:1256624/149‑1 (MQ=255)
                                                                                                          tATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACg                                      >  2:1640177/1‑149 (MQ=255)
                                                                                                           aTGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGa                                     >  2:1909779/1‑149 (MQ=255)
                                                                                                            tGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGAc                                    >  1:1182491/1‑149 (MQ=255)
                                                                                                             gTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGAcc                                   <  1:13211/149‑1 (MQ=255)
                                                                                                                cGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGc                                >  2:188643/1‑149 (MQ=255)
                                                                                                                 gATTCCCTGCCTTACCGCAACGAGGGGGCGTTGCTGTGTCTCCTCCTGATCGGCTCCCAGACGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGcc                               <  1:3565987/149‑1 (MQ=255)
                                                                                                                  aTTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCAGGATTGCGCTGGCCGCGATGCACGACCTGCCg                              >  2:4585012/1‑149 (MQ=255)
                                                                                                                     cccTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACt                           <  1:3848606/149‑1 (MQ=255)
                                                                                                                      ccTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTa                          <  2:2408814/149‑1 (MQ=255)
                                                                                                                       cTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTAt                         <  1:1034186/149‑1 (MQ=255)
                                                                                                                        tGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTAtt                        >  1:3025997/1‑149 (MQ=255)
                                                                                                                             tACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGt                   >  2:1471816/1‑149 (MQ=255)
                                                                                                                             tACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGt                   >  2:1083025/1‑149 (MQ=255)
                                                                                                                             tACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGt                   >  2:4178201/1‑149 (MQ=255)
                                                                                                                                 gCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcg               <  1:3453168/149‑1 (MQ=255)
                                                                                                                                 gCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcg               >  1:1674391/1‑149 (MQ=255)
                                                                                                                                   aaCGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcggg             <  2:595855/149‑1 (MQ=255)
                                                                                                                                    acgacgCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATCAAGGGCTGCCCGCCACCCTGATTGCGCTGGCCGCGATGCACGACCTGCCGACTAGTCTGGTGCcgggc            >  2:2217269/1‑149 (MQ=255)
                                                                                                                                    acgacgCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcgggc            >  2:856421/1‑149 (MQ=255)
                                                                                                                                    acgacgCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcgggc            <  1:266338/149‑1 (MQ=255)
                                                                                                                                       acgCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcgggcggg         <  1:3877681/149‑1 (MQ=255)
                                                                                                                                        cgcgGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCgggcgggg        <  1:4703423/149‑1 (MQ=255)
                                                                                                                                         gcggcgATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCgggcggggc       >  1:1555721/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GCATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGACGC  >  CP009273/279189‑279479

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: ‑