Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,559,991 T→G 6.0% intergenic (+16/‑48) fepA → / → PP_2243 ferric enterobactin transport system outer membrane subunit/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,559,9910TG 6.0% 257.4 / 4.4 80intergenic (+16/‑48)fepA/PP_2243ferric enterobactin transport system outer membrane subunit/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (41/34);  new base G (2/3);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCGGGTAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCC  >  NC_002947/2559847‑2560138
                                                                                                                                                |                                                                                                                                                   
gTGCGGGTAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAACAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCgg                                                                                                                                                 <  4:304614/149‑1 (MQ=255)
    gggTAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGccc                                                                                                                                              >  3:378517/1‑148 (MQ=255)
      gTAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAAt                                                                                                                                           >  7:237681/1‑149 (MQ=255)
       tAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCAGGGCGGcccttgg                                                                                                                                          <  5:335352/149‑5 (MQ=255)
       tAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCAGGGCGGcccttgg                                                                                                                                          <  3:198394/149‑5 (MQ=255)
          gCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATcgcg                                                                                                                                       >  8:123179/1‑149 (MQ=255)
          gCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCAGGGCGGccctt                                                                                                                                            >  7:17461/1‑142 (MQ=255)
          gCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCAGGGCGGccctt                                                                                                                                            <  8:17461/144‑3 (MQ=255)
            gTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGt                                                                                                                                     >  7:171442/1‑149 (MQ=255)
               aacaacGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGt                                                                                                                                     <  5:266458/146‑1 (MQ=255)
               aacaacGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGt                                                                                                                                     >  6:266458/1‑146 (MQ=255)
                 caacGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAAt                                                                                                                                >  6:271380/1‑149 (MQ=255)
                 caacGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAAt                                                                                                                                >  2:93464/1‑149 (MQ=255)
                  aacGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATc                                                                                                                               <  5:249752/149‑1 (MQ=255)
                   acGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCg                                                                                                                              <  8:339276/149‑1 (MQ=255)
                     gTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCggg                                                                                                                            <  1:142590/149‑1 (MQ=255)
                      tCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGt                                                                                                                           <  4:378517/149‑1 (MQ=255)
                        ttCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTcg                                                                                                                         >  5:78536/1‑149 (MQ=255)
                          cGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTcgc                                                                                                                        <  8:237681/148‑1 (MQ=255)
                           gACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGa                                                                                                                      <  5:271380/149‑1 (MQ=255)
                              aaGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATg                                                                                                                   <  1:50507/149‑1 (MQ=255)
                                gCAGATCCTGCGCGGCGGTAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATcgc                                                                                                                                        <  5:246748/126‑1 (MQ=255)
                                gCAGATCCTGCGCGGCGGTAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATcgc                                                                                                                                        >  6:246748/1‑126 (MQ=255)
                                gCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGc                                                                                                                 <  7:346013/149‑1 (MQ=255)
                                gCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGc                                                                                                                 <  3:362421/149‑1 (MQ=255)
                                gCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGGGTTTGCTGCGCTGAAGATCACCTGCGCAGCCCACGGGCCCCCCGCCGGCCCAAGCGCGGTTCAAGCCGGACGCGAATgcc                                                                                                                 >  7:58756/1‑147 (MQ=255)
                                   gATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCcga                                                                                                              <  8:24618/149‑1 (MQ=255)
                                             ggcggcAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCtt                                                                                                    <  6:78536/149‑1 (MQ=255)
                                             ggcggcAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCtt                                                                                                    <  1:137053/149‑1 (MQ=255)
                                                ggcAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCa                                                                                                 <  7:123179/149‑1 (MQ=255)
                                                  cAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGc                                                                                               >  1:68405/1‑149 (MQ=255)
                                                   aaCGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGc                                                                                               >  6:39251/1‑148 (MQ=255)
                                                          gcagcTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGt                                                                                        >  7:51386/1‑148 (MQ=255)
                                                          gcagcTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGt                                                                                        <  8:51386/148‑1 (MQ=255)
                                                            agcTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTTAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGtctc                                                                                     <  8:171442/149‑1 (MQ=255)
                                                            agcTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGtctc                                                                                     <  4:172751/149‑1 (MQ=255)
                                                                     gCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGc                                                                            <  3:381817/149‑1 (MQ=255)
                                                                       aaaCACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCt                                                                           >  5:165753/1‑148 (MQ=255)
                                                                        aaCACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTa                                                                          >  8:341722/1‑148 (MQ=255)
                                                                         acacGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGc                                                                                <  7:377323/141‑1 (MQ=255)
                                                                         acacGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGc                                                                                >  8:377323/1‑141 (MQ=255)
                                                                              tataACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCggg                                                                                                                            <  5:4789/92‑1 (MQ=255)
                                                                              tataACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCggg                                                                                                                            >  6:4789/1‑92 (MQ=255)
                                                                              tataACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAACCGTATCGCAGCGTGCTAAACCagg                                                                   >  6:282928/1‑147 (MQ=255)
                                                                                   ccAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGATCTATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCAccc                                                              <  5:282928/149‑1 (MQ=255)
                                                                                    cAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGATTGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCAccc                                                              >  3:376178/1‑148 (MQ=255)
                                                                                    cAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCAccc                                                              >  3:10606/1‑148 (MQ=255)
                                                                                     aGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGAGCTAAACCATGCACCCCa                                                            >  6:282329/1‑149 (MQ=255)
                                                                                               ggCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGc                                                   >  2:257236/1‑148 (MQ=255)
                                                                                                gCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCga                                                 <  6:165753/149‑1 (MQ=255)
                                                                                                   gtgtTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCgaga                                               >  4:388598/1‑148 (MQ=255)
                                                                                                     gtTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCTGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAAt                                            >  1:333376/1‑149 (MQ=255)
                                                                                                          ctgcGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCt                                        <  2:68405/148‑1 (MQ=255)
                                                                                                           tgcGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTgg                                      >  5:160069/1‑149 (MQ=255)
                                                                                                             cgcTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGt                                     >  7:18762/1‑148 (MQ=255)
                                                                                                             cgcTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGt                                     <  6:160069/148‑1 (MQ=255)
                                                                                                               cTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGat                                  >  8:286881/1‑148 (MQ=255)
                                                                                                               cTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGa                                   >  8:34787/1‑148 (MQ=255)
                                                                                                                  aaTATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCga                                                 >  8:177853/1‑131 (MQ=255)
                                                                                                                  aaTATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCga                                                 <  7:177853/131‑1 (MQ=255)
                                                                                                                  aaTATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAg                                <  4:10606/148‑1 (MQ=255)
                                                                                                                   atatCAACTTCTGATCCCAAGGGCTGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGc                               <  2:100187/148‑1 (MQ=255)
                                                                                                                       cAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATcc                          >  3:130992/1‑149 (MQ=255)
                                                                                                                         aCTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCAGCTGCGAGAAATGCCTGGTGAAAGCCATCCGc                        >  2:69373/1‑149 (MQ=255)
                                                                                                                          cTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGcc                       <  3:388598/149‑1 (MQ=255)
                                                                                                                           ttCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCt                      >  4:378623/1‑149 (MQ=255)
                                                                                                                                aTCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCACGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGAc                   <  5:39251/147‑1 (MQ=255)
                                                                                                                                 tCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCt                >  7:55064/1‑149 (MQ=255)
                                                                                                                                  cccAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCt                >  7:382917/1‑148 (MQ=255)
                                                                                                                                  cccAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCt                >  5:309957/1‑148 (MQ=255)
                                                                                                                                   ccAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGa              >  6:174336/1‑149 (MQ=255)
                                                                                                                                   ccAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGa              >  2:242221/1‑149 (MQ=255)
                                                                                                                                    cAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGAt             >  3:312682/1‑149 (MQ=255)
                                                                                                                                     aaGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACAATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATc            >  8:220809/1‑149 (MQ=255)
                                                                                                                                         gccgccCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGc         <  4:62776/148‑1 (MQ=255)
                                                                                                                                         gccgccCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGc         >  3:62776/1‑148 (MQ=255)
                                                                                                                                             cccTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGc    >  4:195174/1‑149 (MQ=255)
                                                                                                                                             cccTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGc    <  5:174336/149‑1 (MQ=255)
                                                                                                                                              ccTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGcc   >  7:4046/1‑149 (MQ=255)
                                                                                                                                               cTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGccc  <  5:282329/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GTGCGGGTAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCCC  >  NC_002947/2559847‑2560138

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