Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 554,868 G→T 8.0% E147* (GAA→TAA)  sfmC → putative periplasmic pilus chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273554,8680GT 8.0% 692.1 / 23.5 194E147* (GAA→TAA) sfmCputative periplasmic pilus chaperone
Reads supporting (aligned to +/- strand):  ref base G (91/87);  new base T (7/9);  total (98/96)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 199 total aligned reads displayed.

CCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGCTGCACCCAGAG  >  CP009273/554721‑555015
                                                                                                                                                   |                                                                                                                                                   
ccGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAg                                                                                                                                                     >  1:2757895/1‑148 (MQ=255)
ccGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAg                                                                                                                                                     >  1:1761282/1‑148 (MQ=255)
ccGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCAGTCCAATTCAATTACAAg                                                                                                                                                     >  2:690756/1‑148 (MQ=255)
  gCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGaa                                                                                                                                                   <  2:1309405/148‑1 (MQ=255)
  gCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGaa                                                                                                                                                   <  1:877088/148‑1 (MQ=255)
   cTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAAt                                                                                                                                                  <  2:1400555/148‑1 (MQ=255)
    tGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAAtt                                                                                                                                                 <  1:3195516/148‑1 (MQ=255)
      gcagcaGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTAc                                                                                                                                               <  2:358866/148‑1 (MQ=255)
      gcagcaGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTAc                                                                                                                                               <  2:2218671/148‑1 (MQ=255)
         gcagATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCg                                                                                                                                            >  1:985084/1‑148 (MQ=255)
         gcagATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCg                                                                                                                                            >  2:3255803/1‑148 (MQ=255)
          cagATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGc                                                                                                                                           <  2:3303506/148‑1 (MQ=255)
          cagATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGc                                                                                                                                           >  2:2266375/1‑148 (MQ=255)
              tCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGaa                                                                                                                                       <  2:3035434/148‑1 (MQ=255)
              tCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGaa                                                                                                                                       <  1:1257719/148‑1 (MQ=255)
                 tGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAgcg                                                                                                                                    >  1:1485344/1‑148 (MQ=255)
                 tGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAgcg                                                                                                                                    >  1:2970111/1‑148 (MQ=255)
                  gAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAgcgc                                                                                                                                   >  1:944506/1‑148 (MQ=255)
                     tctcTGTTCTGTATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCgg                                                                                                                                >  1:265430/1‑148 (MQ=255)
                     tctcTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCgg                                                                                                                                >  2:3079894/1‑148 (MQ=255)
                       tctGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGac                                                                                                                              >  1:2982260/1‑148 (MQ=255)
                       tctGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGac                                                                                                                              >  1:2982215/1‑148 (MQ=255)
                        ctGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGaca                                                                                                                             >  1:3021525/1‑148 (MQ=255)
                            tCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAATAATTACCCGCAGAAGCGCCGGACACACt                                                                                                                         >  2:2608086/1‑148 (MQ=255)
                                gATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCATCTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAg                                                                                                                     >  2:2104935/1‑148 (MQ=255)
                                gATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAg                                                                                                                     >  1:2100849/1‑148 (MQ=255)
                                gATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAg                                                                                                                     <  1:1323573/148‑1 (MQ=255)
                                 aTGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGt                                                                                                                    <  2:2757222/148‑1 (MQ=255)
                                   gAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAATAATTACCCGCAGAAGCGCCGGACACACTCAAGttt                                                                                                                  <  1:2035904/148‑1 (MQ=255)
                                    aaTGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGtttt                                                                                                                 <  1:1316938/148‑1 (MQ=255)
                                     aTGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTc                                                                                                                >  1:2369422/1‑148 (MQ=255)
                                        ttAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTcgcg                                                                                                             <  2:2154834/148‑1 (MQ=255)
                                        ttAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTcgcg                                                                                                             >  1:1427650/1‑148 (MQ=255)
                                          aaGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGAt                                                                                                           <  2:453253/148‑1 (MQ=255)
                                            gACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATcc                                                                                                         >  1:1194884/1‑148 (MQ=255)
                                            gACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATcc                                                                                                         >  1:3074247/1‑148 (MQ=255)
                                             aCGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCg                                                                                                        >  1:331135/1‑148 (MQ=255)
                                              cGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCgg                                                                                                       <  1:2139617/148‑1 (MQ=255)
                                                 tCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGTTaa                                                                                                    >  1:2582854/1‑148 (MQ=255)
                                                 tCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTaa                                                                                                    >  1:2034539/1‑148 (MQ=255)
                                                   ccTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACt                                                                                                  <  1:3243681/148‑1 (MQ=255)
                                                     ttCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCACATGAAATTATTTCTCCGTCCAATTCAATTACAATAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACtat                                                                                                >  2:228471/1‑148 (MQ=255)
                                                      tCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCTCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACtata                                                                                               >  1:520217/1‑148 (MQ=255)
                                                       cGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTatat                                                                                              <  1:169658/148‑1 (MQ=255)
                                                           aGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAAt                                                                                          <  1:1025041/148‑1 (MQ=255)
                                                             aTAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGt                                                                                        <  1:2222155/148‑1 (MQ=255)
                                                               aaaaaTGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTc                                                                                      >  2:1981323/1‑148 (MQ=255)
                                                                aaaaTGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCa                                                                                     >  1:3151402/1‑148 (MQ=255)
                                                                 aaaTGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCAGGTAACTATATCAATGTTCAt                                                                                    >  1:2666912/1‑148 (MQ=255)
                                                                   aTGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATaa                                                                                  <  1:2165733/148‑1 (MQ=255)
                                                                        ttGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAATAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCAt                                                                             <  2:1319668/148‑1 (MQ=255)
                                                                         tGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATc                                                                            >  2:3122644/1‑148 (MQ=255)
                                                                         tGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATc                                                                            >  2:1665848/1‑148 (MQ=255)
                                                                           aaCGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCAc                                                                          <  1:2547306/148‑1 (MQ=255)
                                                                             cGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCt                                                                        >  1:615583/1‑148 (MQ=255)
                                                                              ggCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCtt                                                                       <  1:2713499/148‑1 (MQ=255)
                                                                              ggCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCtt                                                                       <  2:151414/148‑1 (MQ=255)
                                                                               gCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCttt                                                                      <  2:2099035/148‑1 (MQ=255)
                                                                                  ggAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTa                                                                   <  1:1234697/148‑1 (MQ=255)
                                                                                  ggAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTa                                                                   <  1:2884320/148‑1 (MQ=255)
                                                                                   gAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTAt                                                                  >  1:1903347/1‑148 (MQ=255)
                                                                                     aTGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGt                                                                <  2:441151/148‑1 (MQ=255)
                                                                                        ttttGCGACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCAc                                                             >  2:2871821/1‑148 (MQ=255)
                                                                                                 tGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTaa                                                    <  2:3026343/148‑1 (MQ=255)
                                                                                                 tGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTaa                                                    <  1:338531/148‑1 (MQ=255)
                                                                                                   gCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCNGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACt                                                  <  2:2467296/148‑1 (MQ=255)
                                                                                                   gCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACt                                                  <  2:2625650/148‑1 (MQ=255)
                                                                                                    cGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACtt                                                 >  1:383058/1‑148 (MQ=255)
                                                                                                       ttttATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACa                                              <  1:1167514/148‑1 (MQ=255)
                                                                                                       ttttATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACa                                              >  2:3205221/1‑148 (MQ=255)
                                                                                                       ttttATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACa                                              >  1:558847/1‑148 (MQ=255)
                                                                                                       ttttATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACa                                              <  1:2837632/148‑1 (MQ=255)
                                                                                                        tttATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACaa                                             <  1:2871847/148‑1 (MQ=255)
                                                                                                          tATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAATAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGt                                           <  2:2658332/148‑1 (MQ=255)
                                                                                                           aTCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTg                                          <  2:1936136/148‑1 (MQ=255)
                                                                                                           aTCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTg                                          <  2:265430/148‑1 (MQ=255)
                                                                                                           aTCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTg                                          >  1:1961439/1‑148 (MQ=255)
                                                                                                             cgcgCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTGTTTATGTCACCCTGGTTAACTTACAAGTggg                                        >  1:2116308/1‑148 (MQ=255)
                                                                                                                    gAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCaa                                 >  2:2878163/1‑148 (MQ=255)
                                                                                                                    gAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCaa                                 >  2:793718/1‑148 (MQ=255)
                                                                                                                     aaaTTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCaaa                                >  1:1827644/1‑148 (MQ=255)
                                                                                                                          attTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTg                           >  2:3193239/1‑148 (MQ=255)
                                                                                                                                ccGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAAt                     >  1:1005170/1‑148 (MQ=255)
                                                                                                                                ccGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAAt                     >  1:1005112/1‑148 (MQ=255)
                                                                                                                                    ccAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTa                 <  2:1485344/148‑1 (MQ=255)
                                                                                                                                    ccAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTa                 >  2:104364/1‑148 (MQ=255)
                                                                                                                                    ccAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTa                 <  1:4680/148‑1 (MQ=255)
                                                                                                                                     caattcaattACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTAt                <  2:572293/148‑1 (MQ=255)
                                                                                                                                      aattcaattACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATg               >  2:3136007/1‑148 (MQ=255)
                                                                                                                                        ttcaattACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGc             <  2:2976514/148‑1 (MQ=255)
                                                                                                                                        ttcaattACAAGAATAACCCGCAGAAGCGCCGGACACACTAAAGATATCGCGATCCGCTAACTATATAAAAGTTCATAATCCATCACATTTTTATCTCAACCTGGTTAACTTACAAGTGGGAAACCAAAAGAAGTTGAATGATATgtc             >  1:146184/1‑146 (MQ=255)
                                                                                                                                            attACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGCTGCa         >  2:2972678/1‑148 (MQ=255)
                                                                                                                                            attACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGCTGCa         >  1:847297/1‑148 (MQ=255)
                                                                                                                                            attACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGCTGCa         <  2:170939/148‑1 (MQ=255)
                                                                                                                                             ttACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGCTGCAc        >  2:625165/1‑148 (MQ=255)
                                                                                                                                              tACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGCTGCAcc       >  1:963121/1‑148 (MQ=255)
                                                                                                                                              tACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGCTGCAcc       >  1:1263924/1‑148 (MQ=255)
                                                                                                                                              tACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGCTGCAcc       >  2:904136/1‑148 (MQ=255)
                                                                                                                                                cAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGCTGCACCCa     >  1:1162736/1‑148 (MQ=255)
                                                                                                                                                  aGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGCTGCACCCaga   <  2:1371459/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGCTGCACCCAGAG  >  CP009273/554721‑555015

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: