Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,221,550 G→T 9.3% pseudogene (539/1287 nt) mdtQ ← pseudogene; putative channel/filament proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,221,5500GT 9.3% 1045.6 / 51.8 280pseudogene (539/1287 nt)mdtQpseudogene; putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base G (122/132);  new base T (6/20);  total (128/152)
Fisher's exact test for biased strand distribution p-value = 2.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 291 total aligned reads displayed.

GCAACGCGACCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGG  >  CP009273/2221401‑2221699
                                                                                                                                                     |                                                                                                                                                     
gCAACGCGACCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTg                                                                                                                                                       <  1:4424687/150‑1 (MQ=255)
gCAACGCGACCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTg                                                                                                                                                       >  1:4582055/1‑150 (MQ=255)
     gcgACCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTtaagaa                                                                                                                                                  <  2:141801/150‑1 (MQ=255)
      cgACCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGaagaag                                                                                                                                                 <  2:1935553/150‑1 (MQ=255)
      cgACCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGt                                                                                                                                                >  1:4204296/1‑150 (MQ=255)
       gACCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGt                                                                                                                                                >  2:2370805/1‑150 (MQ=255)
       gACCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGt                                                                                                                                                <  1:2618626/150‑1 (MQ=255)
       gACCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGt                                                                                                                                                >  2:1308256/1‑150 (MQ=255)
        aCCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTa                                                                                                                                               <  1:1931675/150‑1 (MQ=255)
          cGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAAt                                                                                                                                             <  2:2111728/150‑1 (MQ=255)
            gTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATcc                                                                                                                                           <  1:2207059/150‑1 (MQ=255)
              ttAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCg                                                                                                                                         <  1:2932646/150‑1 (MQ=255)
                aaGGTTCAATGTCCTTGTATTGTTATTTGTTAGTGGGCTTTACCGGGCCACAATAATTTTCAATTTCCCCGCCTCATCGGTTCGGCGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGtt                                                                                                                                       <  2:2509192/150‑1 (MQ=255)
                    ttCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGa                                                                                                                                   <  1:4097026/150‑1 (MQ=255)
                      cAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGata                                                                                                                                 >  1:212831/1‑150 (MQ=255)
                       aaTGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGatat                                                                                                                                <  1:2357699/150‑1 (MQ=255)
                       aaTGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGatat                                                                                                                                >  1:4377187/1‑150 (MQ=255)
                          gACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGc                                                                                                                             <  2:3370166/150‑1 (MQ=255)
                            cTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTg                                                                                                                           <  2:3003861/150‑1 (MQ=255)
                                 tCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGAt                                                                                                                      >  1:326621/1‑150 (MQ=255)
                                  cTGGTTATTTGTCAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATc                                                                                                                     <  2:2929188/150‑1 (MQ=255)
                                  cTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATc                                                                                                                     <  2:471248/150‑1 (MQ=255)
                                       tATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTcg                                                                                                                >  2:4374891/1‑150 (MQ=255)
                                         tttGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTcgcg                                                                                                              <  1:3549607/150‑1 (MQ=255)
                                          ttGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGa                                                                                                             >  2:2432935/1‑150 (MQ=255)
                                           tGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGAt                                                                                                            <  1:2483498/150‑1 (MQ=255)
                                            gTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATa                                                                                                           <  2:4145006/150‑1 (MQ=255)
                                               aGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATg                                                                                                        <  2:3683129/150‑1 (MQ=255)
                                                gTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATgg                                                                                                       <  2:4455026/150‑1 (MQ=255)
                                                   cgcTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGgtgt                                                                                                    <  1:3196620/150‑1 (MQ=255)
                                                    gcTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGtt                                                                                                   <  2:493241/150‑1 (MQ=255)
                                                    gcTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGtt                                                                                                   <  2:3521039/150‑1 (MQ=255)
                                                    gcTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGtt                                                                                                   <  2:2603708/150‑1 (MQ=255)
                                                     cTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTc                                                                                                  <  1:2227858/150‑1 (MQ=255)
                                                      ttAACCGTGTCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCt                                                                                                 <  1:2383000/150‑1 (MQ=255)
                                                      ttAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCt                                                                                                 <  2:3389364/150‑1 (MQ=255)
                                                       tAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTg                                                                                                >  1:995162/1‑150 (MQ=255)
                                                       tAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTg                                                                                                <  1:728807/150‑1 (MQ=255)
                                                        aaCCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGc                                                                                               <  1:4508744/150‑1 (MQ=255)
                                                        aaCCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGc                                                                                               >  2:4181250/1‑150 (MQ=255)
                                                           cGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGctct                                                                                            >  2:1236285/1‑150 (MQ=255)
                                                             tGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCttt                                                                                          <  2:1664157/150‑1 (MQ=255)
                                                              gCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCtttt                                                                                         <  1:3082948/150‑1 (MQ=255)
                                                               ccTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTTATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCttttt                                                                                        <  1:568505/150‑1 (MQ=255)
                                                                  cAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTa                                                                                     <  1:2541904/150‑1 (MQ=255)
                                                                        ttttcattttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTtgct                                                                               <  2:4318830/150‑1 (MQ=255)
                                                                        ttttcattttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTtgct                                                                               <  2:4318766/150‑1 (MQ=255)
                                                                          ttcattttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTtgctgc                                                                             <  1:4348359/150‑1 (MQ=255)
                                                                          ttcattttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTtgctgc                                                                             >  2:3957201/1‑150 (MQ=255)
                                                                            cattttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCaa                                                                           <  1:3627418/150‑1 (MQ=255)
                                                                             attttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAg                                                                          >  2:3559547/1‑150 (MQ=255)
                                                                              gtttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGa                                                                         <  1:1192950/149‑1 (MQ=255)
                                                                                ttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGAcc                                                                       <  1:371581/150‑1 (MQ=255)
                                                                                ttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGAcc                                                                       <  2:713506/150‑1 (MQ=255)
                                                                                 tcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCg                                                                      >  1:1076597/1‑150 (MQ=255)
                                                                                 tcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCg                                                                      <  2:888915/150‑1 (MQ=255)
                                                                                  ccccGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCgt                                                                     <  1:1319419/150‑1 (MQ=255)
                                                                                      gcgACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGttt                                                                 <  2:1106027/150‑1 (MQ=255)
                                                                                        gACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTaa                                                               >  2:2580983/1‑150 (MQ=255)
                                                                                        gACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTaa                                                               >  1:207112/1‑150 (MQ=255)
                                                                                         aCATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAAc                                                              >  1:413139/1‑150 (MQ=255)
                                                                                          cATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAAcg                                                             >  1:1116520/1‑150 (MQ=255)
                                                                                           aTCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAAcgc                                                            >  2:674094/1‑150 (MQ=255)
                                                                                             cGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAAcgccg                                                          <  1:697417/150‑1 (MQ=255)
                                                                                             cGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAAcgccg                                                          <  1:2850454/150‑1 (MQ=255)
                                                                                                tGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCt                                                       >  2:4301720/1‑150 (MQ=255)
                                                                                                tGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCt                                                       >  2:1358861/1‑150 (MQ=255)
                                                                                                  aGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTgg                                                     >  2:1215084/1‑150 (MQ=255)
                                                                                                     tgctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGtt                                                  >  2:3247688/1‑150 (MQ=255)
                                                                                                       ctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTg                                                <  2:200722/150‑1 (MQ=255)
                                                                                                         gcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGcc                                              >  1:3447487/1‑150 (MQ=255)
                                                                                                            gggTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCCTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGAACGTGTTTAACGCCGCCTGGGTTTGCCACt                                           >  2:4596247/1‑150 (MQ=255)
                                                                                                              gTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAGGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTcc                                         >  1:2681044/1‑150 (MQ=255)
                                                                                                              gTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTcc                                         >  1:3898368/1‑150 (MQ=255)
                                                                                                              gTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTcc                                         >  1:3523492/1‑150 (MQ=255)
                                                                                                              gTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTcc                                         >  2:2639177/1‑150 (MQ=255)
                                                                                                               ttttGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTccc                                        >  2:4424093/1‑150 (MQ=255)
                                                                                                                 ttGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAg                                      >  1:3915851/1‑150 (MQ=255)
                                                                                                                 ttGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAg                                      >  1:1498635/1‑150 (MQ=255)
                                                                                                                     tGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACa                                  >  1:4033445/1‑150 (MQ=255)
                                                                                                                      ggCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTGCAg                                 <  2:251378/150‑1 (MQ=255)
                                                                                                                           tAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCggg                            <  2:1604842/150‑1 (MQ=255)
                                                                                                                             atatCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCt                          >  2:1541457/1‑150 (MQ=255)
                                                                                                                             atatCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCt                          >  2:55864/1‑150 (MQ=255)
                                                                                                                             atatCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCt                          <  2:3275252/150‑1 (MQ=255)
                                                                                                                              tatCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTa                         <  1:1728429/150‑1 (MQ=255)
                                                                                                                              tatCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTa                         <  1:1412253/150‑1 (MQ=255)
                                                                                                                              tatCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTa                         <  1:81230/150‑1 (MQ=255)
                                                                                                                                  ggTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCt                     <  1:766009/150‑1 (MQ=255)
                                                                                                                                   gTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCGACGCTg                    >  1:511005/1‑150 (MQ=255)
                                                                                                                                     ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTgcc                  >  1:2588145/1‑150 (MQ=255)
                                                                                                                                     ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTgcc                  >  2:2688666/1‑150 (MQ=255)
                                                                                                                                       ccACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTgccag                <  1:3770690/150‑1 (MQ=255)
                                                                                                                                       ccACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTgccag                <  1:3770709/150‑1 (MQ=255)
                                                                                                                                        cacaCCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACTCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTgccagc               >  2:2672412/1‑150 (MQ=255)
                                                                                                                                        cacaCCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTgccagc               >  1:1564150/1‑150 (MQ=255)
                                                                                                                                         acacCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTgccagcc              >  2:3837505/1‑150 (MQ=255)
                                                                                                                                           acCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGccagccag            >  1:2653222/1‑150 (MQ=255)
                                                                                                                                                    tGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCg   >  2:2242941/1‑150 (MQ=255)
                                                                                                                                                     gaagaaGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCgg  >  2:2490473/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAACGCGACCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGG  >  CP009273/2221401‑2221699

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: