Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,014,981 G→T 7.8% R177S (CGT→AGT)  ompA ← outer membrane protein A (3a;II*;G;d)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,014,9810GT 7.8% 601.0 / 17.4 167R177S (CGT→AGT) ompAouter membrane protein A (3a;II*;G;d)
Reads supporting (aligned to +/- strand):  ref base G (81/72);  new base T (7/7);  total (88/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

TGAAGTTGAACAGAACGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTT  >  CP009273/1014834‑1015128
                                                                                                                                                   |                                                                                                                                                   
tgaagttgaaCAGAACGTCAGACTTCAGAGTGAAGTGATTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACg                                                                                                                                                     >  1:1937798/1‑148 (MQ=255)
 gaagttgaaCAGAACGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGa                                                                                                                                                    >  1:1900178/1‑148 (MQ=255)
 gaagttgaaCAGAACGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGa                                                                                                                                                    <  2:2065204/148‑1 (MQ=255)
 gaagttgaaCAGAACGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGa                                                                                                                                                    <  2:2065207/148‑1 (MQ=255)
    gttgaaCAGAACGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTg                                                                                                                                                 <  2:765119/148‑1 (MQ=255)
       gaaCAGAACGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCg                                                                                                                                              <  2:1599772/148‑1 (MQ=255)
       gaaCAGAACGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCg                                                                                                                                              <  1:2832822/148‑1 (MQ=255)
          cAGAACGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATg                                                                                                                                           >  1:1423546/1‑148 (MQ=255)
          cAGAACGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATg                                                                                                                                           >  1:399469/1‑148 (MQ=255)
          cAGAACGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATg                                                                                                                                           >  2:3372840/1‑148 (MQ=255)
           aGAACGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATgg                                                                                                                                          <  2:15181/148‑1 (MQ=255)
           aGAACGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATgg                                                                                                                                          <  2:407434/148‑1 (MQ=255)
             aaCGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGgtg                                                                                                                                        >  2:1928034/1‑148 (MQ=255)
              aCGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGgtgt                                                                                                                                       <  1:2907991/148‑1 (MQ=255)
              aCGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGgtgt                                                                                                                                       <  1:1249807/148‑1 (MQ=255)
              aCGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGgtgt                                                                                                                                       <  1:3210737/148‑1 (MQ=255)
               cGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACTAGTGCCGATGgtgtg                                                                                                                                      >  1:1359709/1‑148 (MQ=255)
               cGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACTAGTGCCGATGgtgtg                                                                                                                                      >  1:1359831/1‑148 (MQ=255)
                  cAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGc                                                                                                                                   >  1:3499013/1‑148 (MQ=255)
                  cAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGc                                                                                                                                   >  2:186358/1‑148 (MQ=255)
                  cAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGc                                                                                                                                   >  1:3498776/1‑148 (MQ=255)
                    gACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGt                                                                                                                                 >  1:2973343/1‑148 (MQ=255)
                      cTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCa                                                                                                                               >  2:662800/1‑148 (MQ=255)
                      cTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCa                                                                                                                               >  2:662761/1‑148 (MQ=255)
                      cTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCa                                                                                                                               <  1:1710996/148‑1 (MQ=255)
                       ttCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCAc                                                                                                                              <  2:2875765/148‑1 (MQ=255)
                         cAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCg                                                                                                                            <  2:2201149/148‑1 (MQ=255)
                           gagTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGAt                                                                                                                          >  2:625416/1‑148 (MQ=255)
                            agTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACTAGTGCCGATGGTGTGTGCGTCACCGAtg                                                                                                                         <  1:601522/148‑1 (MQ=255)
                            agTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACTAGTGCCGATGGTGTGTGCGTCACCGAtg                                                                                                                         <  1:2388347/148‑1 (MQ=255)
                            agTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGAtg                                                                                                                         >  2:3005687/1‑148 (MQ=255)
                            agTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGAtg                                                                                                                         >  1:3332349/1‑148 (MQ=255)
                            agTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGAtg                                                                                                                         >  2:192987/1‑148 (MQ=255)
                              tGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGAtgtt                                                                                                                       >  1:1155017/1‑148 (MQ=255)
                               gAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGAtgttg                                                                                                                      <  2:1452446/148‑1 (MQ=255)
                                   tGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTgg                                                                                                                  >  2:2723899/1‑148 (MQ=255)
                                   tGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTgg                                                                                                                  >  2:1987821/1‑148 (MQ=255)
                                    gCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGt                                                                                                                 <  2:2759352/148‑1 (MQ=255)
                                     cTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTc                                                                                                                >  1:3130037/1‑148 (MQ=255)
                                     cTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTc                                                                                                                >  1:737607/1‑148 (MQ=255)
                                          tCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTg                                                                                                           >  1:3249862/1‑148 (MQ=255)
                                          tCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTg                                                                                                           >  1:1376245/1‑148 (MQ=255)
                                              tACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTAt                                                                                                       >  1:2736756/1‑148 (MQ=255)
                                              tACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCCCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTAt                                                                                                       >  1:962750/1‑148 (MQ=255)
                                                       tGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACgg                                                                                              <  2:1021182/148‑1 (MQ=255)
                                                       tGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACgg                                                                                              <  2:1020159/148‑1 (MQ=255)
                                                         ccGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGt                                                                                            <  2:329938/148‑1 (MQ=255)
                                                             aGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCg                                                                                        >  2:487429/1‑148 (MQ=255)
                                                              gCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGa                                                                                       <  2:627531/148‑1 (MQ=255)
                                                                tGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGTGGGGTAGCGAtt                                                                                     <  2:2444725/148‑1 (MQ=255)
                                                                 ggAGCCGGAGCAACTACTGGCGCAGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTGCGGCCGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGAttt                                                                                    <  2:580653/148‑1 (MQ=255)
                                                                  gAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTc                                                                                   >  2:1498471/1‑148 (MQ=255)
                                                                  gAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTc                                                                                   <  1:2509040/148‑1 (MQ=255)
                                                                   aGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCa                                                                                  >  2:1187952/1‑148 (MQ=255)
                                                                        gAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGt                                                                             <  1:2688218/148‑1 (MQ=255)
                                                                         aGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTg                                                                            >  1:2094742/1‑148 (MQ=255)
                                                                             actactGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATcg                                                                        >  1:2705013/1‑148 (MQ=255)
                                                                             actactGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATcg                                                                        >  1:1703896/1‑148 (MQ=255)
                                                                                   ggAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACt                                                                  >  1:2775348/1‑148 (MQ=255)
                                                                                   ggAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACt                                                                  <  1:775292/148‑1 (MQ=255)
                                                                                   ggAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACt                                                                  <  2:653726/148‑1 (MQ=255)
                                                                                    gAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTc                                                                 >  1:73754/1‑148 (MQ=255)
                                                                                      gctgctTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCaa                                                               >  1:3198122/1‑148 (MQ=255)
                                                                                       ctgctTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAac                                                              >  2:3307418/1‑148 (MQ=255)
                                                                                        tgctTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAaca                                                             >  1:1503295/1‑148 (MQ=255)
                                                                                           ttCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCg                                                          <  2:3151741/148‑1 (MQ=255)
                                                                                           ttCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCg                                                          >  2:1229038/1‑148 (MQ=255)
                                                                                           ttCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCg                                                          <  2:1149021/148‑1 (MQ=255)
                                                                                                 cTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCg                                                    >  2:883447/1‑148 (MQ=255)
                                                                                                       gAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGAcc                                              >  2:2582817/1‑148 (MQ=255)
                                                                                                        aaaCGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCg                                             <  1:487429/148‑1 (MQ=255)
                                                                                                         aaCGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACTAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCgg                                            >  1:1497153/1‑148 (MQ=255)
                                                                                                         aaCGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCgg                                            >  1:191968/1‑148 (MQ=255)
                                                                                                          aCGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACTAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGga                                           >  2:1417224/1‑148 (MQ=255)
                                                                                                             gTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGaa                                        <  2:3217369/148‑1 (MQ=255)
                                                                                                                 gAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGc                                    <  1:3372840/148‑1 (MQ=255)
                                                                                                                   aaCACCCAGGCTCAGCATGCCGTTGTCCGGACTAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCg                                  >  1:1045723/1‑148 (MQ=255)
                                                                                                                      acCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGgtg                               >  1:1210734/1‑148 (MQ=255)
                                                                                                                        ccAGGCTCAGCATGCCGTTGTCCGGACTAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTc                             <  2:1359709/148‑1 (MQ=255)
                                                                                                                         cAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCg                            <  2:2740855/148‑1 (MQ=255)
                                                                                                                           ggCTCAGCATGCCGGTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTg                          <  2:3498776/148‑1 (MQ=255)
                                                                                                                            gCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTgg                         <  2:2982276/148‑1 (MQ=255)
                                                                                                                             cTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGt                        >  1:159741/1‑148 (MQ=255)
                                                                                                                              tCAGCATGCCGTTGTCCGGACTAGGGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGtt                       <  2:2451568/148‑1 (MQ=255)
                                                                                                                              tCAGCATGCCGTTGTCCGGACTAGGGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGtt                       <  2:2561056/148‑1 (MQ=255)
                                                                                                                               cAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGttt                      >  2:3360224/1‑148 (MQ=255)
                                                                                                                               cAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGttt                      >  2:348704/1‑148 (MQ=255)
                                                                                                                                    tGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTAcc                 <  2:3538342/148‑1 (MQ=255)
                                                                                                                                     gCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCa                >  1:694981/1‑148 (MQ=255)
                                                                                                                                     gCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCa                <  1:2072122/148‑1 (MQ=255)
                                                                                                                                     gCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCa                >  2:1282709/1‑148 (MQ=255)
                                                                                                                                      ccGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCAt               >  1:1012544/1‑148 (MQ=255)
                                                                                                                                        gTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATaa             <  2:3146797/148‑1 (MQ=255)
                                                                                                                                        gTTGGCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATaa             <  2:483905/148‑1 (MQ=255)
                                                                                                                                            tCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGt         <  2:1423546/148‑1 (MQ=255)
                                                                                                                                                gACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGa     >  2:3015994/1‑148 (MQ=255)
                                                                                                                                                gACGAGAGCCGTTGGTGTGTGCGTCACCGTTGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTTATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGa     <  1:70661/148‑1 (MQ=255)
                                                                                                                                                 aCGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGTTTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGAt    <  1:2284111/148‑1 (MQ=255)
                                                                                                                                                   gAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGAttt  <  1:2292723/148‑1 (MQ=255)
                                                                                                                                                   gAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGAttt  <  1:1094637/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGAAGTTGAACAGAACGTCAGACTTCAGAGTGAAGTGCTTGGTCTGTACTTCCGGTGCCGGAGCTGGAGCCGGAGCAACTACTGGAGCTGCTTCGCCCTGACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTT  >  CP009273/1014834‑1015128

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: