Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 325,161 C→G 8.5% P149A (CCA→GCA)  oprQ → outer‑membrane porin D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947325,1610CG 8.5% 247.1 / 12.5 82P149A (CCA→GCA) oprQouter‑membrane porin D
Reads supporting (aligned to +/- strand):  ref base C (34/41);  new base G (2/5);  total (36/46)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

TGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGG  >  NC_002947/325013‑325309
                                                                                                                                                    |                                                                                                                                                    
tggCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGcc                                                                                                                                                     <  2:2028233/150‑1 (MQ=255)
 ggcggcGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCa                                                                                                                                                    <  1:200859/150‑1 (MQ=255)
   cggcgTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGt                                                                                                                                                  <  1:1028215/150‑1 (MQ=255)
    ggcgTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCa                                                                                                                                                    >  2:1823890/1‑147 (MQ=255)
    ggcgTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCa                                                                                                                                                    <  1:1823890/147‑1 (MQ=255)
     gcgTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCACATTTCAAACACCGTGCTGAAGTACGGTGACCAGAGGCCAGCACTGCCAGtgc                                                                                                                                                >  2:314681/1‑150 (MQ=255)
       gTCGACTTCTTCAAGCCGCACGAAGCCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGtgctg                                                                                                                                              <  1:1395313/150‑1 (MQ=255)
       gTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGGCAGtgctg                                                                                                                                              <  2:231405/150‑1 (MQ=255)
        tCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGc                                                                                                                                                      >  2:1295967/1‑141 (MQ=255)
        tCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGc                                                                                                                                                      <  1:1295967/141‑1 (MQ=255)
        tCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGGCAGtgctgg                                                                                                                                             <  2:1490148/150‑2 (MQ=255)
        tCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGtgctgc                                                                                                                                             >  2:651309/1‑150 (MQ=255)
         cGACTTCTTCAAGCCAGACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCa                                                                                                                                            <  1:1442580/150‑1 (MQ=255)
          gACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAg                                                                                                                                           >  1:1139561/1‑150 (MQ=255)
           aCTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGGCAGtgctggcat                                                                                                                                          <  1:1057707/150‑5 (MQ=255)
            cttcttCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATG‑CAGCACTGGCAGtgctggcatct                                                                                                                                        <  1:719986/150‑7 (MQ=255)
             ttcttcAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCTCTGCCAGTGCTGCAGTac                                                                                                                                        <  2:757766/150‑1 (MQ=255)
              tcttcAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCTAGCACTGGCAGtgctggcatctg                                                                                                                                       <  2:939337/150‑8 (MQ=255)
                ttcAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTacgac                                                                                                                                     >  1:1165926/1‑150 (MQ=255)
                   aaGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGAt                                                                                                                                  <  2:2245811/150‑1 (MQ=255)
                     gCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGc                                                                                                                                <  2:732426/150‑1 (MQ=255)
                         cACGAAAACCACAATGCCGAAGCCGAGGCCAGCCAGCTGCACAACCTGGAGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAATTACGGTGACCAGATGCCAGCCCTGCCAGTGCTGCAGTACGACGATGCAcaa                                                                                                                            >  2:594008/1‑148 (MQ=255)
                          aCGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTc                                                                                                                           >  1:689762/1‑150 (MQ=255)
                           cGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTct                                                                                                                          <  2:2375928/150‑1 (MQ=255)
                           cGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTct                                                                                                                          <  1:2200877/150‑1 (MQ=255)
                           cGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTct                                                                                                                          <  1:2387337/150‑1 (MQ=255)
                                ccaccaATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCtgctgc                                                                                                                     <  2:1139561/150‑1 (MQ=255)
                                  accaATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCa                                                                                                                   >  1:1127617/1‑150 (MQ=255)
                                  accaATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCa                                                                                                                   >  1:208401/1‑150 (MQ=255)
                                    caATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGa                                                                                                                 >  1:1760433/1‑150 (MQ=255)
                                    caATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGa                                                                                                                 >  1:2180870/1‑150 (MQ=255)
                                      aTGCCGAAGTCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGaaa                                                                                                               <  1:103627/150‑1 (MQ=255)
                                      aTGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCACTACGACGATGCACGTCTGCTGCCAGaaa                                                                                                               >  2:985618/1‑150 (MQ=255)
                                       ggCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAg                                                                                                              <  2:533629/149‑1 (MQ=255)
                                        gCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGc                                                                                                             <  1:637949/150‑1 (MQ=255)
                                        gCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGc                                                                                                             <  2:1334721/150‑1 (MQ=255)
                                         ccGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCt                                                                                                            >  1:2021037/1‑150 (MQ=255)
                                             aGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGtgc                                                                                                                                                <  2:230314/110‑1 (MQ=255)
                                             aGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGtgc                                                                                                                                                >  1:230314/1‑110 (MQ=255)
                                               gCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCg                                                                                                      <  2:1760433/150‑1 (MQ=255)
                                                  aCGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCa                                                                                                   >  1:289383/1‑150 (MQ=255)
                                                            gccgcACAACCTGGCGCGTGGTGGTGCTGCTGTGACGATGCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCAc                                                                                         <  2:803457/150‑1 (MQ=255)
                                                               gcACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTc                                                                                      <  1:269659/150‑1 (MQ=255)
                                                                   aaCCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCaaa                                                                                  <  2:208401/150‑1 (MQ=255)
                                                                   aaCCTGGCGAGCCGTGGTGCTTCATCGAACCTCCGCATTTCCGACTCCGTTTTGAATTACTGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCAccaagggc                                                                                  <  1:1289577/150‑7 (MQ=255)
                                                                         gcgcGTGGTGTTGCTGCTGTGAAGTTCCGCTTTTCCTACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATc                                                                            <  2:749924/150‑1 (MQ=255)
                                                                           gcgTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCGGAAAGCTTCACCGGCACCCTGATCACCTCCaaa                                                                                  >  2:1650053/1‑142 (MQ=255)
                                                                           gcgTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCaaa                                                                                  <  1:1650053/142‑1 (MQ=255)
                                                                           gcgTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCa                                                                          <  2:1284517/150‑1 (MQ=255)
                                                                                       tgctgTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACt                                                              <  1:651309/150‑1 (MQ=255)
                                                                                         ctgTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGATGCCTTATAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGa                                                            >  1:24821/1‑150 (MQ=255)
                                                                                          tgtgAAGTTCCGCCTTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGaa                                                           >  2:2153848/1‑150 (MQ=255)
                                                                                          tgtgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGaa                                                           <  2:623520/150‑1 (MQ=255)
                                                                                                 ttCCGCATTTCCAACACCGTACTGAAGGACGGTGACTAGATTCCAGCACTGCCATTGCTGCAGTACGACGATGAACGTCTACTGCCAGAAAGCCTCACCGCCAAGCCGACGACCTGCAAATATAACCAGTGCCTGGAACTGAACgccggc                                                    >  2:1130879/1‑150 (MQ=255)
                                                                                                  tCCGCATTTCCAACACCCTGCTGAATTACGGTGACCAGATGCCAGCACGGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACgccggcc                                                   <  1:262999/150‑1 (MQ=255)
                                                                                                    cGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGc                                                                      <  2:1552047/129‑1 (MQ=255)
                                                                                                    cGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGc                                                                      >  1:1552047/1‑129 (MQ=255)
                                                                                                    cGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGAACGTATCACACGTTTGAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGc                                                 >  1:1354442/1‑150 (MQ=255)
                                                                                                     gCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCtt                                                                                                           >  2:2041745/1‑91 (MQ=255)
                                                                                                     gCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCtt                                                                                                           <  1:2041745/91‑1 (MQ=255)
                                                                                                     gCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCt                                                >  1:382304/1‑150 (MQ=255)
                                                                                                       aTTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTc                                              <  1:1416342/150‑1 (MQ=255)
                                                                                                          tCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCAcc                                           <  2:1165926/150‑1 (MQ=255)
                                                                                                             aaCACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGAGGATGCACGTCTGCTGACAGAAAGCTTCACCGGCACCCGGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCGGCTTCACCCAg                                        >  2:297075/1‑150 (MQ=255)
                                                                                                             aaCACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAg                                        <  2:1885630/150‑1 (MQ=255)
                                                                                                                    tGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAgcgc                                 <  2:1127617/150‑1 (MQ=255)
                                                                                                                        gAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCaa                             <  2:382304/150‑1 (MQ=255)
                                                                                                                         aaGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAag                            >  1:2281820/1‑150 (MQ=255)
                                                                                                                              cGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACCACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAgagcgc                       <  1:1831159/150‑1 (MQ=255)
                                                                                                                                   aCCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTc                                                                                                          <  2:1517466/62‑1 (MQ=255)
                                                                                                                                   aCCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTc                                                                                                          >  1:1517466/1‑62 (MQ=255)
                                                                                                                                    ccAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAgcgc                                 >  1:1646271/1‑134 (MQ=255)
                                                                                                                                    ccAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAgcgc                                 <  2:1646271/134‑1 (MQ=255)
                                                                                                                                     cAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAgagcgccgagcgc                >  2:1256369/1‑150 (MQ=255)
                                                                                                                                     cAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAgagcgccgagcgc                >  1:957183/1‑150 (MQ=255)
                                                                                                                                        aTGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGt             >  1:429420/1‑150 (MQ=255)
                                                                                                                                          tgcAGCACTGGCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGa           >  1:533808/3‑150 (MQ=255)
                                                                                                                                          tgcAGCACTGGCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGa           >  1:850355/3‑150 (MQ=255)
                                                                                                                                          gCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGa           >  1:2250821/1‑150 (MQ=255)
                                                                                                                                            cAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAtg         <  2:689762/150‑1 (MQ=255)
                                                                                                                                            cAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAtg         >  1:897853/1‑150 (MQ=255)
                                                                                                                                             aGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAtgg        <  2:289383/150‑1 (MQ=255)
                                                                                                                                              ccACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCCCCAAAGAGCTCCAGGGCCTGGAACGGAACGCCGGCCGCGTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAtggt       <  1:297075/149‑1 (MQ=255)
                                                                                                                                               cACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAtggtg      >  2:1637726/1‑150 (MQ=255)
                                                                                                                                                   gCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAtggtggtgg  >  1:357070/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGG  >  NC_002947/325013‑325309

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑