Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,221,550 G→T 23.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,5500GT23.3% 694.9 / 154.9 265pseudogene (539/1287 nt)mdtQpseudogene; putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base G (94/109);  new base T (26/36);  total (120/145)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 281 total aligned reads displayed.

GCAACGCGACCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGG  >  CP009273/2221401‑2221699
                                                                                                                                                     |                                                                                                                                                     
gCAACGCGACCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTg                                                                                                                                                       >  1:4108362/1‑150 (MQ=255)
    cgcgACCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTgaaga                                                                                                                                                   <  1:4235488/150‑1 (MQ=255)
     gcgACCGGTTTAAGCTTCAATGACGTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTgaagaa                                                                                                                                                  >  2:2622353/1‑150 (MQ=255)
      cgACCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGaagaag                                                                                                                                                 >  2:886484/1‑150 (MQ=255)
       gACCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGt                                                                                                                                                <  2:3258806/150‑1 (MQ=255)
        aCCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTa                                                                                                                                               >  1:403429/1‑150 (MQ=255)
          cGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAAt                                                                                                                                             >  2:1283903/1‑150 (MQ=255)
          cGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAAt                                                                                                                                             <  1:3608525/150‑1 (MQ=255)
          cGGTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAAt                                                                                                                                             <  2:4221980/150‑1 (MQ=255)
           ggTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATc                                                                                                                                            <  1:3506409/150‑1 (MQ=255)
            gTTTAAGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATcc                                                                                                                                           >  1:1475810/1‑150 (MQ=255)
                aaGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGtt                                                                                                                                       <  1:541312/150‑1 (MQ=255)
                 aGCTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTc                                                                                                                                      <  1:657429/150‑1 (MQ=255)
                   cTTCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTg                                                                                                                                    <  1:2983822/150‑1 (MQ=255)
                     tCAATGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGat                                                                                                                                  >  2:977776/1‑150 (MQ=255)
                       aaTGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGatat                                                                                                                                >  2:537470/1‑150 (MQ=255)
                       aaTGACTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGatat                                                                                                                                >  2:1523114/1‑150 (MQ=255)
                           aCTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCt                                                                                                                            <  2:4536974/150‑1 (MQ=255)
                            cTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTg                                                                                                                           <  1:246853/150‑1 (MQ=255)
                            cTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTg                                                                                                                           <  1:814134/150‑1 (MQ=255)
                            cTTTGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTg                                                                                                                           <  2:724362/150‑1 (MQ=255)
                              ttGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGc                                                                                                                         <  1:1823227/150‑1 (MQ=255)
                               tGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCg                                                                                                                        >  1:778138/1‑150 (MQ=255)
                               tGTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCg                                                                                                                        <  2:3977321/150‑1 (MQ=255)
                                gTCTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGa                                                                                                                       <  2:3209341/150‑1 (MQ=255)
                                  cTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATc                                                                                                                     <  2:554787/150‑1 (MQ=255)
                                  cTGGTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATc                                                                                                                     <  2:1535325/150‑1 (MQ=255)
                                    ggTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCgg                                                                                                                   >  1:3179249/1‑150 (MQ=255)
                                    ggTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCgg                                                                                                                   >  1:4465567/1‑150 (MQ=255)
                                     gTTATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGt                                                                                                                  <  2:3357326/150‑1 (MQ=255)
                                       tATTTGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTACTATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTcg                                                                                                                <  1:3124042/150‑1 (MQ=255)
                                           tGTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGAt                                                                                                            >  1:2948276/1‑150 (MQ=255)
                                            gTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATa                                                                                                           <  2:753156/150‑1 (MQ=255)
                                            gTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATa                                                                                                           <  2:753189/150‑1 (MQ=255)
                                            gTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATa                                                                                                           <  2:1000062/150‑1 (MQ=255)
                                            gTAAGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATa                                                                                                           <  1:1708563/150‑1 (MQ=255)
                                              aaGTGCGCTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAAt                                                                                                         <  2:392179/150‑1 (MQ=255)
                                                  gcgcTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGgtg                                                                                                     >  2:3755360/1‑150 (MQ=255)
                                                    gcTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGtt                                                                                                   <  2:2627590/150‑1 (MQ=255)
                                                    gcTTAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGtt                                                                                                   <  2:112789/150‑1 (MQ=255)
                                                      ttAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCt                                                                                                 <  2:4073144/150‑1 (MQ=255)
                                                      ttAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCt                                                                                                 <  2:225574/150‑1 (MQ=255)
                                                      ttAACCGTGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCt                                                                                                 <  1:2477572/150‑1 (MQ=255)
                                                             tGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCttt                                                                                          >  1:2593150/1‑150 (MQ=255)
                                                             tGCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCttt                                                                                          <  1:3701025/150‑1 (MQ=255)
                                                              gCCTCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCtttt                                                                                         >  1:1124649/1‑150 (MQ=255)
                                                                 tCAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCt                                                                                      <  1:3245822/150‑1 (MQ=255)
                                                                  cAATAATTTTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTa                                                                                     >  1:3805432/1‑150 (MQ=255)
                                                                      aatTTTCATTTTCCACGCGACATCGTTGATCTGCTGCCTGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTtg                                                                                 <  1:305407/150‑1 (MQ=255)
                                                                       attttcattttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTtgc                                                                                >  1:1979622/1‑150 (MQ=255)
                                                                        ttttcattttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTtgct                                                                               <  1:79965/150‑1 (MQ=255)
                                                                        ttttcattttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTtgct                                                                               <  2:1198449/150‑1 (MQ=255)
                                                                        ttttcattttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTtgct                                                                               >  1:4412278/1‑150 (MQ=255)
                                                                           tcattttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCa                                                                            <  2:1048843/150‑1 (MQ=255)
                                                                             attttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAg                                                                          >  1:3426774/1‑150 (MQ=255)
                                                                                   cccGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCgtg                                                                    <  2:4626960/150‑1 (MQ=255)
                                                                                   cccGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCgtg                                                                    >  1:4363815/1‑150 (MQ=255)
                                                                                     cgcgACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGtt                                                                  >  2:2551812/1‑150 (MQ=255)
                                                                                     cgcgACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGtt                                                                  <  1:456503/150‑1 (MQ=255)
                                                                                     cgcgACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGtt                                                                  <  1:2773809/150‑1 (MQ=255)
                                                                                      gcgACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGttt                                                                 <  2:1964649/150‑1 (MQ=255)
                                                                                      gcgACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGttt                                                                 <  2:2155743/150‑1 (MQ=255)
                                                                                      gcgACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGttt                                                                 >  1:1395030/1‑150 (MQ=255)
                                                                                      gcgACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGttt                                                                 <  1:2555988/150‑1 (MQ=255)
                                                                                        gACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTaa                                                               <  2:1997407/150‑1 (MQ=255)
                                                                                        gACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTaa                                                               >  2:2252088/1‑150 (MQ=255)
                                                                                        gACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTaa                                                               >  1:603796/1‑150 (MQ=255)
                                                                                         aCATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAAc                                                              >  2:341301/1‑150 (MQ=255)
                                                                                          cATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAAcg                                                             <  2:1412443/150‑1 (MQ=255)
                                                                                          cATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAAcg                                                             <  2:3489913/150‑1 (MQ=255)
                                                                                          cATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAAcg                                                             >  1:2088084/1‑150 (MQ=255)
                                                                                            tCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAAcgcc                                                           <  2:267620/150‑1 (MQ=255)
                                                                                             cGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAAcgccg                                                          >  2:4166382/1‑150 (MQ=255)
                                                                                              gTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAAcgccgc                                                         >  2:2753896/1‑150 (MQ=255)
                                                                                              gTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAAcgccgc                                                         >  1:434794/1‑150 (MQ=255)
                                                                                                      gctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGttt                                                 <  1:2057114/150‑1 (MQ=255)
                                                                                                      gctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGttt                                                 <  1:1056752/150‑1 (MQ=255)
                                                                                                      gctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGttt                                                 <  1:1042218/150‑1 (MQ=255)
                                                                                                      gctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGttt                                                 <  2:3071081/150‑1 (MQ=255)
                                                                                                        tgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGc                                               <  1:2414967/150‑1 (MQ=255)
                                                                                                           cGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCAc                                            <  2:832879/150‑1 (MQ=255)
                                                                                                             ggTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTc                                          >  1:3189899/1‑150 (MQ=255)
                                                                                                                   gCTGGCATTAATATCGGTTTCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTa                                    >  2:975751/1‑150 (MQ=255)
                                                                                                                   gCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTa                                    <  2:108773/150‑1 (MQ=255)
                                                                                                                     tGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACa                                  >  2:2516080/1‑150 (MQ=255)
                                                                                                                           tAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCggg                            <  2:484347/150‑1 (MQ=255)
                                                                                                                                tCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTGTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACg                       >  2:2329493/1‑150 (MQ=255)
                                                                                                                                  ggTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCt                     <  2:1755705/150‑1 (MQ=255)
                                                                                                                                    tttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTgc                   >  2:4414847/1‑150 (MQ=255)
                                                                                                                                      tCCACACCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTgcca                 <  2:3566973/150‑1 (MQ=255)
                                                                                                                                      tCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTgcca                 >  1:3907685/1‑150 (MQ=255)
                                                                                                                                      tCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTgcca                 <  2:663379/150‑1 (MQ=255)
                                                                                                                                       ccACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTgccag                >  1:4193512/1‑150 (MQ=255)
                                                                                                                                        cacaCCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTgccagc               >  1:3404212/1‑150 (MQ=255)
                                                                                                                                         acacCTTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTgccagcc              <  1:1683408/150‑1 (MQ=255)
                                                                                                                                            ccTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCcagccagc           >  1:2667875/1‑150 (MQ=255)
                                                                                                                                              ttCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCaa         >  2:4448618/1‑150 (MQ=255)
                                                                                                                                                     gaagaaGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCgg  >  1:2019048/1‑150 (MQ=255)
                                                                                                                                                     gaagaaGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCgg  <  1:1273216/150‑1 (MQ=255)
                                                                                                                                                     gaagaaGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCgg  >  1:2019086/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: