Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,446,543 G→A 5.0% pseudogene (1938/2646 nt) yfcU ← pseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,446,5430GA 5.0% 1081.1 / 28.3 297pseudogene (1938/2646 nt)yfcUpseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (140/142);  new base A (7/8);  total (147/150)
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.91e-01

Only 100 of 297 total aligned reads displayed.

CGGCAACGCCATCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAG  >  CP009273/2446395‑2446690
                                                                                                                                                    |                                                                                                                                                   
 ggCAACGCCATCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTc                                                                                                                                                   <  1:5041098/150‑1 (MQ=255)
        ccATCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTgc                                                                                                                                            >  1:2413115/1‑150 (MQ=255)
        ccATCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTgc                                                                                                                                            >  1:2413315/1‑150 (MQ=255)
          aTCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTCGTTCGCACTGAGGTCAACCTgcgc                                                                                                                                          >  2:4526876/1‑150 (MQ=255)
          aTCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTgcgc                                                                                                                                          >  2:3247585/1‑150 (MQ=255)
                atcaatcaGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGa                                                                                                                                    >  2:3406231/1‑150 (MQ=255)
                   aatcaGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAAcc                                                                                                                                 >  2:1528087/1‑150 (MQ=255)
                            gCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCt                                                                                                                        <  2:4740408/150‑1 (MQ=255)
                             cGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGATCAACCTGCGCCAGCGAACCATCATGGCTg                                                                                                                       <  1:225139/150‑1 (MQ=255)
                               tgtCCCGCCCATATTCTGGTTACGGTGAAGTGCGCCGCCGTGGTGAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTgta                                                                                                                     <  2:2688201/150‑1 (MQ=255)
                               tgtCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTgta                                                                                                                     <  1:2836146/150‑1 (MQ=255)
                                gtCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTgtag                                                                                                                    <  2:1203712/150‑1 (MQ=255)
                                 tCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTgtagt                                                                                                                   <  2:2159006/150‑1 (MQ=255)
                                  cccGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTgtagta                                                                                                                  <  1:2987397/150‑1 (MQ=255)
                                    cGCCCATATTCTGGGTACGGTGGAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGc                                                                                                                >  1:686964/1‑150 (MQ=255)
                                       ccATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCAt                                                                                                             <  1:4308821/150‑1 (MQ=255)
                                       ccATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCAt                                                                                                             <  1:4640766/150‑1 (MQ=255)
                                         atatTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCa                                                                                                           <  2:480606/150‑1 (MQ=255)
                                         atatTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCa                                                                                                           >  1:78840/1‑150 (MQ=255)
                                            ttCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACg                                                                                                        <  1:1673463/150‑1 (MQ=255)
                                            ttCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACg                                                                                                        <  1:4877834/150‑1 (MQ=255)
                                              cTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCt                                                                                                      >  2:374011/1‑150 (MQ=255)
                                               tGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTg                                                                                                     >  2:363511/1‑150 (MQ=255)
                                                  gTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGgtg                                                                                                  >  2:3284888/1‑150 (MQ=255)
                                                     cGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGgtgtgt                                                                                               <  1:2701471/150‑1 (MQ=255)
                                                     cGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGgtgtgt                                                                                               >  1:4419752/1‑150 (MQ=255)
                                                     cGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGgtgtgt                                                                                               <  2:580520/150‑1 (MQ=255)
                                                     cGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGATCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGgtgtgt                                                                                               >  2:3491926/1‑150 (MQ=255)
                                                      ggTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGtt                                                                                              >  1:846331/1‑150 (MQ=255)
                                                      ggTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGtt                                                                                              <  1:4322782/150‑1 (MQ=255)
                                                       gTGAAGTGCGCCTCCGTGGGGAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGATCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGttt                                                                                             <  2:2834448/150‑1 (MQ=255)
                                                           aGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCa                                                                                         >  1:3873299/1‑150 (MQ=255)
                                                            gTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCAc                                                                                        <  1:266011/150‑1 (MQ=255)
                                                             tGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACt                                                                                       <  1:4729053/150‑1 (MQ=255)
                                                              gcgcCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTg                                                                                      <  1:3496343/150‑1 (MQ=255)
                                                               cgccgcCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTgg                                                                                     >  1:3933869/1‑150 (MQ=255)
                                                               cgccgcCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTgg                                                                                     <  1:2940299/150‑1 (MQ=255)
                                                                gccgccGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGt                                                                                    >  1:5143004/1‑150 (MQ=255)
                                                                 ccgccgTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTg                                                                                   >  2:4284941/1‑150 (MQ=255)
                                                                  cgccgTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGc                                                                                  >  2:4698017/1‑150 (MQ=255)
                                                                   gccgTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGcc                                                                                 <  1:275791/150‑1 (MQ=255)
                                                                    ccgTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCg                                                                                <  2:1272618/150‑1 (MQ=255)
                                                                    ccgTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCg                                                                                <  2:4386435/150‑1 (MQ=255)
                                                                    ccgTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCg                                                                                <  2:2927436/150‑1 (MQ=255)
                                                                     cgTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGa                                                                               <  1:1317109/150‑1 (MQ=255)
                                                                     cgTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGa                                                                               >  1:1010969/1‑150 (MQ=255)
                                                                             gtaAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCACATCATGGTATTTCGCACTGAGGACAACCATCGCCAGCGAACCATCATGCCAGTAGTAGCCATCAACGCTGGTGTGTTAGTCACTGGTGCCGATGTTCAAc                                                                       >  1:897119/1‑150 (MQ=255)
                                                                              taAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACt                                                                      >  2:800880/1‑150 (MQ=255)
                                                                               aaGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGATGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTg                                                                     <  2:1953823/150‑1 (MQ=255)
                                                                               aaGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTg                                                                     <  1:4997593/150‑1 (MQ=255)
                                                                               aaGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTg                                                                     >  2:1873979/1‑150 (MQ=255)
                                                                                 gCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGAt                                                                   <  1:1583066/150‑1 (MQ=255)
                                                                                 gCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGAt                                                                   <  2:2786551/150‑1 (MQ=255)
                                                                                 gCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGATCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGAt                                                                   <  1:3431647/150‑1 (MQ=255)
                                                                                  cGTTGCGCCGCCCTGTAACGACACGCCCGCAGAGGTGTAGTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATa                                                                  <  2:4106041/150‑1 (MQ=255)
                                                                                  cGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATa                                                                  >  2:4629071/1‑150 (MQ=255)
                                                                                  cGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATa                                                                  <  1:1806651/150‑1 (MQ=255)
                                                                                        gccgcaCTGTAACGACAAGCCCGCAGTGGTGTATTGCCCTTCATGGTTGTTCGCACTGAGGTCAACCTGCGCCTGCGAACCATCATGGCTGTAGTAGCCATGAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGt                                                            <  1:2635579/150‑1 (MQ=255)
                                                                                          cgccCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTcg                                                          >  1:5038096/1‑150 (MQ=255)
                                                                                            cccTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTcgcg                                                        >  1:3175836/1‑150 (MQ=255)
                                                                                             ccTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGt                                                       >  1:4201744/1‑150 (MQ=255)
                                                                                                 tAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATc                                                   >  1:495325/1‑150 (MQ=255)
                                                                                                  aaCGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCg                                                  <  2:1442031/150‑1 (MQ=255)
                                                                                                   aCGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATCGTGAGTCGCGTCATCGa                                                 >  2:3422650/1‑150 (MQ=255)
                                                                                                     gACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGaca                                               >  2:4910704/1‑150 (MQ=255)
                                                                                                        aaGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACgg                                            >  2:2380661/1‑150 (MQ=255)
                                                                                                         aGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGc                                           <  1:2476627/150‑1 (MQ=255)
                                                                                                           cccGCAGAGGTGTATGGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTa                                         <  1:4817528/150‑1 (MQ=255)
                                                                                                            ccGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTaa                                        <  2:5152156/150‑1 (MQ=255)
                                                                                                            ccGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGATCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTaa                                        >  2:2436631/1‑150 (MQ=255)
                                                                                                             cGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTaaa                                       >  2:3526143/1‑150 (MQ=255)
                                                                                                             cGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTaaa                                       >  2:901039/1‑150 (MQ=255)
                                                                                                             cGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTaaa                                       >  2:2591005/1‑150 (MQ=255)
                                                                                                                agagGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATa                                    <  1:19257/150‑1 (MQ=255)
                                                                                                                  agGGGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATTTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAAc                                  >  2:1326883/1‑150 (MQ=255)
                                                                                                                   ggTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAAcc                                 <  1:2151908/150‑1 (MQ=255)
                                                                                                                    gtgtATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCg                                <  2:3987886/150‑1 (MQ=255)
                                                                                                                    gtgtATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCg                                >  1:560226/1‑150 (MQ=255)
                                                                                                                    gtgtATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCg                                >  2:3079270/1‑150 (MQ=255)
                                                                                                                       tATTGCCCTTCATGGTAGTTCGCACTGAGATCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGAcc                             >  2:5078246/1‑150 (MQ=255)
                                                                                                                        aTTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCt                            <  1:1346202/150‑1 (MQ=255)
                                                                                                                        aTTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCt                            >  1:1499710/1‑150 (MQ=255)
                                                                                                                            cccTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGgct                        >  1:1739117/1‑150 (MQ=255)
                                                                                                                             ccTTCATGGTAGTTCGCACTGAGGTCAACCTGCGACAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGgctg                       <  2:2308164/150‑1 (MQ=255)
                                                                                                                               ttCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGgctgct                     <  2:4467807/150‑1 (MQ=255)
                                                                                                                               ttCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGgctgct                     <  1:1147579/150‑1 (MQ=255)
                                                                                                                               ttCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGgctgct                     <  1:1410358/150‑1 (MQ=255)
                                                                                                                                   tGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTgtcc                 >  2:4228621/1‑150 (MQ=255)
                                                                                                                                    ggTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTgtccg                <  1:695114/150‑1 (MQ=255)
                                                                                                                                     gTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTgtccgt               <  2:2489832/150‑1 (MQ=255)
                                                                                                                                     gTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTgtccgt               <  2:5167183/150‑1 (MQ=255)
                                                                                                                                        gTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTccc            >  2:1659161/1‑150 (MQ=255)
                                                                                                                                        gTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTccc            <  2:5216755/150‑1 (MQ=255)
                                                                                                                                          tCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCAc          <  2:1740453/150‑1 (MQ=255)
                                                                                                                                          tCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCAc          >  2:2989802/1‑150 (MQ=255)
                                                                                                                                           cGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACt         >  2:2876853/1‑150 (MQ=255)
                                                                                                                                                tGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCAt    >  1:2299593/1‑150 (MQ=255)
                                                                                                                                                 gAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTCAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATa   >  1:3094946/1‑150 (MQ=255)
                                                                                                                                                  aGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAg  >  1:148434/1‑150 (MQ=255)
                                                                                                                                                  aGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAg  >  1:1381253/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGCAACGCCATCGGCATCAATCAGCAGGCGTGTCCCGCCCATATTCTGGGTACGGTGAAGTGCGCCGCCGTGGGTAGTAAGCGTTGCGCCGCCCTGTAACGACAAGCCCGCAGAGGTGTATTGCCCTTCATGGTAGTTCGCACTGAGGTCAACCTGCGCCAGCGAACCATCATGGCTGTAGTAGCCATCAACGCTGGTGTGTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAG  >  CP009273/2446395‑2446690

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: