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.9% 610.1 / 2.6 189W71* (TGG→TAG) ydiYacid‑inducible putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (85/100);  total (88/101)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 189 total aligned reads displayed.

GCTTGCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTG  >  CP009273/1799984‑1800273
                                                                                                                                                 |                                                                                                                                                
        gACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCAtggt                                                                                                                                        >  1:1593250/1‑148 (MQ=255)
        gACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCAtggt                                                                                                                                        >  1:1595918/1‑148 (MQ=255)
        gACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCAtggt                                                                                                                                        >  1:2202112/1‑148 (MQ=255)
         aCCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCAtggtg                                                                                                                                       <  2:2313416/148‑1 (MQ=255)
            acataaataGTCATAAGCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGgtggtg                                                                                                                                    <  2:3177808/146‑1 (MQ=255)
               taaataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTcg                                                                                                                                 >  1:2271709/1‑147 (MQ=255)
                  ataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGcc                                                                                                                              >  2:3565046/1‑148 (MQ=255)
                   taGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCg                                                                                                                             >  1:2422490/1‑148 (MQ=255)
                   taGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCg                                                                                                                             <  1:1455564/148‑1 (MQ=255)
                    aGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGt                                                                                                                            >  2:1764631/1‑148 (MQ=255)
                      tCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTaa                                                                                                                          <  2:907495/148‑1 (MQ=255)
                      tCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTaa                                                                                                                          <  1:3147094/148‑1 (MQ=255)
                      tCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTaa                                                                                                                          <  1:898129/148‑1 (MQ=255)
                      tCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTaa                                                                                                                          <  2:1460337/148‑1 (MQ=255)
                      tCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTaa                                                                                                                          >  2:2337402/1‑148 (MQ=255)
                        aTAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGt                                                                                                                        <  2:1718310/148‑1 (MQ=255)
                              atttaaattaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCCCAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTgaggag                                                                                                                  <  1:1309841/146‑1 (MQ=255)
                               gTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGc                                                                                                                 <  2:2880765/148‑1 (MQ=255)
                               gTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGc                                                                                                                 <  2:202864/148‑1 (MQ=255)
                                   aattaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCtttt                                                                                                             <  1:3098765/148‑1 (MQ=255)
                                   aattaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCtttt                                                                                                             <  1:678951/148‑1 (MQ=255)
                                   aattaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCtttt                                                                                                             >  1:2367191/1‑148 (MQ=255)
                                    attaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTc                                                                                                            >  1:122757/1‑148 (MQ=255)
                                     ttaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCg                                                                                                           >  1:550079/1‑148 (MQ=255)
                                      taaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGt                                                                                                          <  1:3063400/148‑1 (MQ=255)
                                      taaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGt                                                                                                          >  1:2851854/1‑148 (MQ=255)
                                       aacaCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTa                                                                                                         <  2:1014462/145‑1 (MQ=255)
                                        aaaCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTAt                                                                                                        >  2:3461972/1‑148 (MQ=255)
                                          aCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTg                                                                                                      <  1:2200205/148‑1 (MQ=255)
                                           cGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGc                                                                                                     <  1:2337402/148‑1 (MQ=255)
                                           cGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGc                                                                                                     >  1:2815599/1‑148 (MQ=255)
                                           cGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGc                                                                                                     <  1:3233346/148‑1 (MQ=255)
                                              cTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCAACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTg                                                                                                  >  2:3319707/1‑148 (MQ=255)
                                               tccgcccgcccgccGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGa                                                                                                 <  2:2017768/146‑1 (MQ=255)
                                               tACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGa                                                                                                 <  1:1921717/148‑1 (MQ=255)
                                                aCGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGAt                                                                                                <  1:572040/148‑1 (MQ=255)
                                                    ccgcccgccGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATtgtg                                                                                            <  1:1082528/148‑1 (MQ=255)
                                                        ccgccgCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAg                                                                                        <  2:1533532/148‑1 (MQ=255)
                                                         cgccgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGa                                                                                       >  2:2580340/1‑148 (MQ=255)
                                                          gccgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCTTGGTGGTGTCAGCCGTATGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGat                                                                                      <  2:3429398/148‑1 (MQ=255)
                                                          gccgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGat                                                                                      >  1:1041916/1‑148 (MQ=255)
                                                           ccgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGata                                                                                     <  1:1764631/148‑1 (MQ=255)
                                                           ccgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGata                                                                                     <  1:2736467/148‑1 (MQ=255)
                                                            cgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGatat                                                                                    >  2:3099758/1‑148 (MQ=255)
                                                             gcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATc                                                                                   <  2:320963/148‑1 (MQ=255)
                                                                 atTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGc                                                                               >  2:2122320/1‑148 (MQ=255)
                                                                  tttttCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCg                                                                              >  1:3155572/1‑148 (MQ=255)
                                                                  tttttCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCg                                                                              <  2:467012/148‑1 (MQ=255)
                                                                  tttttCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCg                                                                              >  1:1440321/1‑148 (MQ=255)
                                                                   ttttCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGt                                                                             >  1:949411/1‑148 (MQ=255)
                                                                       cTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAg                                                                         <  1:2865229/148‑1 (MQ=255)
                                                                        tGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGa                                                                        <  2:2822182/148‑1 (MQ=255)
                                                                        tGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGATGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGa                                                                        >  2:1786550/1‑148 (MQ=255)
                                                                           agaaCGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCTACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTTTTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTa                                                                     <  2:1176707/148‑1 (MQ=255)
                                                                           agaaCGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTa                                                                     >  1:2961858/1‑148 (MQ=255)
                                                                             aaCGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTAcc                                                                   <  1:2310094/148‑1 (MQ=255)
                                                                               cGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCtt                                                                 >  2:796684/1‑148 (MQ=255)
                                                                               cGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCtt                                                                 >  2:3419504/1‑148 (MQ=255)
                                                                                   catcatTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGaa                                                             <  2:3340206/148‑1 (MQ=255)
                                                                                   catcatTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGaa                                                             >  1:1687956/1‑148 (MQ=255)
                                                                                     tcatTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAg                                                           <  1:1293491/148‑1 (MQ=255)
                                                                                     tcatTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAg                                                           <  1:1293494/148‑1 (MQ=255)
                                                                                       atTAGAAGAAGTATTACTGGCTTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGt                                                         >  2:1425735/1‑148 (MQ=255)
                                                                                         tAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGttt                                                       <  1:500641/148‑1 (MQ=255)
                                                                                         tAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGttt                                                       <  2:1244729/148‑1 (MQ=255)
                                                                                         tAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGttt                                                       >  1:2990223/1‑148 (MQ=255)
                                                                                           gaagaaGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGttttt                                                     >  1:753556/1‑148 (MQ=255)
                                                                                           gaagaaGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGttttt                                                     <  2:1527726/148‑1 (MQ=255)
                                                                                            aagaagTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGtttttt                                                    <  2:1709240/148‑1 (MQ=255)
                                                                                            aagaagTAGTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGtttttt                                                    <  2:2896032/148‑1 (MQ=255)
                                                                                              gaagTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAg                                                  <  2:3274003/148‑1 (MQ=255)
                                                                                              gaagTATTACTGGCGTGCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAg                                                  <  2:3214132/148‑1 (MQ=255)
                                                                                               aagTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAgc                                                 >  1:153038/1‑148 (MQ=255)
                                                                                               aagTATTACTGGCGGTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAgc                                                 >  1:1956714/1‑148 (MQ=255)
                                                                                                 gTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAgcgc                                               >  1:446805/1‑148 (MQ=255)
                                                                                                     tACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGc                                           >  2:1639035/1‑148 (MQ=255)
                                                                                                        tGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAgg                                        >  1:7154/1‑148 (MQ=255)
                                                                                                         ggCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGa                                       <  1:1700604/148‑1 (MQ=255)
                                                                                                         ggCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGa                                       <  1:398200/148‑1 (MQ=255)
                                                                                                          gCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGat                                      >  2:2963621/1‑148 (MQ=255)
                                                                                                          gCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGat                                      <  2:724415/148‑1 (MQ=255)
                                                                                                          gCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGat                                      <  1:258403/148‑1 (MQ=255)
                                                                                                          gCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGat                                      <  1:257961/148‑1 (MQ=255)
                                                                                                          gCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGat                                      >  2:1412137/1‑148 (MQ=255)
                                                                                                                ccccACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCa                                <  2:2664605/148‑1 (MQ=255)
                                                                                                                  ccACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATg                              >  2:2795683/1‑148 (MQ=255)
                                                                                                                     caGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACa                           <  1:2937620/148‑1 (MQ=255)
                                                                                                                      aGCGACCAGGCAGTGGTGTGTCCATTCTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAg                          <  2:578663/148‑1 (MQ=255)
                                                                                                                        cGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTa                        >  2:3045271/1‑148 (MQ=255)
                                                                                                                           ccAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTaaaa                     <  2:2450999/148‑1 (MQ=255)
                                                                                                                            cAGGCAGTGGTGTGTCCATACTAGGTCATGGTTGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCCGGATCATCCATGACAGTaaaaa                    >  2:2511071/1‑148 (MQ=255)
                                                                                                                            cAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTaaaaa                    <  1:3565046/148‑1 (MQ=255)
                                                                                                                               gCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACgg                 >  1:2714772/1‑148 (MQ=255)
                                                                                                                                    ggtgTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAatca            <  1:173665/148‑1 (MQ=255)
                                                                                                                                     gtgtgtCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAatcat           >  2:1944077/1‑148 (MQ=255)
                                                                                                                                      tgtgtCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAatcatc          <  1:432606/148‑1 (MQ=255)
                                                                                                                                       gtgtCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCg         >  2:3531746/1‑148 (MQ=255)
                                                                                                                                          tCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTCTTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATTACAGTAAAAACGGAATCATCGGCa      <  1:2543192/148‑1 (MQ=255)
                                                                                                                                             aTACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCt   >  1:854428/1‑148 (MQ=255)
                                                                                                                                              tACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTg  >  1:3463722/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCTTGCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTG  >  CP009273/1799984‑1800273

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: