Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,448,088 G→T 5.1% pseudogene (393/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,448,0880GT 5.1% 777.1 / 12.2 203pseudogene (393/2646 nt)yfcUpseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (114/78);  new base T (6/5);  total (120/83)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

CAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCATA  >  CP009273/2447943‑2448231
                                                                                                                                                 |                                                                                                                                               
 aGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCt                                                                                                                                              <  2:319323/148‑1 (MQ=255)
  gTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCtt                                                                                                                                             <  1:2642973/148‑1 (MQ=255)
   tGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTc                                                                                                                                            >  2:173317/1‑148 (MQ=255)
    gATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTcc                                                                                                                                           <  1:2709128/148‑1 (MQ=255)
           tAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGAc                                                                                                                                    >  2:279624/1‑148 (MQ=255)
              tCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCgg                                                                                                                                 >  1:3210811/1‑148 (MQ=255)
              tCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCgg                                                                                                                                 >  1:3137926/1‑148 (MQ=255)
                  cAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCtt                                                                                                                             >  2:2319493/1‑148 (MQ=255)
                   aaTGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTc                                                                                                                            <  1:1747652/148‑1 (MQ=255)
                     tGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAg                                                                                                                          >  1:84140/1‑148 (MQ=255)
                     tGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAg                                                                                                                          <  1:272987/148‑1 (MQ=255)
                     tGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAg                                                                                                                          >  1:2584462/1‑148 (MQ=255)
                      gATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAgg                                                                                                                         >  2:2828209/1‑148 (MQ=255)
                         cccGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCAt                                                                                                                      >  1:1665540/1‑148 (MQ=255)
                           cGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCAttt                                                                                                                    >  1:1667689/1‑148 (MQ=255)
                            ggAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTa                                                                                                                   <  2:2890895/148‑1 (MQ=255)
                             gagaTGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTAc                                                                                                                  <  1:553050/148‑1 (MQ=255)
                             gagaTGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTAc                                                                                                                  <  2:2528557/148‑1 (MQ=255)
                              agaTGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTAcc                                                                                                                 >  1:2996487/1‑148 (MQ=255)
                              agaTGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTAcc                                                                                                                 >  1:1528670/1‑148 (MQ=255)
                              agaTGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTAcc                                                                                                                 >  1:1484079/1‑148 (MQ=255)
                               gaTGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCa                                                                                                                >  1:2173377/1‑148 (MQ=255)
                               gaTGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCa                                                                                                                <  2:711557/148‑1 (MQ=255)
                                aTGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCag                                                                                                               >  1:2903004/1‑147 (MQ=255)
                                aTGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCAt                                                                                                               >  2:3163716/1‑148 (MQ=255)
                                     atcatcCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTgg                                                                                                          >  2:2427602/1‑148 (MQ=255)
                                         tcCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTcc                                                                                                      <  2:2859522/148‑1 (MQ=255)
                                          cccAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCa                                                                                                     >  1:2121208/1‑148 (MQ=255)
                                           ccAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCAt                                                                                                    <  2:1694538/148‑1 (MQ=255)
                                               cGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCa                                                                                                >  1:250965/1‑148 (MQ=255)
                                                gagaAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAg                                                                                               <  1:3123208/148‑1 (MQ=255)
                                                  gaAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGAt                                                                                             >  2:2380562/1‑148 (MQ=255)
                                                   aaGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAACCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGAtt                                                                                            <  2:2068265/148‑1 (MQ=255)
                                                     ggTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGAtttt                                                                                          >  2:804966/1‑148 (MQ=255)
                                                      gTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGAttttt                                                                                         >  2:888245/1‑148 (MQ=255)
                                                        ggATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACTCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCg                                                                                       >  2:39692/1‑148 (MQ=255)
                                                         gATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGc                                                                                      >  1:650343/1‑148 (MQ=255)
                                                         gATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGc                                                                                      >  1:2114300/1‑148 (MQ=255)
                                                         gATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGc                                                                                      >  1:2310279/1‑148 (MQ=255)
                                                          aTCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGcc                                                                                     >  2:2028719/1‑148 (MQ=255)
                                                           tCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCa                                                                                    >  1:1255080/1‑148 (MQ=255)
                                                           tCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACCAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCa                                                                                    >  1:669772/1‑148 (MQ=255)
                                                            cccAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCAc                                                                                   >  2:2165290/1‑148 (MQ=255)
                                                                gTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGtct                                                                               >  2:2350446/1‑148 (MQ=255)
                                                                gTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGtct                                                                               >  1:2439660/1‑148 (MQ=255)
                                                                 tCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGtctt                                                                              >  1:3159009/1‑148 (MQ=255)
                                                                 tCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGtctt                                                                              >  2:1863162/1‑148 (MQ=255)
                                                                       cAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTa                                                                        >  2:1771018/1‑148 (MQ=255)
                                                                       cAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTa                                                                        >  2:1771095/1‑148 (MQ=255)
                                                                        aGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTaa                                                                       <  1:3275195/148‑1 (MQ=255)
                                                                        aGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTaa                                                                       >  2:2646832/1‑148 (MQ=255)
                                                                         ggTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTaaa                                                                      >  2:2970517/1‑148 (MQ=255)
                                                                          gtagaggcGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAAc                                                                     <  1:11974/141‑1 (MQ=255)
                                                                           tataTTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAAcc                                                                    >  2:491724/1‑148 (MQ=255)
                                                                            atatTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCa                                                                   >  1:3168185/1‑148 (MQ=255)
                                                                             tatTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCaa                                                                  >  2:16175/1‑148 (MQ=255)
                                                                             tatTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCaa                                                                  <  1:2334282/148‑1 (MQ=255)
                                                                              atTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACTCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCaaa                                                                 >  2:1424749/1‑148 (MQ=255)
                                                                              atTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACTCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCaaa                                                                 >  2:616936/1‑148 (MQ=255)
                                                                              atTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCaaa                                                                 >  1:1372768/1‑148 (MQ=255)
                                                                                tCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACt                                                               >  2:1416401/1‑148 (MQ=255)
                                                                                    ggTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTgcgc                                                           <  1:279624/148‑1 (MQ=255)
                                                                                     gTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGcc                                                          >  1:1584272/1‑148 (MQ=255)
                                                                                       aaGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGccac                                                        <  2:701907/148‑1 (MQ=255)
                                                                                         gCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGccacca                                                      >  2:1183803/1‑148 (MQ=255)
                                                                                         gCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGccacca                                                      >  2:3124424/1‑148 (MQ=255)
                                                                                         gCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGccacca                                                      >  1:1859713/1‑148 (MQ=255)
                                                                                                gTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCgg                                               <  1:2044464/148‑1 (MQ=255)
                                                                                                   gTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGgtgt                                            >  1:1955277/1‑148 (MQ=255)
                                                                                                   gTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGgtgt                                            <  2:2772766/148‑1 (MQ=255)
                                                                                                      aaaTGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAg                                         <  2:2121208/148‑1 (MQ=255)
                                                                                                        aTGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGAc                                       >  1:2891865/1‑148 (MQ=255)
                                                                                                        aTGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGAc                                       >  1:808485/1‑148 (MQ=255)
                                                                                                         tGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACa                                      >  2:286867/1‑148 (MQ=255)
                                                                                                            cTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGc                                   >  1:1253833/1‑148 (MQ=255)
                                                                                                            cTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGc                                   >  1:355831/1‑148 (MQ=255)
                                                                                                              aaTGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCAt                                 <  2:552667/148‑1 (MQ=255)
                                                                                                                tGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAg                               >  1:996372/1‑148 (MQ=255)
                                                                                                                tGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAg                               >  1:1914356/1‑148 (MQ=255)
                                                                                                                  cGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCt                             <  2:2583475/148‑1 (MQ=255)
                                                                                                                    gACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCttt                           >  2:2668770/1‑148 (MQ=255)
                                                                                                                      cTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTg                         <  2:2899380/148‑1 (MQ=255)
                                                                                                                         gCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTa                      <  1:1668831/148‑1 (MQ=255)
                                                                                                                            tAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACg                   >  2:1231209/1‑148 (MQ=255)
                                                                                                                            tAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACg                   <  2:1255080/148‑1 (MQ=255)
                                                                                                                             aaaTCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGt                  >  1:913643/1‑148 (MQ=255)
                                                                                                                             aaaTCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGt                  <  1:2502397/148‑1 (MQ=255)
                                                                                                                              aaTCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTc                 <  2:2343351/148‑1 (MQ=255)
                                                                                                                              aaTCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTc                 <  2:3159009/148‑1 (MQ=255)
                                                                                                                              aaTCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTc                 <  2:1517718/148‑1 (MQ=255)
                                                                                                                                 cAGCCTTAATTTCCACTCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATc              <  1:2949140/148‑1 (MQ=255)
                                                                                                                                  aGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCt             <  2:2722337/148‑1 (MQ=255)
                                                                                                                                  aGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCt             <  1:2160632/148‑1 (MQ=255)
                                                                                                                                     cTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCa          <  2:1372768/148‑1 (MQ=255)
                                                                                                                                     cTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCa          >  1:935883/1‑148 (MQ=255)
                                                                                                                                       tAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCAcc        >  2:899312/1‑148 (MQ=255)
                                                                                                                                          tttCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGc     <  2:2722936/148‑1 (MQ=255)
                                                                                                                                            tCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCAt   <  1:2190/148‑1 (MQ=255)
                                                                                                                                             ccACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCATa  <  2:1667676/148‑1 (MQ=255)
                                                                                                                                             ccACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCATa  >  2:718942/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCATA  >  CP009273/2447943‑2448231

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: