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,6080CA98.4% 1243.6 / 7.2 386R177I (AGA→ATA) yjiCuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base A (180/200);  total (184/202)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 386 total aligned reads displayed.

GGCTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAATTTAGC  >  CP009273/4545460‑4545754
                                                                                                                                                    |                                                                                                                                                  
 gCTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATAtt                                                                                                                                                  <  1:2970473/150‑1 (MQ=255)
  cTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCTCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATAttt                                                                                                                                                 <  2:3696637/150‑1 (MQ=255)
  cTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATAttt                                                                                                                                                 <  1:370280/150‑1 (MQ=255)
     cTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTg                                                                                                                                              <  1:4729896/150‑1 (MQ=255)
      tataCCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGc                                                                                                                                             >  1:2629366/1‑150 (MQ=255)
      tataCCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGc                                                                                                                                             <  1:1166225/150‑1 (MQ=255)
          ccAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTga                                                                                                                                         >  1:835514/1‑150 (MQ=255)
           cAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTgag                                                                                                                                        <  1:384868/150‑1 (MQ=255)
              ttattGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAg                                                                                                                                     >  2:4231821/1‑150 (MQ=255)
                     gTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAg                                                                                                                              >  1:1417167/1‑150 (MQ=255)
                     gTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAg                                                                                                                              <  2:3514945/150‑1 (MQ=255)
                       cATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAggg                                                                                                                            >  1:1239165/1‑150 (MQ=255)
                          tACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTg                                                                                                                         <  1:2729328/150‑1 (MQ=255)
                             ccAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGaaa                                                                                                                      >  1:343973/1‑150 (MQ=255)
                              cAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAat                                                                                                                     >  1:2469780/1‑150 (MQ=255)
                              cAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGATACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAat                                                                                                                     >  1:2320404/1‑150 (MQ=255)
                                aCTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAatat                                                                                                                   <  2:5038869/150‑1 (MQ=255)
                                aCTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAatat                                                                                                                   >  1:2215286/1‑150 (MQ=255)
                                   cAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATc                                                                                                                <  1:1949900/150‑1 (MQ=255)
                                    aaTATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATcc                                                                                                               <  2:1455715/150‑1 (MQ=255)
                                     atatTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCt                                                                                                              >  1:668734/1‑150 (MQ=255)
                                     atatTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCt                                                                                                              <  2:553852/150‑1 (MQ=255)
                                      tatTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTa                                                                                                             >  1:2827272/1‑150 (MQ=255)
                                        ttAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATc                                                                                                           <  1:4218470/150‑1 (MQ=255)
                                         tAGCATTAAGTTATGCGTCAAGGGCGAATCCGGAAGTATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCg                                                                                                          <  2:3619740/150‑1 (MQ=255)
                                            cATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTgg                                                                                                       >  2:84687/1‑150 (MQ=255)
                                            cATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTgg                                                                                                       >  2:5244204/1‑150 (MQ=255)
                                             aTTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGt                                                                                                      >  2:4682343/1‑150 (MQ=255)
                                              ttAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGtt                                                                                                     <  2:3324872/150‑1 (MQ=255)
                                              ttAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGtt                                                                                                     <  2:3373673/150‑1 (MQ=255)
                                              ttAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGtt                                                                                                     <  1:4636119/150‑1 (MQ=255)
                                                 aGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCg                                                                                                  <  1:106596/150‑1 (MQ=255)
                                                       gCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATaa                                                                                            <  2:1612830/150‑1 (MQ=255)
                                                         ggCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAAcc                                                                                          >  1:3610674/1‑150 (MQ=255)
                                                         ggCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAAcc                                                                                          <  2:4132573/150‑1 (MQ=255)
                                                         ggCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAAcc                                                                                          >  2:2482117/1‑150 (MQ=255)
                                                            aaGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAgg                                                                                       >  1:3350988/1‑150 (MQ=255)
                                                             aGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGa                                                                                      >  2:2286981/1‑150 (MQ=255)
                                                              gggCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGaa                                                                                     <  2:1434536/150‑1 (MQ=255)
                                                               ggCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAAc                                                                                    <  1:4959634/150‑1 (MQ=255)
                                                                  gAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACAtt                                                                                 >  1:1882840/1‑150 (MQ=255)
                                                                   aaTCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTg                                                                                >  2:3282697/1‑150 (MQ=255)
                                                                    aTCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGc                                                                               <  2:1853889/150‑1 (MQ=255)
                                                                       cGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCACTCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAAc                                                                            <  2:4476090/150‑1 (MQ=255)
                                                                       cGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAAc                                                                            <  2:5279508/150‑1 (MQ=255)
                                                                       cGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAAc                                                                            <  2:2449444/150‑1 (MQ=255)
                                                                           agaaTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGACCATTGCAACCAgg                                                                        >  2:5058636/1‑150 (MQ=255)
                                                                              aTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGcc                                                                     <  2:2228385/150‑1 (MQ=255)
                                                                              aTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGcc                                                                     >  1:1376499/1‑150 (MQ=255)
                                                                                ttGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTa                                                                   >  2:1085791/1‑150 (MQ=255)
                                                                                  gCTACCCCTGTTACCGTCATTATCTGATCATTTTAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAAt                                                                 <  2:2850751/150‑1 (MQ=255)
                                                                                  gCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAAt                                                                 <  1:2237476/150‑1 (MQ=255)
                                                                                   cTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATg                                                                <  2:2394816/150‑1 (MQ=255)
                                                                                   cTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATg                                                                <  2:3910986/150‑1 (MQ=255)
                                                                                    tACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGt                                                               >  1:3973717/1‑150 (MQ=255)
                                                                                    tACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGt                                                               >  2:3742210/1‑150 (MQ=255)
                                                                                       cccTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTAt                                                            >  1:5165115/1‑150 (MQ=255)
                                                                                       cccTGTTACCGTCATTATCTGATCAATATAGATATCAATCATTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTAt                                                            >  1:4047348/1‑150 (MQ=255)
                                                                                         cTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTAttt                                                          >  1:2579854/1‑150 (MQ=255)
                                                                                          tGTTACCGTCATTATCTGATCAATATAGATATCACTCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTc                                                         <  2:2682724/150‑1 (MQ=255)
                                                                                            ttACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCAt                                                       <  1:3147825/150‑1 (MQ=255)
                                                                                             tACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCGCAGTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATg                                                      >  1:184345/1‑150 (MQ=255)
                                                                                              aCCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGTAAGGGTTGAAATATATCATATCGTGGTTCCGGCATAACCAGGAAAATTGCAACCAGGGCCTAATGTTATTTCTTGa                                                     >  2:2820778/1‑150 (MQ=255)
                                                                                               ccGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGaa                                                    >  2:1734429/1‑150 (MQ=255)
                                                                                                cGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAAc                                                   <  1:1654839/150‑1 (MQ=255)
                                                                                                 gTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTCTCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACt                                                  >  2:2012370/1‑150 (MQ=255)
                                                                                                 gTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACt                                                  >  2:4324536/1‑150 (MQ=255)
                                                                                                   cATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGt                                                >  1:3719369/1‑150 (MQ=255)
                                                                                                     ttATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGtat                                              >  1:5072239/1‑150 (MQ=255)
                                                                                                       aTCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAg                                            <  1:927922/150‑1 (MQ=255)
                                                                                                       aTCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAg                                            <  2:598119/150‑1 (MQ=255)
                                                                                                       aTCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAg                                            >  1:73009/1‑150 (MQ=255)
                                                                                                       aTCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTATGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAAACAGGAACATTGCAACCAGGGCCTAATGTTAATTCATGAACTGTATAg                                            >  1:1760008/1‑150 (MQ=255)
                                                                                                          tGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTcc                                         <  2:769787/150‑1 (MQ=255)
                                                                                                          tGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTcc                                         >  2:1203553/1‑150 (MQ=255)
                                                                                                           gATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCt                                        <  1:2918749/150‑1 (MQ=255)
                                                                                                            aTCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTc                                       >  2:2539598/1‑150 (MQ=255)
                                                                                                            aTCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTc                                       >  1:524449/1‑150 (MQ=255)
                                                                                                               aaTATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAAc                                    >  2:4888650/1‑150 (MQ=255)
                                                                                                                  ataGATATCAATCCTTGGCTGCTCATTGGTCCATATTTATGCCTGAGAAGCTGTAAGTGTTGAAATAAATCCTATCTTGGTTCCGGAATAACCAGTAAAATAGCAACCAGGGCCTAATGTTATTTCAATAACTTTATAGTCCTCAACatg                                 >  1:2018245/1‑148 (MQ=255)
                                                                                                                     gATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATc                              >  2:3022515/1‑150 (MQ=255)
                                                                                                                     gATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATc                              <  1:428078/150‑1 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCg                             >  2:2702983/1‑150 (MQ=255)
                                                                                                                       tatCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATGTCATGAACTGTATAGTCCTCAACAGGATCGa                            >  2:179208/1‑150 (MQ=255)
                                                                                                                        atCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGaa                           <  1:4127071/150‑1 (MQ=255)
                                                                                                                          cAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAg                         <  2:4442658/150‑1 (MQ=255)
                                                                                                                          cAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTACAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAg                         >  1:4413360/1‑150 (MQ=255)
                                                                                                                           aaTCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAgg                        <  1:786743/150‑1 (MQ=255)
                                                                                                                              ccTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGct                     >  2:1333721/1‑150 (MQ=255)
                                                                                                                                  ggCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTc                 >  2:3456948/1‑150 (MQ=255)
                                                                                                                                   gctgctCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCa                >  2:1521916/1‑150 (MQ=255)
                                                                                                                                       ctCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTTTATAGTCCTCAACAGGATCGAAAGGGCTCTTCataaa            <  1:3780956/150‑1 (MQ=255)
                                                                                                                                       ctCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCataaa            >  2:1004962/1‑150 (MQ=255)
                                                                                                                                          aTTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCataaataa         >  2:3210515/1‑150 (MQ=255)
                                                                                                                                            tGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATaaataaat       <  1:4825991/150‑1 (MQ=255)
                                                                                                                                            tGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATaaataaat       <  1:558523/150‑1 (MQ=255)
                                                                                                                                            tGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATaaataaat       <  1:1144016/150‑1 (MQ=255)
                                                                                                                                             ggTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAAtt      <  2:4052645/150‑1 (MQ=255)
                                                                                                                                             ggTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAAtt      <  1:1457489/150‑1 (MQ=255)
                                                                                                                                              gTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAAttt     >  1:3895503/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGCTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAATTTAGC  >  CP009273/4545460‑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: