Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,545,608 C→A 100% R177I (AGA→ATA)  yjiC ← uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,545,6080CA97.7% 702.9 / 4.6 221R177I (AGA→ATA) yjiCuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base A (117/99);  total (117/104)
Fisher's exact test for biased strand distribution p-value = 2.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 221 total aligned reads displayed.

GCTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAATTTAGC  >  CP009273/4545461‑4545754
                                                                                                                                                   |                                                                                                                                                  
gCTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATATATATCAATCCTTGGCTGCTCATTGGTCCATAtt                                                                                                                                                  >  1:3452788/1‑150 (MQ=255)
gCTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATAtt                                                                                                                                                  >  2:1588081/1‑150 (MQ=255)
   tCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATAttttt                                                                                                                                               <  1:513909/150‑1 (MQ=255)
    cTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCTATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTg                                                                                                                                              <  2:556864/150‑1 (MQ=255)
    cTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTg                                                                                                                                              >  2:1233197/1‑150 (MQ=255)
     tataCCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGACTCCGGCAGAATGTGCTACCCCTGTTCCCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTga                                                                                                                                             >  2:5079485/1‑149 (MQ=255)
     tataCCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGATGATCATTGGTCCATATTTTTGc                                                                                                                                             >  2:135534/1‑150 (MQ=255)
     tataACAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCATTGGCTGCTCATTGGTCCATATTTTTGc                                                                                                                                             >  2:4577823/1‑150 (MQ=255)
      ataCCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGcc                                                                                                                                            <  1:4784038/150‑1 (MQ=255)
        aCCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTg                                                                                                                                          <  2:3083886/150‑1 (MQ=255)
         ccAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTga                                                                                                                                         <  2:4352442/150‑1 (MQ=255)
           aaTTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTgaga                                                                                                                                       >  1:2420441/1‑150 (MQ=255)
            attattGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATNCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGaa                                                                                                                                      <  1:2680508/150‑1 (MQ=255)
            attattGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGaa                                                                                                                                      <  1:756994/150‑1 (MQ=255)
            attattGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGaa                                                                                                                                      <  2:2696321/150‑1 (MQ=255)
               attGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCt                                                                                                                                   <  2:1175649/150‑1 (MQ=255)
                 tGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTgg                                                                                                                                 <  1:12453/150‑1 (MQ=255)
                  gggTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGa                                                                                                                                >  1:4871795/1‑150 (MQ=255)
                   ggTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGaa                                                                                                                               <  2:754289/150‑1 (MQ=255)
                   ggTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGaa                                                                                                                               >  2:2487717/1‑150 (MQ=255)
                    gTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAg                                                                                                                              >  1:1001025/1‑150 (MQ=255)
                     tCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAgg                                                                                                                             <  1:4245686/150‑1 (MQ=255)
                      cATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAggg                                                                                                                            >  2:2710769/1‑150 (MQ=255)
                      cATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAggg                                                                                                                            >  2:1290642/1‑150 (MQ=255)
                      cATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAggg                                                                                                                            >  1:4731742/1‑150 (MQ=255)
                        ttACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGtt                                                                                                                          >  1:439739/1‑150 (MQ=255)
                          aCCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGa                                                                                                                        >  1:3640467/1‑150 (MQ=255)
                           cccAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGaa                                                                                                                       <  1:3018295/150‑1 (MQ=255)
                           cccAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATGTTTGCCTGAGAAGCTGGAAGGGTTGaa                                                                                                                       >  2:2638007/1‑150 (MQ=255)
                            ccAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGaaa                                                                                                                      >  2:2577750/1‑150 (MQ=255)
                               aCTCAATATTAGCATTAAGTAATGCGGAAATGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATAAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAATGGTTGAAatat                                                                                                                   >  1:150895/1‑150 (MQ=255)
                                cTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAatata                                                                                                                  >  1:3593994/1‑150 (MQ=255)
                                cTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAatata                                                                                                                  <  1:1657876/150‑1 (MQ=255)
                                cTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAatata                                                                                                                  >  1:4163292/1‑150 (MQ=255)
                                 tCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATGTTCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAatatat                                                                                                                 <  2:586220/150‑1 (MQ=255)
                                  cAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATc                                                                                                                >  1:337538/1‑150 (MQ=255)
                                   aaTATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATcc                                                                                                               >  2:4527596/1‑150 (MQ=255)
                                       ttAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATc                                                                                                           <  2:28857/150‑1 (MQ=255)
                                       ttAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATc                                                                                                           <  2:792370/150‑1 (MQ=255)
                                          gCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTg                                                                                                        <  1:4931898/150‑1 (MQ=255)
                                             ttAAGTAATGCGGCAAGGGCGACTCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGtt                                                                                                     >  1:3808408/1‑150 (MQ=255)
                                                aGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCg                                                                                                  >  1:1456306/1‑150 (MQ=255)
                                                 gTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCgg                                                                                                 >  2:33439/1‑150 (MQ=255)
                                                  tAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGc                                                                                                >  1:1687094/1‑150 (MQ=255)
                                                  tAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGc                                                                                                >  1:1190027/1‑150 (MQ=255)
                                                   aaTGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCa                                                                                               <  1:3997736/150‑1 (MQ=255)
                                                      gCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATaa                                                                                            >  1:3915379/1‑150 (MQ=255)
                                                       cGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAAc                                                                                           >  2:160482/1‑150 (MQ=255)
                                                       cGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAAc                                                                                           <  2:626419/150‑1 (MQ=255)
                                                          cAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAg                                                                                        >  2:87386/1‑150 (MQ=255)
                                                           aaGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAgg                                                                                       <  1:96885/150‑1 (MQ=255)
                                                                cGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACAt                                                                                  >  2:2209755/1‑150 (MQ=255)
                                                                 taaTCCGGAAGAATGTGCTACACCTGTTACCGTCATTATCTGATAAATATAGATATCACTCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACAtt                                                                                 <  1:2870764/149‑1 (MQ=255)
                                                                    tCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCa                                                                              >  1:682092/1‑150 (MQ=255)
                                                                     ccGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCaa                                                                             <  2:3789633/150‑1 (MQ=255)
                                                                     ccGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCaa                                                                             <  2:3789932/150‑1 (MQ=255)
                                                                      cGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCCGGAACATTGCAAc                                                                            >  2:3551435/1‑150 (MQ=255)
                                                                          agaaTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAgg                                                                        <  2:3575776/150‑1 (MQ=255)
                                                                          agaaTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAgg                                                                        >  2:2912069/1‑150 (MQ=255)
                                                                           gaaTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAggg                                                                       <  2:3529214/150‑1 (MQ=255)
                                                                           gaaTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAggg                                                                       <  2:4230577/150‑1 (MQ=255)
                                                                              tttGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCt                                                                    <  1:2945614/150‑1 (MQ=255)
                                                                                 gCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAAt                                                                 >  1:4608156/1‑150 (MQ=255)
                                                                                 gCTACCCCTGTGACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAAt                                                                 <  2:3479434/150‑1 (MQ=255)
                                                                                 gCTACCACTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAAt                                                                 <  2:2773540/150‑1 (MQ=255)
                                                                                  cTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATg                                                                >  1:827584/1‑150 (MQ=255)
                                                                                  cTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATg                                                                >  2:4151315/1‑150 (MQ=255)
                                                                                   tACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGt                                                               >  1:2664917/1‑150 (MQ=255)
                                                                                                 tCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTg                                                 <  1:2963819/150‑1 (MQ=255)
                                                                                                        cTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTc                                          <  1:3509779/150‑1 (MQ=255)
                                                                                                             cAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCaa                                     >  1:1800011/1‑150 (MQ=255)
                                                                                                             cAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCaa                                     <  1:1991317/150‑1 (MQ=255)
                                                                                                              aaTATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAAc                                    >  2:4750671/1‑150 (MQ=255)
                                                                                                                 ataGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAgg                                 >  1:2580007/1‑150 (MQ=255)
                                                                                                                  taGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGa                                >  2:2169691/1‑150 (MQ=255)
                                                                                                                    gATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATc                              <  2:1394286/150‑1 (MQ=255)
                                                                                                                     atatCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCg                             <  2:4398805/150‑1 (MQ=255)
                                                                                                                     atatCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCg                             <  1:1072058/150‑1 (MQ=255)
                                                                                                                        tCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGaaa                          >  1:4483288/1‑150 (MQ=255)
                                                                                                                        tCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGaaa                          <  2:2284242/150‑1 (MQ=255)
                                                                                                                        tCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGaaa                          <  1:1217094/150‑1 (MQ=255)
                                                                                                                         cAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAg                         >  2:1691575/1‑150 (MQ=255)
                                                                                                                          aaTCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAgg                        >  2:4478055/1‑150 (MQ=255)
                                                                                                                          aaTCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAgg                        >  1:925345/1‑150 (MQ=255)
                                                                                                                               ttGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGctct                   >  2:2394483/1‑150 (MQ=255)
                                                                                                                               ttGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGctct                   >  2:3895470/1‑150 (MQ=255)
                                                                                                                               ttGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGctct                   >  2:1991700/1‑150 (MQ=255)
                                                                                                                               ttGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGctct                   >  2:1183946/1‑150 (MQ=255)
                                                                                                                                tGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCtt                  >  2:3622434/1‑150 (MQ=255)
                                                                                                                                tGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCtt                  <  2:3293332/150‑1 (MQ=255)
                                                                                                                                 ggCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCTAAAGGGCTCTTc                 <  1:1588081/150‑1 (MQ=255)
                                                                                                                                    tgctCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCata              >  1:394963/1‑150 (MQ=255)
                                                                                                                                       tCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCataaat           >  1:2545877/1‑150 (MQ=255)
                                                                                                                                         aTTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGTTACTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGTCTATTCaaaaataa         >  2:3016658/1‑150 (MQ=255)
                                                                                                                                         aTTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCataaataa         >  2:4936400/1‑150 (MQ=255)
                                                                                                                                            ggTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAAtt      >  1:1082910/1‑150 (MQ=255)
                                                                                                                                             gTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAAttt     >  2:1118199/1‑150 (MQ=255)
                                                                                                                                             gTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAAttt     <  2:1456306/150‑1 (MQ=255)
                                                                                                                                              tCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAATTTa    >  1:1398842/1‑150 (MQ=255)
                                                                                                                                                cATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAATTTAGc  >  1:2224768/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAATTTAGC  >  CP009273/4545461‑4545754

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: