Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,976,807 C→A 13.3% H164N (CAT→AAT)  aslB → putative AslA‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,976,8070CA13.3% 514.8 / 38.6 161H164N (CAT→AAT) aslBputative AslA‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base C (72/67);  new base A (6/16);  total (78/83)
Fisher's exact test for biased strand distribution p-value = 3.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

TGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGCCTGGCTGATGAAACAACTGCCCGCG  >  CP009273/3976660‑3976954
                                                                                                                                                   |                                                                                                                                                   
tGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAAcatc                                                                                                                                                     <  1:2069308/148‑1 (MQ=255)
tGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAAcatc                                                                                                                                                     <  1:2069311/148‑1 (MQ=255)
tGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAAcatc                                                                                                                                                     >  2:473875/1‑148 (MQ=255)
  gTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAAcataat                                                                                                                                                   <  2:2166815/148‑1 (MQ=255)
   tGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATTCGTGCCCTGACGCTCCTGCAAAAACATAATg                                                                                                                                                  <  2:1706717/148‑1 (MQ=255)
   tGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATAATg                                                                                                                                                  <  2:1706691/148‑1 (MQ=255)
    gcgcGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGt                                                                                                                                                 <  1:2434016/148‑1 (MQ=255)
    gcgcGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATAATGt                                                                                                                                                 <  2:2023804/148‑1 (MQ=255)
      gcgTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATAATGTCg                                                                                                                                               <  1:2719894/148‑1 (MQ=255)
          ttCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACta                                                                                                                                           >  2:2852237/1‑148 (MQ=255)
            ctcGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACtata                                                                                                                                         <  1:3210962/148‑1 (MQ=255)
             tcGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACtatat                                                                                                                                        <  1:3042482/148‑2 (MQ=255)
                 ggAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTg                                                                                                                                    <  1:136466/148‑1 (MQ=255)
                     catcatTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTCTAACGTGCTgg                                                                                                                                >  2:1687083/1‑148 (MQ=255)
                      atcatTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGt                                                                                                                               >  2:194663/1‑148 (MQ=255)
                        catTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGAGTGCCCTCACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCt                                                                                                                             >  1:2825789/1‑148 (MQ=255)
                         atTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTg                                                                                                                            >  2:2325737/1‑148 (MQ=255)
                          ttttCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGc                                                                                                                           >  1:2312074/1‑148 (MQ=255)
                           tttCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCg                                                                                                                          >  1:1527894/1‑148 (MQ=255)
                           tttCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCg                                                                                                                          >  2:700149/1‑148 (MQ=255)
                           tttCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCg                                                                                                                          <  2:2109413/148‑1 (MQ=255)
                           tttCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATAATGTCGACTATAACGTGCTGGTCTGCg                                                                                                                          <  1:1673997/148‑1 (MQ=255)
                            ttCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGt                                                                                                                         >  2:600535/1‑148 (MQ=255)
                              cTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTa                                                                                                                       <  1:2024610/148‑1 (MQ=255)
                              cTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTa                                                                                                                       <  1:2884737/148‑1 (MQ=255)
                              cTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTa                                                                                                                       >  2:194080/1‑148 (MQ=255)
                              cTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTa                                                                                                                       <  1:1234779/148‑1 (MQ=255)
                              cTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTa                                                                                                                       <  2:292769/148‑1 (MQ=255)
                                   tgGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGc                                                                                                                  >  1:1924890/1‑148 (MQ=255)
                                    ggTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCa                                                                                                                 >  2:352561/1‑148 (MQ=255)
                                     gTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCAc                                                                                                                >  1:1166585/1‑148 (MQ=255)
                                     gTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATAATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCAc                                                                                                                <  2:1807125/148‑1 (MQ=255)
                                          tCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAgcg                                                                                                           >  1:51590/1‑148 (MQ=255)
                                                 aTGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCgcagca                                                                                                    >  2:1478463/1‑148 (MQ=255)
                                                 aTGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCgcagca                                                                                                    <  2:1241784/148‑1 (MQ=255)
                                                  tGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCaa                                                                                                   <  1:2582948/148‑1 (MQ=255)
                                                  tGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATAATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCaa                                                                                                   >  1:843731/1‑148 (MQ=255)
                                                     cccGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCg                                                                                                <  1:2451377/148‑1 (MQ=255)
                                                        gCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTg                                                                                             >  2:3138901/1‑148 (MQ=255)
                                                          cTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCa                                                                                           <  2:567911/148‑1 (MQ=255)
                                                            gagaTCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAgg                                                                                         >  1:267093/1‑148 (MQ=255)
                                                            gagaTCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAgg                                                                                         >  1:266351/1‑148 (MQ=255)
                                                             agaTCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGt                                                                                        >  2:3012395/1‑148 (MQ=255)
                                                               aTCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGtct                                                                                      >  1:956012/1‑146 (MQ=255)
                                                                tcaacaATCAATATCTCGTGAATAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTACTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGtata                                                                                     <  2:1851312/145‑1 (MQ=255)
                                                                  cacaATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATg                                                                                   <  1:1997804/148‑1 (MQ=255)
                                                                    caatcaatATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATAATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGAt                                                                                 <  1:1265855/148‑1 (MQ=255)
                                                                      atcaatATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGAttt                                                                               <  1:2291866/148‑1 (MQ=255)
                                                                       tcaatATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGAtttt                                                                              >  2:1943747/1‑148 (MQ=255)
                                                                       tcaatATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGAtttt                                                                              >  1:2715817/1‑148 (MQ=255)
                                                                       tcaatATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGAtttt                                                                              >  2:1943578/1‑148 (MQ=255)
                                                                           tatCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATAATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTtgt                                                                          >  2:1903335/1‑148 (MQ=255)
                                                                            atCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTtgtg                                                                         <  1:352564/148‑1 (MQ=255)
                                                                            atCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTtgtg                                                                         >  2:1475091/1‑148 (MQ=255)
                                                                            atCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATAATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTtgtg                                                                         <  2:295915/148‑1 (MQ=255)
                                                                              cgcgTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCg                                                                       >  1:477970/1‑148 (MQ=255)
                                                                              cgcgTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCg                                                                       >  1:2056620/1‑148 (MQ=255)
                                                                                cgTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGAt                                                                     >  1:1968147/1‑148 (MQ=255)
                                                                                     cTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATAATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCggg                                                                <  1:2594320/148‑1 (MQ=255)
                                                                                      tAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGa                                                               >  2:965730/1‑148 (MQ=255)
                                                                                      tAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGa                                                               >  2:965697/1‑148 (MQ=255)
                                                                                          ggTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCg                                                           >  1:1625079/1‑148 (MQ=255)
                                                                                          ggTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCg                                                           <  1:2008806/148‑1 (MQ=255)
                                                                                           gTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGa                                                          <  2:1883776/148‑1 (MQ=255)
                                                                                           gTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATAATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGa                                                          <  2:356216/148‑1 (MQ=255)
                                                                                                aGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCa                                                     >  1:220322/1‑148 (MQ=255)
                                                                                                  aCCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATc                                                   >  2:2917168/1‑148 (MQ=255)
                                                                                                     cACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAg                                                >  1:1087023/1‑148 (MQ=255)
                                                                                                     cACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAg                                                >  1:3096402/1‑148 (MQ=255)
                                                                                                      aCGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATAATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGt                                               <  1:1107068/148‑1 (MQ=255)
                                                                                                         cATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAACACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTa                                            <  1:993696/148‑1 (MQ=255)
                                                                                                         cATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTa                                            <  1:994466/148‑1 (MQ=255)
                                                                                                         cATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTa                                            <  2:2579081/148‑1 (MQ=255)
                                                                                                           tAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTAtt                                          >  2:1606653/1‑148 (MQ=255)
                                                                                                           tAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTAtt                                          <  2:2253391/148‑1 (MQ=255)
                                                                                                            aaGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCACCCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTc                                         <  2:614366/148‑1 (MQ=255)
                                                                                                            aaGCTGGTGATGCGTGCCCAGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGAGCGAGTCGAATTCATCCAGTTTATTc                                         >  2:2380780/1‑148 (MQ=255)
                                                                                                               cTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCgg                                      <  1:600535/148‑1 (MQ=255)
                                                                                                               cTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATAATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCgg                                      >  1:3119709/1‑148 (MQ=255)
                                                                                                                tGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCggt                                     <  2:1994165/148‑1 (MQ=255)
                                                                                                                tGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCggt                                     <  2:1819097/148‑1 (MQ=255)
                                                                                                                   tGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCggtggt                                  >  2:2349717/1‑148 (MQ=255)
                                                                                                                         gTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAgcg                            <  2:982603/148‑1 (MQ=255)
                                                                                                                         gTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAgcg                            <  2:982506/148‑1 (MQ=255)
                                                                                                                          tGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAgcgc                           <  2:248366/148‑1 (MQ=255)
                                                                                                                           gCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGcc                          <  2:1166585/148‑1 (MQ=255)
                                                                                                                           gCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGcc                          <  2:2801559/148‑1 (MQ=255)
                                                                                                                              cTGACGCTCCTGCAAAAACATAATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGCCTgg                       <  2:843731/148‑1 (MQ=255)
                                                                                                                              cTGACGCTCCTGCAAAAACATAATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGCCTgg                       <  2:844676/148‑1 (MQ=255)
                                                                                                                                  cGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGCCTGGCtga                   <  2:94850/148‑1 (MQ=255)
                                                                                                                                  cGCTCCTGCAAAAACATAATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGCCTGGCtga                   <  2:2235619/148‑1 (MQ=255)
                                                                                                                                    cTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGCCTGGCtgatg                 >  1:2370543/1‑148 (MQ=255)
                                                                                                                                         gCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGCCTGGCTGATGAaaca            <  2:396194/148‑1 (MQ=255)
                                                                                                                                         gCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGCCTGGCTGATGAaaca            <  2:397531/148‑1 (MQ=255)
                                                                                                                                         gCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGCCTGGCTGATGAaaca            <  1:1186604/148‑1 (MQ=255)
                                                                                                                                          cAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGCCTGGCTGATGAaacaa           <  2:306486/148‑1 (MQ=255)
                                                                                                                                          cAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGCCTGGCTGATGAaacaa           >  1:2193966/1‑148 (MQ=255)
                                                                                                                                              aaCATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGCCTGGCTGATGAAACAACTGc       >  2:938252/1‑148 (MQ=255)
                                                                                                                                                catcatGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGCCTGGCTGATGAAACAACTGccc     >  2:630757/1‑148 (MQ=255)
                                                                                                                                                  tcatGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTATTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGCCTGGCTGATTAAACAACTGCCcg  >  1:2407113/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGCCTGGCTGATGAAACAACTGCCCGCG  >  CP009273/3976660‑3976954

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: