Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,800,129 C→T 100% W71* (TGG→TAG)  ydiY ← acid‑inducible putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,800,1290CT97.5% 657.0 / 5.4 205W71* (TGG→TAG) ydiYacid‑inducible putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (87/112);  total (91/114)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 205 total aligned reads displayed.

CTTGCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAG  >  CP009273/1799985‑1800270
                                                                                                                                                |                                                                                                                                             
cTTGCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAgg                                                                                                                                            >  1:656605/1‑148 (MQ=255)
  tGCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCCGTGGTGTGTCCATACTAGGTc                                                                                                                                          >  1:224311/1‑148 (MQ=255)
  tGCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTc                                                                                                                                          >  1:926335/1‑148 (MQ=255)
      tGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCAtgg                                                                                                                                      <  2:1683875/148‑1 (MQ=255)
          cAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATggtggt                                                                                                                                  <  1:2007182/148‑1 (MQ=255)
            aataaataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGgtgt                                                                                                                                >  2:3125253/1‑148 (MQ=255)
             ataaataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTc                                                                                                                               >  2:3229426/1‑148 (MQ=255)
                aataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGc                                                                                                                            >  2:1679403/1‑148 (MQ=255)
                aataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGc                                                                                                                            >  2:2437346/1‑148 (MQ=255)
                 atcgTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGcc                                                                                                                           <  2:715603/145‑1 (MQ=255)
                 ataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGcc                                                                                                                           >  1:835362/1‑148 (MQ=255)
                   aGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGt                                                                                                                         >  1:612794/1‑148 (MQ=255)
                    gTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTa                                                                                                                        <  1:1997386/148‑1 (MQ=255)
                     tCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTaa                                                                                                                       >  1:1640860/1‑148 (MQ=255)
                         aaTCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTga                                                                                                                   <  1:1827354/148‑1 (MQ=255)
                         aaTCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTga                                                                                                                   <  1:1827418/148‑1 (MQ=255)
                          aTCAGTTAAATTAACACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTgag                                                                                                                  <  2:2030110/148‑1 (MQ=255)
                          aTCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTgag                                                                                                                  >  1:735443/1‑148 (MQ=255)
                           tCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTgagg                                                                                                                 >  1:2706137/1‑148 (MQ=255)
                            cAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTgagga                                                                                                                >  2:3407742/1‑148 (MQ=255)
                              gTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGc                                                                                                              <  2:1579358/148‑1 (MQ=255)
                              gTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGc                                                                                                              <  1:820965/148‑1 (MQ=255)
                                 aaattaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCttt                                                                                                           >  1:1266477/1‑148 (MQ=255)
                                  aattaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCtttt                                                                                                          >  2:1907577/1‑148 (MQ=255)
                                   attaaaACGGCTACGCCCGCCCGCCGAAAATTTTTCTGAAGAACCCTCATCATAAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTc                                                                                                         >  1:1132759/1‑148 (MQ=255)
                                     taaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGt                                                                                                       <  1:3335305/148‑1 (MQ=255)
                                     taaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGt                                                                                                       <  1:973332/148‑1 (MQ=255)
                                      aaaaCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTa                                                                                                      <  2:408809/148‑1 (MQ=255)
                                       aaaCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTAt                                                                                                     <  2:2975311/148‑1 (MQ=255)
                                        aaCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTAtt                                                                                                    >  2:2075709/1‑148 (MQ=255)
                                         aCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTg                                                                                                   >  2:1583901/1‑148 (MQ=255)
                                            gCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCt                                                                                                <  1:2798836/148‑1 (MQ=255)
                                             cTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTg                                                                                               <  1:110481/148‑1 (MQ=255)
                                             cTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTg                                                                                               >  1:3504349/1‑148 (MQ=255)
                                              tACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGa                                                                                              >  1:1357611/1‑148 (MQ=255)
                                                       ccgccgCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAg                                                                                     <  1:2617241/148‑1 (MQ=255)
                                                        cgccgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGa                                                                                    <  1:1967963/148‑1 (MQ=255)
                                                        cgccgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGa                                                                                    <  1:53178/148‑1 (MQ=255)
                                                         gccgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTNGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGat                                                                                   <  1:1779858/148‑1 (MQ=255)
                                                         gccgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGat                                                                                   <  1:2275420/148‑1 (MQ=255)
                                                         gccgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGat                                                                                   <  1:694026/148‑1 (MQ=255)
                                                           cgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGatat                                                                                 <  1:40607/148‑1 (MQ=255)
                                                           cgcATATTTTTCTGAAGAACGATCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGatat                                                                                 <  2:3469560/148‑1 (MQ=255)
                                                             cATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCCAGATATcc                                                                               >  1:20123/1‑148 (MQ=255)
                                                             cATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATcc                                                                               >  1:2231630/1‑148 (MQ=255)
                                                              atatTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCt                                                                              <  2:2323166/148‑1 (MQ=255)
                                                              atatTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCt                                                                              <  2:25410/148‑1 (MQ=255)
                                                              atatTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCt                                                                              <  2:558614/148‑1 (MQ=255)
                                                               tatTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTg                                                                             >  1:2942622/1‑148 (MQ=255)
                                                                  ttttCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGt                                                                          <  2:1540942/148‑1 (MQ=255)
                                                                  ttttCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGt                                                                          <  2:1125973/148‑1 (MQ=255)
                                                                   tttCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGtt                                                                         <  1:2082179/148‑1 (MQ=255)
                                                                    ttCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTc                                                                        <  1:1807753/148‑1 (MQ=255)
                                                                    ttCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTc                                                                        <  1:27678/148‑1 (MQ=255)
                                                                      cTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAg                                                                      <  2:247381/148‑1 (MQ=255)
                                                                      cTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAg                                                                      <  1:1009847/148‑1 (MQ=255)
                                                                        gaagaaCGCTCATCAGTAGAAGAAGTATTACTGGCGTTCCCCCCCAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGAt                                                                    <  2:487680/148‑1 (MQ=255)
                                                                            aaCGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTAcc                                                                >  2:1267098/1‑148 (MQ=255)
                                                                            aaCGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTAcc                                                                >  1:401949/1‑148 (MQ=255)
                                                                            aaCGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTAcc                                                                >  1:2169984/1‑148 (MQ=255)
                                                                              cGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCtt                                                              <  1:3400412/148‑1 (MQ=255)
                                                                                   atcatTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGaaa                                                         >  2:1607603/1‑148 (MQ=255)
                                                                                     catTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAgg                                                       >  2:3261923/1‑148 (MQ=255)
                                                                                      atTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGt                                                      >  1:2588688/1‑148 (MQ=255)
                                                                                      atTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGt                                                      <  1:2454943/148‑1 (MQ=255)
                                                                                       ttAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGtt                                                     <  1:273859/148‑1 (MQ=255)
                                                                                       ttAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGtt                                                     >  2:3337426/1‑148 (MQ=255)
                                                                                        tAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGttt                                                    >  2:3185627/1‑148 (MQ=255)
                                                                                        tAGAAGAAGTATTACTGGCGGTCCCCCACAGCGACCAGGCAGTGGGGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTGCGAAAGGttt                                                    >  2:3183930/1‑148 (MQ=255)
                                                                                          gaagaaGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGttttt                                                  >  2:1187757/1‑148 (MQ=255)
                                                                                          gaagaaGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGttttt                                                  <  1:2650098/148‑1 (MQ=255)
                                                                                           aagaagTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGtttttt                                                 <  1:231224/148‑1 (MQ=255)
                                                                                            agaagTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTa                                                <  1:3125253/148‑1 (MQ=255)
                                                                                            agaagTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTa                                                <  2:636450/148‑1 (MQ=255)
                                                                                             gaagTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAg                                               >  1:1903488/1‑148 (MQ=255)
                                                                                             gaagTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAg                                               >  1:2520322/1‑148 (MQ=255)
                                                                                                    tACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGc                                        <  1:3014610/148‑1 (MQ=255)
                                                                                                    tACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGc                                        <  1:1550490/148‑1 (MQ=255)
                                                                                                       tGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAgg                                     <  2:1360407/148‑1 (MQ=255)
                                                                                                            ttCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGatcat                                >  1:2855140/1‑148 (MQ=255)
                                                                                                             tCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAgg                                     <  2:2260076/142‑1 (MQ=255)
                                                                                                              cccccACATCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTTATTGTGCATTATATCCTGCGTTCAGATTACCTTCGTTAGGTTTTTTAGCGCTGGCAGGTTCATcc                              <  2:1980420/148‑1 (MQ=255)
                                                                                                              cccccACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATcc                              >  1:2426034/1‑148 (MQ=255)
                                                                                                                cccACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCAt                            >  2:1288025/1‑148 (MQ=255)
                                                                                                                 ccACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATg                           >  1:952093/1‑148 (MQ=255)
                                                                                                                 ccACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATg                           >  2:285279/1‑148 (MQ=255)
                                                                                                                 ccACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATg                           >  2:2162000/1‑148 (MQ=255)
                                                                                                                   acaGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGAc                         >  1:1585845/1‑148 (MQ=255)
                                                                                                                    caGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACa                        <  1:596844/148‑1 (MQ=255)
                                                                                                                    caGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACa                        <  2:109420/148‑1 (MQ=255)
                                                                                                                     cgCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAg                       <  2:2846171/147‑1 (MQ=255)
                                                                                                                     aGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAg                       <  1:3330767/148‑1 (MQ=255)
                                                                                                                          ccAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTaaaa                  <  2:2894040/148‑1 (MQ=255)
                                                                                                                             ggCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACg               >  2:703716/1‑148 (MQ=255)
                                                                                                                               cAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGa             <  1:2200350/148‑1 (MQ=255)
                                                                                                                                aGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGaa            <  1:1969555/148‑1 (MQ=255)
                                                                                                                                  tggtgTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAatc          <  2:2771959/148‑1 (MQ=255)
                                                                                                                                  tggtgTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAatc          >  2:1969402/1‑148 (MQ=255)
                                                                                                                                     tgtgtCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAatcatc       <  2:2991658/148‑1 (MQ=255)
                                                                                                                                      gtgtCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCg      >  1:1533202/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                             
CTTGCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAG  >  CP009273/1799985‑1800270

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: