Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,146,030 T→A 16.9% S194C (AGT→TGT)  yegL ← VMA domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,146,0300TA16.9% 732.3 / 132.6 273S194C (AGT→TGT) yegLVMA domain protein
Reads supporting (aligned to +/- strand):  ref base T (79/148);  new base A (17/29);  total (96/177)
Fisher's exact test for biased strand distribution p-value = 8.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 276 total aligned reads displayed.

CATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCGAAACACTT  >  CP009273/2145881‑2146179
                                                                                                                                                     |                                                                                                                                                     
 aTCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTg                                                                                                                                                      <  1:1732655/150‑1 (MQ=255)
  tCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTgg                                                                                                                                                     <  2:1955628/150‑1 (MQ=255)
    cTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGAt                                                                                                                                                   >  2:246927/1‑150 (MQ=255)
       aCAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAAc                                                                                                                                                <  2:3599825/150‑1 (MQ=255)
          aGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAg                                                                                                                                             <  2:772649/150‑1 (MQ=255)
           ggCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGc                                                                                                                                            >  2:3973971/1‑150 (MQ=255)
           ggCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGc                                                                                                                                            >  2:3629664/1‑150 (MQ=255)
            gCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCt                                                                                                                                           <  1:196902/150‑1 (MQ=255)
            gCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCt                                                                                                                                           <  1:4101082/150‑1 (MQ=255)
             cAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTa                                                                                                                                          <  1:1056174/150‑1 (MQ=255)
             cAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTa                                                                                                                                          <  1:2560161/150‑1 (MQ=255)
                  ccGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAg                                                                                                                                     <  2:85897/150‑1 (MQ=255)
                  ccGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAg                                                                                                                                     <  2:1090683/150‑1 (MQ=255)
                   cgcgCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGt                                                                                                                                    <  1:3575217/150‑1 (MQ=255)
                    gcgcTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGtt                                                                                                                                   <  2:4555567/150‑1 (MQ=255)
                    gcgcTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGtt                                                                                                                                   <  2:4555459/150‑1 (MQ=255)
                     cgcTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTc                                                                                                                                  <  2:155558/150‑1 (MQ=255)
                       cTGATGTGCGACGGTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCAc                                                                                                                                <  1:2006071/150‑1 (MQ=255)
                            gtgCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAAt                                                                                                                           >  2:2258566/1‑150 (MQ=255)
                            gtgCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAAt                                                                                                                           <  2:3544540/150‑1 (MQ=255)
                             tgCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAAtt                                                                                                                          <  2:611821/150‑1 (MQ=255)
                              gCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTg                                                                                                                         <  1:1284305/150‑1 (MQ=255)
                                  cGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAAt                                                                                                                     <  2:3525635/150‑1 (MQ=255)
                                   gTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATc                                                                                                                    >  2:1912890/1‑150 (MQ=255)
                                     tCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCt                                                                                                                  >  1:3147890/1‑150 (MQ=255)
                                          cGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCag                                                                                                             <  1:3680676/150‑1 (MQ=255)
                                              tGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGc                                                                                                         >  2:3410837/1‑150 (MQ=255)
                                               gATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCa                                                                                                        <  1:3428144/150‑1 (MQ=255)
                                                aTGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCaa                                                                                                       <  1:2191772/150‑1 (MQ=255)
                                                   cATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAgg                                                                                                    <  1:2709889/150‑1 (MQ=255)
                                                     tAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCt                                                                                                  <  1:2432679/150‑1 (MQ=255)
                                                        aCCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGAc                                                                                               <  2:3218136/150‑1 (MQ=255)
                                                         ccAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACg                                                                                              <  2:105001/150‑1 (MQ=255)
                                                          cAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGACTTGTAATCCTTGCAGAGGCAAAGGCTGACGa                                                                                             >  1:4549302/1‑150 (MQ=255)
                                                             aCGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACg                                                                                          <  2:4401986/150‑1 (MQ=255)
                                                               gCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCt                                                                                        <  1:1547843/150‑1 (MQ=255)
                                                                   gCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAAtt                                                                                    <  2:2397425/150‑1 (MQ=255)
                                                                   gCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAAtt                                                                                    <  1:3277786/150‑1 (MQ=255)
                                                                      cacaCTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTtgt                                                                                 >  1:896357/1‑150 (MQ=255)
                                                                      cacaCTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTtgt                                                                                 <  2:2057607/150‑1 (MQ=255)
                                                                        cacTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGc                                                                               <  1:300434/150‑1 (MQ=255)
                                                                            gACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGc                                                                           >  1:3939506/1‑150 (MQ=255)
                                                                            gACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGc                                                                           >  1:4259007/1‑150 (MQ=255)
                                                                              cGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGt                                                                         <  2:3141357/150‑1 (MQ=255)
                                                                              cGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGGTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGt                                                                         <  2:449751/150‑1 (MQ=255)
                                                                               gTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTc                                                                        <  2:952/150‑1 (MQ=255)
                                                                                tCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCt                                                                       >  1:1010349/1‑150 (MQ=255)
                                                                                  cAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTc                                                                     >  2:3315688/1‑150 (MQ=255)
                                                                                   aGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCa                                                                    >  2:2973652/1‑150 (MQ=255)
                                                                                    gCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCat                                                                   <  1:4054374/150‑1 (MQ=255)
                                                                                      cTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCatat                                                                 >  2:713881/1‑150 (MQ=255)
                                                                                           ggCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGcc                                                            <  1:3804205/150‑1 (MQ=255)
                                                                                               cTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCt                                                        <  2:2620208/150‑1 (MQ=255)
                                                                                                    aGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAaca                                                   <  2:4434994/150‑1 (MQ=255)
                                                                                                     gCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAacac                                                  <  1:483597/150‑1 (MQ=255)
                                                                                                      cacaACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACAcc                                                 <  1:1586331/150‑1 (MQ=255)
                                                                                                       acaACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTACCAGCGTCTTCATATCCGCCCCCTGAACACCg                                                <  2:579718/150‑1 (MQ=255)
                                                                                                        caACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGa                                               >  1:4336261/1‑150 (MQ=255)
                                                                                                          aCTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATg                                             <  1:3973999/150‑1 (MQ=255)
                                                                                                          aCTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATg                                             <  1:1804040/150‑1 (MQ=255)
                                                                                                          aCTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATg                                             <  1:3629664/150‑1 (MQ=255)
                                                                                                            ttCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGa                                           <  1:3755755/150‑1 (MQ=255)
                                                                                                             tCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGaa                                          <  1:4500859/150‑1 (MQ=255)
                                                                                                               cGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGaaaa                                        <  2:2822940/150‑1 (MQ=255)
                                                                                                                gTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAaaag                                       >  1:4201779/1‑150 (MQ=255)
                                                                                                                gTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAaaag                                       <  1:1199452/150‑1 (MQ=255)
                                                                                                                 ttCCCGTCGTGTAGCTGGAAACCGAACTCAGCGAACTGGATAATCATATAATCATTTCACTGAATTGTAATCCTTGTAGAGGCATAGTCTTACGAAAGCTATTTTGTGGCAGCGTCTTCATATCCGCCCCCTGAACTCCGATGGAaaaga                                      <  2:3066180/150‑1 (MQ=255)
                                                                                                                  tCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAaaagaa                                     <  1:970865/150‑1 (MQ=255)
                                                                                                                  tCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAaaagaa                                     >  1:2205158/1‑150 (MQ=255)
                                                                                                                   cccGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAaaagaaa                                    <  2:2234970/150‑1 (MQ=255)
                                                                                                                    ccGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAaaagaaag                                   >  1:1226942/1‑150 (MQ=255)
                                                                                                                    ccGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAaaagaaag                                   >  2:840239/1‑150 (MQ=255)
                                                                                                                     cGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGc                                  >  1:3469761/1‑150 (MQ=255)
                                                                                                                        cGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCaaa                               <  2:462022/150‑1 (MQ=255)
                                                                                                                         gTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAAc                              <  2:279505/150‑1 (MQ=255)
                                                                                                                           ggAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCt                            <  1:2904617/150‑1 (MQ=255)
                                                                                                                           ggAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCt                            <  1:2278784/150‑1 (MQ=255)
                                                                                                                           ggAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCt                            >  1:670732/1‑150 (MQ=255)
                                                                                                                            gAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCtt                           <  1:1750040/150‑1 (MQ=255)
                                                                                                                              gCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCtttt                         <  2:2121089/150‑1 (MQ=255)
                                                                                                                               cGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATCTCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTa                        <  2:4061103/150‑1 (MQ=255)
                                                                                                                                gggAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTAt                       <  1:1445311/150‑1 (MQ=255)
                                                                                                                                gggAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTAt                       >  2:1174615/1‑150 (MQ=255)
                                                                                                                                 ggAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTAtc                      >  1:3413481/1‑150 (MQ=255)
                                                                                                                                   aaaCCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTAtctt                    >  2:624360/1‑150 (MQ=255)
                                                                                                                                   aaaCCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTAtctt                    <  2:1007526/150‑1 (MQ=255)
                                                                                                                                   aaaCCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTAtctt                    <  2:479082/150‑1 (MQ=255)
                                                                                                                                    aaCCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTAtcttc                   <  2:3477017/150‑1 (MQ=255)
                                                                                                                                     aCCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTAtcttct                  <  2:1363738/150‑1 (MQ=255)
                                                                                                                                      ccGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATcttctt                 >  1:4415277/1‑150 (MQ=255)
                                                                                                                                      ccGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATcttctt                 >  1:330230/1‑150 (MQ=255)
                                                                                                                                       cGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCttcttc                >  2:4143191/1‑150 (MQ=255)
                                                                                                                                        gAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTcc               >  1:837470/1‑150 (MQ=255)
                                                                                                                                           cGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTccccc            <  2:662138/150‑1 (MQ=255)
                                                                                                                                            gCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTcccccc           >  2:2720623/1‑150 (MQ=255)
                                                                                                                                            gCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTcccccc           <  1:4482953/150‑1 (MQ=255)
                                                                                                                                             cAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCg          <  2:3940072/150‑1 (MQ=255)
                                                                                                                                              aGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCGa         >  2:1469199/1‑150 (MQ=255)
                                                                                                                                                 gAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCGAAac      >  1:261309/1‑150 (MQ=255)
                                                                                                                                                     tGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCGAAACACtt  <  2:4303282/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCGAAACACTT  >  CP009273/2145881‑2146179

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: