Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,895,069 C→T 11.3% G258G (GGG→GGA yoaE ← putative membrane protein/conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,895,0690CT11.3% 1007.3 / 85.3 317G258G (GGG→GGAyoaEputative membrane protein/conserved protein
Reads supporting (aligned to +/- strand):  ref base C (133/148);  new base T (18/18);  total (151/166)
Fisher's exact test for biased strand distribution p-value = 8.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 325 total aligned reads displayed.

CGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAAT  >  CP009273/1894921‑1895218
                                                                                                                                                    |                                                                                                                                                     
 gcgcgGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTccc                                                                                                                                                     <  1:69198/150‑1 (MQ=255)
  cgcgGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCg                                                                                                                                                    <  1:3463435/150‑1 (MQ=255)
  cgcgGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCg                                                                                                                                                    >  1:1320407/1‑150 (MQ=255)
  cgcgGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCg                                                                                                                                                    >  1:1751071/1‑150 (MQ=255)
    cggcgTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGcc                                                                                                                                                  <  2:193712/150‑1 (MQ=255)
     ggcgTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGccc                                                                                                                                                 >  1:787388/1‑150 (MQ=255)
      gcgTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTTCCGCCCa                                                                                                                                                <  1:2838991/150‑1 (MQ=255)
      gcgTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCa                                                                                                                                                <  2:479525/150‑1 (MQ=255)
            tGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAAc                                                                                                                                          <  2:112459/150‑1 (MQ=255)
                cccGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTACTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTTCCGCCCATCAAACGCAg                                                                                                                                      >  1:466611/1‑150 (MQ=255)
                   gcgcAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTTCCGCCCATCAAACGCAGGAt                                                                                                                                   <  1:2965687/150‑1 (MQ=255)
                   gcgcAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGAt                                                                                                                                   <  1:2157989/150‑1 (MQ=255)
                     gcAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTTCCGCCCATCAAACGCAGGATgg                                                                                                                                 <  2:2596044/150‑1 (MQ=255)
                       agagAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCa                                                                                                                               >  2:391029/1‑150 (MQ=255)
                        gagaACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCAt                                                                                                                              <  1:1197212/150‑1 (MQ=255)
                        gagaACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCAt                                                                                                                              >  2:1247243/1‑150 (MQ=255)
                         agaACGCGACGCCAGCGTCAGTACGCCGTTAATCAAGTAACGATATTCTTCGGCAAAAGCACCTTCCGGGATCGGCATCGTCATCGGGAGATCGGCATCGAGCTTAACAATGGCCTGACGTTTCCCGCCCATCAAACACAGGATGGCaac                                                                                                                             >  1:2824731/1‑148 (MQ=255)
                              gcgACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAg                                                                                                                        <  2:2662583/150‑1 (MQ=255)
                                 aCGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTAc                                                                                                                     >  1:3572688/1‑150 (MQ=255)
                                  cGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTAcg                                                                                                                    >  2:3290213/1‑150 (MQ=255)
                                      aGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCt                                                                                                                <  2:855006/150‑1 (MQ=255)
                                           cAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCTGTACGCGCTCGCAg                                                                                                           <  1:3395762/150‑1 (MQ=255)
                                           cAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAg                                                                                                           >  1:4502879/1‑150 (MQ=255)
                                            aGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGc                                                                                                          >  2:1306963/1‑150 (MQ=255)
                                             gTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCg                                                                                                         <  1:3089940/150‑1 (MQ=255)
                                             gTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCg                                                                                                         <  1:3489022/150‑1 (MQ=255)
                                              tACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCgg                                                                                                        <  2:4426202/150‑1 (MQ=255)
                                               aCGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGc                                                                                                       <  1:3985404/150‑1 (MQ=255)
                                               aCGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGc                                                                                                       <  1:3856625/150‑1 (MQ=255)
                                                cGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCa                                                                                                      >  1:2891809/1‑150 (MQ=255)
                                                   cGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAg                                                                                                   <  2:3045482/150‑1 (MQ=255)
                                                    gTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGt                                                                                                  <  1:4315611/150‑1 (MQ=255)
                                                       aaTCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGa                                                                                               >  1:3660243/1‑150 (MQ=255)
                                                        aTCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGAc                                                                                              <  1:2363388/150‑1 (MQ=255)
                                                         tCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACt                                                                                             <  1:3329317/150‑1 (MQ=255)
                                                          cATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACtg                                                                                            >  1:3231045/1‑150 (MQ=255)
                                                               aaCGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTTCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGc                                                                                       >  2:2025056/1‑150 (MQ=255)
                                                               aaCGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGc                                                                                       <  2:4503410/150‑1 (MQ=255)
                                                               aaCGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGc                                                                                       <  2:2172088/150‑1 (MQ=255)
                                                                    tcttcttcGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTTCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATa                                                                                  <  2:1209749/150‑1 (MQ=255)
                                                                     cttcttcGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTTCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATaa                                                                                 >  2:3793012/1‑150 (MQ=255)
                                                                      ttcttcGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATaaa                                                                                <  2:4608253/150‑1 (MQ=255)
                                                                         ttcGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGtt                                                                             >  2:4097755/1‑150 (MQ=255)
                                                                          tcGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTg                                                                            >  1:3183603/1‑150 (MQ=255)
                                                                          tcGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTg                                                                            <  2:1246577/150‑1 (MQ=255)
                                                                           cGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGc                                                                           <  2:3469901/150‑1 (MQ=255)
                                                                              cAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGAc                                                                        <  2:4193971/150‑1 (MQ=255)
                                                                                aaTGCACCTTCCGGGATCGGCATCGGCATCNGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGAcgc                                                                      >  1:4280021/1‑150 (MQ=255)
                                                                                aaTGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGAcgc                                                                      >  1:1276246/1‑150 (MQ=255)
                                                                                   gCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCa                                                                   >  2:3346014/1‑150 (MQ=255)
                                                                                     aCCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAAt                                                                 >  2:3569811/1‑150 (MQ=255)
                                                                                     aCCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAAt                                                                 <  2:3150557/150‑1 (MQ=255)
                                                                                       cTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCt                                                               <  2:1491045/150‑1 (MQ=255)
                                                                                          ccGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGt                                                            <  2:4036754/150‑1 (MQ=255)
                                                                                          ccGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGt                                                            >  2:3458671/1‑150 (MQ=255)
                                                                                            gggATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTa                                                          >  1:2531936/1‑150 (MQ=255)
                                                                                               atcggcatcggcATCGGGTTATCGGCATAGTGCTGAACATAGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGAACGCGCTCGAAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAac                                                       >  2:1634978/1‑150 (MQ=255)
                                                                                                tcggcatcggcaTCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAaca                                                      <  1:19200/150‑1 (MQ=255)
                                                                                                tcggcatcggcaTCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAaca                                                      >  1:1455127/1‑150 (MQ=255)
                                                                                                 cggcatcggcatCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAacac                                                     <  2:558785/150‑1 (MQ=255)
                                                                                                  ggcatcggcatcGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACt                                                    >  2:2937789/1‑150 (MQ=255)
                                                                                                    caacggcaacggGTTATCGGCATCGTGCAGAACATTGGCCTGACGATTTCCGCCCATCAAACGCAGGAAGGCATCGGCAGTACGCGCTCGAAGCGCCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTc                                                  >  1:4587994/4‑150 (MQ=255)
                                                                                                     atcggcatcggGTTATCAGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCg                                                 >  1:100413/1‑150 (MQ=255)
                                                                                                      tcggcatcggGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCga                                                <  1:3744619/150‑1 (MQ=255)
                                                                                                       cggcatcggGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCgat                                               <  1:2306339/150‑1 (MQ=255)
                                                                                                           atcggGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCgatgatg                                           <  1:55120/150‑1 (MQ=255)
                                                                                                           atcggGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCgatgatg                                           >  1:2194710/1‑150 (MQ=255)
                                                                                                             cggGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCgatgatgat                                         <  2:387595/150‑1 (MQ=255)
                                                                                                              gggTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATc                                        <  1:4169689/150‑1 (MQ=255)
                                                                                                              gggTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCACGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGAAGCGCAATCTGGTTAAACACTTCGATGATGATc                                        >  1:2050223/1‑150 (MQ=255)
                                                                                                                  tATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCgaga                                    <  1:1492594/150‑1 (MQ=255)
                                                                                                                  tATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCgaga                                    >  2:3701332/1‑150 (MQ=255)
                                                                                                                     cGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGc                                 >  1:1415887/1‑150 (MQ=255)
                                                                                                                     cGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGc                                 >  2:305260/1‑150 (MQ=255)
                                                                                                                      ggCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGcc                                <  1:2125116/150‑1 (MQ=255)
                                                                                                                       gCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCa                               <  2:4280316/150‑1 (MQ=255)
                                                                                                                          tCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATc                            <  2:2238730/150‑1 (MQ=255)
                                                                                                                           cGTGCTGAACATTGGCCTGACGTTTTCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATcg                           <  1:4450862/150‑1 (MQ=255)
                                                                                                                              gCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCgg                        >  1:4157093/1‑150 (MQ=255)
                                                                                                                                tGAACATTGGCCTGACGTTTTCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCa                      >  1:2428396/1‑150 (MQ=255)
                                                                                                                                tGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCa                      <  2:3207351/150‑1 (MQ=255)
                                                                                                                                tGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCa                      >  2:1954913/1‑150 (MQ=255)
                                                                                                                                tGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATAGAGAAGCCAATCGCGGCa                      >  2:2447830/1‑150 (MQ=255)
                                                                                                                                  aaCATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATa                    >  1:960263/1‑150 (MQ=255)
                                                                                                                                  aaCATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATa                    <  2:3471022/150‑1 (MQ=255)
                                                                                                                                   aCATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCACAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATAc                   >  1:1635014/1‑150 (MQ=255)
                                                                                                                                   aCATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATAc                   >  1:3025012/1‑150 (MQ=255)
                                                                                                                                     aTTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAg                 <  1:1946980/150‑1 (MQ=255)
                                                                                                                                      ttGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATTATGATCGAGAAGCCAATCGCGGCATACAgg                <  2:227678/150‑1 (MQ=255)
                                                                                                                                         gCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTaa             >  1:2577326/1‑150 (MQ=255)
                                                                                                                                           cTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAAcc           <  2:1519747/150‑1 (MQ=255)
                                                                                                                                            tGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCt          >  2:1362639/1‑150 (MQ=255)
                                                                                                                                            tGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCt          >  1:425795/1‑150 (MQ=255)
                                                                                                                                             gACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCtt         >  2:961298/1‑150 (MQ=255)
                                                                                                                                                gTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCg      <  2:2901177/150‑1 (MQ=255)
                                                                                                                                                gTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCg      <  2:2210235/150‑1 (MQ=255)
                                                                                                                                                gTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCg      >  2:1809463/1‑150 (MQ=255)
                                                                                                                                                  ttCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGa    >  2:2840955/1‑150 (MQ=255)
                                                                                                                                                  ttCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGa    >  1:4083462/1‑150 (MQ=255)
                                                                                                                                                    cccGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAAt  >  2:610284/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAAT  >  CP009273/1894921‑1895218

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: