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,1290CT99.1% 1049.3 / ‑0.8 323W71* (TGG→TAG) ydiYacid‑inducible putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (181/138);  total (183/140)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 323 total aligned reads displayed.

CAGCTTGCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGC  >  CP009273/1799982‑1800274
                                                                                                                                                   |                                                                                                                                                 
     tGCCTGACCAAATAAATAGTCATAATCATTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCAt                                                                                                                                            <  1:4086063/150‑1 (MQ=255)
      gCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCAtg                                                                                                                                           <  1:292110/150‑1 (MQ=255)
          gACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCAtggtgg                                                                                                                                       >  2:271774/1‑150 (MQ=255)
           aCCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATggtggt                                                                                                                                      <  1:423475/150‑1 (MQ=255)
            ccAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGgtggtg                                                                                                                                     >  1:482389/1‑150 (MQ=255)
             cAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGgtgt                                                                                                                                    >  2:5312528/1‑150 (MQ=255)
              aaataaatAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTc                                                                                                                                   >  1:1012297/1‑150 (MQ=255)
               aataaataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCa                                                                                                                                  <  1:3537744/150‑1 (MQ=255)
               aataaataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCa                                                                                                                                  <  1:697365/150‑1 (MQ=255)
                  aaataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGcc                                                                                                                               <  1:2598666/150‑1 (MQ=255)
                   aataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCg                                                                                                                              >  2:1377722/1‑150 (MQ=255)
                    ataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGt                                                                                                                             <  2:3601908/150‑1 (MQ=255)
                    ataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGt                                                                                                                             <  2:927022/150‑1 (MQ=255)
                      aGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTaa                                                                                                                           >  2:3196239/1‑150 (MQ=255)
                      aGTCATAATCAGATAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTaa                                                                                                                           >  2:2295244/1‑150 (MQ=255)
                       gTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAg                                                                                                                          >  1:2332032/1‑150 (MQ=255)
                       gTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAg                                                                                                                          <  1:2506585/150‑1 (MQ=255)
                         cATAATCAGTTAAAATAACACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTg                                                                                                                        <  1:3660820/150‑1 (MQ=255)
                          aTAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACCGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTga                                                                                                                       <  2:411025/150‑1 (MQ=255)
                          aTAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTga                                                                                                                       <  1:1281734/150‑1 (MQ=255)
                           tAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTgag                                                                                                                      <  1:5244381/150‑1 (MQ=255)
                           tAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTgag                                                                                                                      <  2:3860579/150‑1 (MQ=255)
                            aaTCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTgagg                                                                                                                     <  2:1790963/150‑1 (MQ=255)
                            aaTCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTgagg                                                                                                                     <  1:5035104/150‑1 (MQ=255)
                            aaTCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTgagg                                                                                                                     >  1:2072644/1‑150 (MQ=255)
                            aaTCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTgagg                                                                                                                     >  1:4928784/1‑150 (MQ=255)
                             aTCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTgagga                                                                                                                    <  2:437943/150‑1 (MQ=255)
                              tCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTgaggag                                                                                                                   >  2:2353127/1‑150 (MQ=255)
                               cAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGc                                                                                                                  >  2:4953315/1‑150 (MQ=255)
                               cAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGc                                                                                                                  >  2:350866/1‑150 (MQ=255)
                                aGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTCGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCCGCCGTGAGTGAGGAGCt                                                                                                                 >  2:4501295/1‑150 (MQ=255)
                                  ttaaattaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCttt                                                                                                               <  2:4665071/150‑1 (MQ=255)
                                    aaattaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTc                                                                                                             >  2:1579428/1‑150 (MQ=255)
                                     aattaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCg                                                                                                            >  1:3120836/1‑150 (MQ=255)
                                     aattaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCg                                                                                                            >  1:900441/1‑150 (MQ=255)
                                      attaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGt                                                                                                           <  2:4771833/150‑1 (MQ=255)
                                        taaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTAt                                                                                                         <  2:1864450/150‑1 (MQ=255)
                                        taaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTAt                                                                                                         >  1:15983/1‑150 (MQ=255)
                                         aaaaCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTAtt                                                                                                        >  1:3218600/1‑150 (MQ=255)
                                          aaaCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTg                                                                                                       >  1:4431497/1‑150 (MQ=255)
                                          aaaCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTg                                                                                                       <  1:2012577/150‑1 (MQ=255)
                                           aaCGGCTACGCCCGCCCGCCGCATATTTTTCAGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGc                                                                                                      >  2:5164831/1‑150 (MQ=255)
                                            aCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGcc                                                                                                     <  2:645503/150‑1 (MQ=255)
                                              ggCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTg                                                                                                   >  2:4185745/1‑150 (MQ=255)
                                              ggCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTg                                                                                                   <  1:57567/150‑1 (MQ=255)
                                              ggCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTg                                                                                                   >  1:2377942/1‑150 (MQ=255)
                                                 tACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGAtt                                                                                                >  2:5220475/1‑150 (MQ=255)
                                                 tACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGAtt                                                                                                >  2:844377/1‑150 (MQ=255)
                                                        gcccgccGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAg                                                                                         >  2:3145577/1‑150 (MQ=255)
                                                          ccgccgCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGat                                                                                       <  1:135099/150‑1 (MQ=255)
                                                          ccgccgCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGat                                                                                       <  1:3945010/150‑1 (MQ=255)
                                                             ccgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATc                                                                                    >  2:3767300/1‑150 (MQ=255)
                                                              cgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATcc                                                                                   >  2:3521128/1‑150 (MQ=255)
                                                                   atTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGt                                                                              >  1:4509891/1‑150 (MQ=255)
                                                                        tCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGa                                                                         <  2:4692275/150‑1 (MQ=255)
                                                                          tGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGAtt                                                                       <  2:3076453/150‑1 (MQ=255)
                                                                          tGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGAtt                                                                       <  1:1195670/150‑1 (MQ=255)
                                                                            aagaaCGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTTAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTAc                                                                     <  2:3330144/150‑1 (MQ=255)
                                                                                  gCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGa                                                               <  1:1983478/150‑1 (MQ=255)
                                                                                   cTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGAATTGACTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGaa                                                              >  1:855461/1‑150 (MQ=255)
                                                                                     catcatTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAg                                                            >  1:2429189/1‑150 (MQ=255)
                                                                                     catcatTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAg                                                            >  1:1710970/1‑150 (MQ=255)
                                                                                       tcatTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGt                                                          >  1:4622168/1‑150 (MQ=255)
                                                                                        catTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGtt                                                         <  1:3370058/150‑1 (MQ=255)
                                                                                        catTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGtt                                                         >  1:1575393/1‑150 (MQ=255)
                                                                                        catTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGtt                                                         <  1:5174615/150‑1 (MQ=255)
                                                                                         atTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGttt                                                        <  2:2945283/150‑1 (MQ=255)
                                                                                          ttAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGtttt                                                       <  1:389288/150‑1 (MQ=255)
                                                                                          ttAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGtttt                                                       <  1:5136669/150‑1 (MQ=255)
                                                                                           tAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGttttt                                                      >  2:4051736/1‑150 (MQ=255)
                                                                                             gaagaaGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTa                                                    >  1:1551868/1‑150 (MQ=255)
                                                                                             gaagaaGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTa                                                    >  1:249059/1‑150 (MQ=255)
                                                                                              aagaagTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAg                                                   >  2:2801397/1‑150 (MQ=255)
                                                                                               agaagTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAgc                                                  >  2:2184204/1‑150 (MQ=255)
                                                                                                 aagTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAgcgc                                                <  1:2167679/150‑1 (MQ=255)
                                                                                                  agTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCt                                               <  1:4451999/150‑1 (MQ=255)
                                                                                                     aTTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGc                                            <  2:4587229/150‑1 (MQ=255)
                                                                                                     aTTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGc                                            <  1:2084073/150‑1 (MQ=255)
                                                                                                     aTTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGc                                            <  2:1142183/150‑1 (MQ=255)
                                                                                                           ggCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGatc                                      >  1:4688109/1‑150 (MQ=255)
                                                                                                           ggCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGatc                                      >  2:2617733/1‑150 (MQ=255)
                                                                                                            gCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGCGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGatca                                     <  1:1514595/150‑1 (MQ=255)
                                                                                                              gTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGatcatc                                   <  2:10248/150‑1 (MQ=255)
                                                                                                                 cccccACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCAt                                >  2:2425581/1‑150 (MQ=255)
                                                                                                                   cccACAGCGACCAGGCCGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTGTCGTAGTGCCTGATTGTGCCAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCCTGa                              >  2:4219438/1‑150 (MQ=255)
                                                                                                                    ccACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGAc                             <  2:3667127/150‑1 (MQ=255)
                                                                                                                      acaGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAg                           <  1:872065/150‑1 (MQ=255)
                                                                                                                        aGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTa                         >  2:3260815/1‑150 (MQ=255)
                                                                                                                          cGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTaaa                       >  2:1817983/1‑150 (MQ=255)
                                                                                                                            aCCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTaaaaa                     <  1:3404781/150‑1 (MQ=255)
                                                                                                                            aCCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTaaaaa                     >  2:1246001/1‑150 (MQ=255)
                                                                                                                               aGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACgg                  <  2:2025596/150‑1 (MQ=255)
                                                                                                                               aGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACgg                  >  1:2842681/1‑150 (MQ=255)
                                                                                                                                ggCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGa                 >  1:297714/1‑150 (MQ=255)
                                                                                                                                ggCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGa                 <  1:1833196/150‑1 (MQ=255)
                                                                                                                                   aGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGActc              >  2:347047/1‑150 (MQ=255)
                                                                                                                                      ggtgTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAatcatc           >  1:3545018/1‑150 (MQ=255)
                                                                                                                                         gtttCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGc        >  1:112834/4‑150 (MQ=255)
                                                                                                                                             ccATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCt    >  1:438028/1‑150 (MQ=255)
                                                                                                                                             ccATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCt    <  2:2752284/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CAGCTTGCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGC  >  CP009273/1799982‑1800274

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: