Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,189,766 C→A 7.0% T112K (ACG→AAG)  yjaG → DUF416 domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,189,7660CA 7.0% 1374.7 / 48.7 358T112K (ACG→AAG) yjaGDUF416 domain protein
Reads supporting (aligned to +/- strand):  ref base C (170/163);  new base A (15/10);  total (185/173)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 367 total aligned reads displayed.

AGACGCAAAAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATATTCC  >  CP009273/4189617‑4189915
                                                                                                                                                     |                                                                                                                                                     
aGACGCAAAAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAc                                                                                                                                                       <  2:2319025/150‑1 (MQ=255)
aGACGCAAAAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGAGACCTGGTCCATTCGCGTTTGAGTGGTGAAAc                                                                                                                                                       >  2:1019118/1‑150 (MQ=255)
   cGCAAAAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCt                                                                                                                                                    >  1:1262280/1‑150 (MQ=255)
   cGCAAAAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCt                                                                                                                                                    <  2:549464/150‑1 (MQ=255)
   cGCAAAAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCt                                                                                                                                                    <  2:3164400/150‑1 (MQ=255)
    gCAATAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTAGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTc                                                                                                                                                   <  1:3319090/150‑1 (MQ=255)
      aaacGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGATTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGa                                                                                                                                                 <  1:201240/150‑1 (MQ=255)
      aaaaGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGa                                                                                                                                                 <  1:201202/150‑1 (MQ=255)
       aaaGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGaa                                                                                                                                                >  2:2288681/1‑150 (MQ=255)
       aaaGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGaa                                                                                                                                                <  1:76026/150‑1 (MQ=255)
       aaaGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGaa                                                                                                                                                >  2:2812553/1‑150 (MQ=255)
        aaGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAac                                                                                                                                               <  1:292064/150‑1 (MQ=255)
           tAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACAcg                                                                                                                                            >  2:1142682/1‑150 (MQ=255)
               gttCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGt                                                                                                                                        <  2:1250784/149‑1 (MQ=255)
                ttCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTg                                                                                                                                       >  2:2968121/1‑150 (MQ=255)
                   gACAGCCAACTGGAGAAATTTGAAGAAGCGAGTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGaa                                                                                                                                    <  2:3997333/150‑1 (MQ=255)
                    aCAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCTACGATTTCGATCTGTACGGCGTTTATCCGGCAATCTAAGCCTGAGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAg                                                                                                                                   >  2:2084709/1‑150 (MQ=255)
                    aCAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAg                                                                                                                                   >  2:4967574/1‑150 (MQ=255)
                     cAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGt                                                                                                                                  >  2:3673326/1‑150 (MQ=255)
                     cAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGt                                                                                                                                  <  2:2280050/150‑1 (MQ=255)
                        ccAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAg                                                                                                                               >  1:672600/1‑150 (MQ=255)
                        ccAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAg                                                                                                                               >  1:1556684/1‑150 (MQ=255)
                         cAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGt                                                                                                                              <  1:3063682/150‑1 (MQ=255)
                          aaCTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTa                                                                                                                             <  2:3011416/150‑1 (MQ=255)
                           aCTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTaa                                                                                                                            >  2:5126933/1‑150 (MQ=255)
                            cTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAg                                                                                                                           >  1:4087968/1‑150 (MQ=255)
                             tGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGa                                                                                                                          >  1:3988276/1‑150 (MQ=255)
                             tGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGa                                                                                                                          <  2:1782444/150‑1 (MQ=255)
                              ggAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGAc                                                                                                                         <  2:4556344/150‑1 (MQ=255)
                              ggAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGAc                                                                                                                         <  2:2002901/150‑1 (MQ=255)
                                agaAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCt                                                                                                                       >  1:3717092/1‑150 (MQ=255)
                                   aaTTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCa                                                                                                                    >  1:4205485/1‑150 (MQ=255)
                                        gaagaaGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGGTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACg                                                                                                               >  2:3754825/1‑150 (MQ=255)
                                          agaagCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGAc                                                                                                             <  1:2387763/150‑1 (MQ=255)
                                           gaagCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGAcc                                                                                                            >  2:1059008/1‑150 (MQ=255)
                                            aagCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCg                                                                                                           >  1:2683316/1‑150 (MQ=255)
                                               cGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGTCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTg                                                                                                        <  1:4544823/150‑1 (MQ=255)
                                               cGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTg                                                                                                        >  2:4063007/1‑150 (MQ=255)
                                                gATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGc                                                                                                       >  2:1373829/1‑150 (MQ=255)
                                                gATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGc                                                                                                       <  2:2893880/150‑1 (MQ=255)
                                                   tCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGAt                                                                                                    >  2:3213133/1‑150 (MQ=255)
                                                     cTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGc                                                                                                  >  1:4704619/1‑150 (MQ=255)
                                                      ttCAGCCGACGATTTCGATCTGTTCGTCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCt                                                                                                 <  2:4231105/150‑1 (MQ=255)
                                                      ttCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCt                                                                                                 >  2:2506691/1‑150 (MQ=255)
                                                            cgacgaTTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAAt                                                                                           <  1:2278427/150‑1 (MQ=255)
                                                               cgaTTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGAc                                                                                        >  2:1222184/1‑150 (MQ=255)
                                                                gaTTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACAGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACt                                                                                       >  2:3703627/1‑150 (MQ=255)
                                                                   ttCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAg                                                                                    >  1:1405863/1‑150 (MQ=255)
                                                                   ttCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAg                                                                                    <  1:3005681/150‑1 (MQ=255)
                                                                   ttCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAg                                                                                    <  1:632573/150‑1 (MQ=255)
                                                                     cggtCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGGTGGGGGGGTGAAACGCTCGCACCCGCGGTGGCAGTGAGTAAGACCGCCATCACGACCGGTGCGATGCTGGACCTGACTCAGGc                                                                                  >  1:4001098/4‑150 (MQ=255)
                                                                      gATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCt                                                                                 >  1:459637/1‑150 (MQ=255)
                                                                       aTCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTg                                                                                >  2:3479646/1‑150 (MQ=255)
                                                                           gTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTcg                                                                            >  1:3210435/1‑150 (MQ=255)
                                                                            tACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTcgc                                                                           >  2:3346194/1‑150 (MQ=255)
                                                                             aCGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTcgcg                                                                          >  1:3454712/1‑150 (MQ=255)
                                                                                 cGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAAt                                                                      >  1:5082959/1‑150 (MQ=255)
                                                                                  gTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAAta                                                                     <  2:2684567/150‑2 (MQ=255)
                                                                                  gTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATg                                                                     <  2:4278317/150‑1 (MQ=255)
                                                                                       tCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGa                                                                <  2:4770361/150‑1 (MQ=255)
                                                                                         cGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATg                                                              >  2:440845/1‑150 (MQ=255)
                                                                                         cGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATg                                                              >  2:1692847/1‑150 (MQ=255)
                                                                                          ggCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGa                                                             <  2:1336701/150‑1 (MQ=255)
                                                                                             aaTCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAaga                                                          >  1:579675/1‑150 (MQ=255)
                                                                                             aaTCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAaga                                                          >  1:2148938/1‑150 (MQ=255)
                                                                                             aaTCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAaga                                                          >  2:1823562/1‑150 (MQ=255)
                                                                                             aaTCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAaga                                                          >  1:3591836/1‑150 (MQ=255)
                                                                                             aaTCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGATCGAACAAGCGTTGGAAGTGATTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAaga                                                          >  1:2811015/1‑150 (MQ=255)
                                                                                                cGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCt                                                       >  1:1695440/1‑150 (MQ=255)
                                                                                                  aTGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTa                                                     >  2:2192742/1‑150 (MQ=255)
                                                                                                  aTGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTa                                                     <  2:4349253/150‑1 (MQ=255)
                                                                                                         cGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGaaaa                                              <  1:4159818/150‑1 (MQ=255)
                                                                                                          gTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAAc                                             >  2:1176618/1‑150 (MQ=255)
                                                                                                          gTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAAc                                             >  1:4022943/1‑150 (MQ=255)
                                                                                                           tGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAAcc                                            >  2:167636/1‑150 (MQ=255)
                                                                                                             gCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCa                                          >  1:3350861/1‑150 (MQ=255)
                                                                                                               gTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGc                                        <  2:3286430/150‑1 (MQ=255)
                                                                                                                 tAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTg                                      <  2:3389878/150‑1 (MQ=255)
                                                                                                                 tAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTg                                      <  2:3436411/150‑1 (MQ=255)
                                                                                                                 tAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTg                                      >  2:3922992/1‑150 (MQ=255)
                                                                                                                 tAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTg                                      >  2:624690/1‑150 (MQ=255)
                                                                                                                  aaGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGt                                     <  2:160431/150‑1 (MQ=255)
                                                                                                                       aaCTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGc                                >  2:2093272/1‑150 (MQ=255)
                                                                                                                        aCTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCa                               >  2:2439308/1‑150 (MQ=255)
                                                                                                                           ggTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGa                            >  1:1005333/1‑150 (MQ=255)
                                                                                                                           ggTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGGAGAGCAAGa                            <  2:3824542/150‑1 (MQ=255)
                                                                                                                            gTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGaa                           >  1:3929542/1‑150 (MQ=255)
                                                                                                                             tCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAAt                          <  2:4638275/150‑1 (MQ=255)
                                                                                                                              ccATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATg                         >  1:4776365/1‑150 (MQ=255)
                                                                                                                                aTTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATggg                       <  1:4466704/150‑1 (MQ=255)
                                                                                                                                     cgTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACAtt                  <  1:2506691/150‑1 (MQ=255)
                                                                                                                                     cgTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACAtt                  <  1:3253310/150‑1 (MQ=255)
                                                                                                                                            gtggtgAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTggg           >  1:805236/1‑150 (MQ=255)
                                                                                                                                                tgAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAat       >  2:1454455/1‑150 (MQ=255)
                                                                                                                                                 gAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGACata      >  2:787619/1‑150 (MQ=255)
                                                                                                                                                 gAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAata      >  1:4826873/1‑150 (MQ=255)
                                                                                                                                                   aaCGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATAtt    >  2:4718151/1‑150 (MQ=255)
                                                                                                                                                   aaCGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATAtt    <  1:4246502/150‑1 (MQ=255)
                                                                                                                                                   aaCGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATAtt    >  1:895978/1‑150 (MQ=255)
                                                                                                                                                     cGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATATTcc  >  1:1746106/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGACGCAAAAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATATTCC  >  CP009273/4189617‑4189915

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: