Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 325,161 C→G 9.8% 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 9.8% 252.2 / 16.2 82P149A (CCA→GCA) oprQouter‑membrane porin D
Reads supporting (aligned to +/- strand):  ref base C (35/39);  new base G (4/4);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

GTGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGG  >  NC_002947/325012‑325309
                                                                                                                                                      |                                                                                                                                                    
gtggCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGACAGACGGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGc                                                                                                                                                      <  1:1637558/150‑1 (MQ=255)
   gcggcgTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGAGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAg                                                                                                                                                   <  2:368082/150‑1 (MQ=255)
        gTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGGCAGtgctg                                                                                                                                              <  1:333841/150‑1 (MQ=255)
         tCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGGCAGtgctgg                                                                                                                                             <  2:1693061/150‑2 (MQ=255)
          cGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCa                                                                                                                                            <  1:1170908/150‑1 (MQ=255)
            aCTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGt                                                                                                                                          <  2:222458/150‑1 (MQ=255)
              ttcttcAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGGGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTac                                                                                                                                        >  1:1727150/1‑150 (MQ=255)
               tcttcAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTacg                                                                                                                                       <  2:909736/150‑1 (MQ=255)
                    aaGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGc                                                                                                                                                      <  2:843917/130‑1 (MQ=255)
                    aaGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGc                                                                                                                                                      >  1:843917/1‑130 (MQ=255)
                    aaGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACAACGTGCTGAAGTACGGTGACCAGATTCCAGCAGTGCCAGCGCTGCAGTCCGCCGAt                                                                                                                                  >  2:2230523/1‑150 (MQ=255)
                       ccGCACGAACCCACCAATGCCGAAGGCACCCCCCCCACGTGCATATACATCGCCTGTGGTTGTGCTGCTGTGTAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCa                                                                                                                               <  1:1115801/150‑1 (MQ=255)
                       ccGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCa                                                                                                                               <  2:85472/150‑1 (MQ=255)
                       ccGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCa                                                                                                                               >  2:1588288/1‑150 (MQ=255)
                            cGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGTGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTct                                                                                                                          <  2:717076/150‑1 (MQ=255)
                             gAAACCACCAATGCCGAAGGCAACGCCTGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTctg                                                                                                                         <  2:1267784/150‑1 (MQ=255)
                             gAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTctg                                                                                                                         >  2:1568392/1‑150 (MQ=255)
                                  caccaATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGGCAGtgctgg                                                                                                                                             >  2:1923800/1‑124 (MQ=255)
                                  caccaATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGGCAGtgctgg                                                                                                                                             <  1:1923800/125‑2 (MQ=255)
                                     caATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGa                                                                                                                 >  1:1111808/1‑150 (MQ=255)
                                      aaTGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGaa                                                                                                                >  2:1642960/1‑150 (MQ=255)
                                             aaGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCa                                                                                                         >  1:1671795/1‑150 (MQ=255)
                                              aGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCAc                                                                                                        >  1:2382108/1‑150 (MQ=255)
                                               ggCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCAcc                                                                                                       <  1:1535063/150‑1 (MQ=255)
                                                       cAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGCTGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCCCCGGCACCCt                                                                                               <  2:470475/150‑1 (MQ=255)
                                                          cTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGAt                                                                                            <  2:1821840/150‑1 (MQ=255)
                                                          cTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGAt                                                                                            <  2:1671795/150‑1 (MQ=255)
                                                           tCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGCACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATc                                                                                           <  1:58733/150‑1 (MQ=255)
                                                            cgccgcACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGGGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCa                                                                                          >  1:887814/1‑150 (MQ=255)
                                                            cgccgcACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAAAACACATCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCa                                                                                          >  1:887862/1‑150 (MQ=255)
                                                             gccgcACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCAc                                                                                         <  2:2382108/150‑1 (MQ=255)
                                                                    aaCCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCaaa                                                                                  <  1:1642960/150‑1 (MQ=255)
                                                                    aaCCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCAAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCaaa                                                                                  <  1:1233578/150‑1 (MQ=255)
                                                                       ccTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCAc                                                                                         <  1:936695/141‑1 (MQ=255)
                                                                       ccTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCAc                                                                                         >  2:936695/1‑141 (MQ=255)
                                                                             gcgTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCa                                                                          <  2:2361717/150‑1 (MQ=255)
                                                                                   tggtgCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGTCACCCTGATCACCTCCAAAGAGATCCAGGGCCt                                                                    <  1:318888/150‑1 (MQ=255)
                                                                                   tggtgCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCt                                                                    >  2:727871/1‑150 (MQ=255)
                                                                                      tgctgcTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTctgat                                                                                                                       <  2:2143856/96‑3 (MQ=255)
                                                                                      tgctgcTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTctgat                                                                                                                       >  1:2143856/1‑94 (MQ=255)
                                                                                       gctgctGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGaa                                                                <  1:1588288/150‑1 (MQ=255)
                                                                                            tgtgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGaa                                                           <  1:727871/150‑1 (MQ=255)
                                                                                            tgtgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGaa                                                           >  2:1662514/1‑150 (MQ=255)
                                                                                            tgtgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGaa                                                           >  2:1483359/1‑150 (MQ=255)
                                                                                             gtgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCa                                                                                                                                            <  2:2098054/68‑1 (MQ=255)
                                                                                             gtgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCa                                                                                                                                            >  1:2098054/1‑68 (MQ=255)
                                                                                             gtgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCCTCACCTGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAAc                                                          <  1:1134078/150‑1 (MQ=255)
                                                                                                   ttCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTctgct                                                                                                                       <  2:1854485/83‑1 (MQ=255)
                                                                                                   ttCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTctgct                                                                                                                       >  1:1854485/1‑83 (MQ=255)
                                                                                                      cGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATg                                                                                                                                 >  2:773306/1‑70 (MQ=255)
                                                                                                      cGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATg                                                                                                                                 <  1:773306/70‑1 (MQ=255)
                                                                                                      cGCATTGCCAACGCCGTGCTGAAGGACGGGGACCAGAAGCCAGCACTGCCAGTGCTGGAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGc                                                 <  2:1156902/150‑1 (MQ=255)
                                                                                                       gCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCt                                                >  1:2071200/1‑150 (MQ=255)
                                                                                                           ttCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGaa                                                                >  1:947200/1‑130 (MQ=255)
                                                                                                           ttCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGaa                                                                <  2:947200/130‑1 (MQ=255)
                                                                                                           ttCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCAc                                            >  1:1547363/1‑150 (MQ=255)
                                                                                                             ccAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCAccc                                          >  2:1435669/1‑150 (MQ=255)
                                                                                                              cAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGCTGCACGTCTGCTGCCAGAAAGCTTCACAGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCa                                         >  2:1435689/1‑150 (MQ=255)
                                                                                                              cAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGc                                                                      >  2:1792885/1‑121 (MQ=255)
                                                                                                              cAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGc                                                                      <  1:1792885/121‑1 (MQ=255)
                                                                                                               aaCACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAg                                        >  2:941040/1‑150 (MQ=255)
                                                                                                                acacCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCGTCACCCAgg                                       >  2:1635824/1‑150 (MQ=255)
                                                                                                                   ccGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAg                                    <  1:349603/150‑1 (MQ=255)
                                                                                                                    cGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAgc                                   >  1:2245609/1‑150 (MQ=255)
                                                                                                                      tGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAgcgg                                 <  1:1662514/150‑2 (MQ=255)
                                                                                                                       gCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTc                                                                                      <  2:758939/96‑1 (MQ=255)
                                                                                                                       gCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTc                                                                                      >  1:758939/1‑96 (MQ=255)
                                                                                                                        cTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAgcgcgc                               >  1:548341/1‑150 (MQ=255)
                                                                                                                        cTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGACACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAgcgcgc                               <  1:166105/150‑1 (MQ=255)
                                                                                                                                       cAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAgagcgccgagcgc                >  1:2036446/1‑150 (MQ=255)
                                                                                                                                        aGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAgagcgccgagcgcc               <  2:887814/150‑1 (MQ=255)
                                                                                                                                        aGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAgagcgccgagcgcc               <  2:887862/150‑1 (MQ=255)
                                                                                                                                         tactgcAGCACTGGCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCa                                         >  2:294711/6‑123 (MQ=255)
                                                                                                                                         tactgcAGCACTGGCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCa                                         <  1:294711/118‑1 (MQ=255)
                                                                                                                                          aTGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCGCCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGt             <  1:2163239/150‑1 (MQ=255)
                                                                                                                                            tgcAGCACTGGCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGa           >  1:374505/3‑150 (MQ=255)
                                                                                                                                            tgcAGCACTGGCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGa           >  1:374314/3‑150 (MQ=255)
                                                                                                                                             ccAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCAc                                                                                         <  2:512862/71‑1 (MQ=255)
                                                                                                                                             ccAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCAc                                                                                         >  1:512862/1‑71 (MQ=255)
                                                                                                                                              cAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAtg         >  1:2353863/1‑150 (MQ=255)
                                                                                                                                               aGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAtgg        >  1:1735477/1‑150 (MQ=255)
                                                                                                                                                   cTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAtggtggt    <  1:2230523/150‑1 (MQ=255)
                                                                                                                                                    tGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAtggtggtg   >  2:348682/1‑150 (MQ=255)
                                                                                                                                                     gCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAtggtggtgg  >  1:448445/1‑150 (MQ=255)
                                                                                                                                                     gCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCCCCGAGCGCCGTGAtggtggtgg  <  1:2265304/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GTGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGG  >  NC_002947/325012‑325309

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