Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,608,994 G→C 6.7% R202G (CGC→GGC)  waaA ← 3‑deoxy‑D‑manno‑octulosonic acid transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,608,9940GC 6.7% 195.0 / 4.2 60R202G (CGC→GGC) waaA3‑deoxy‑D‑manno‑octulosonic acid transferase
Reads supporting (aligned to +/- strand):  ref base G (22/34);  new base C (2/2);  total (24/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

AAGGCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCC  >  NC_002947/5608849‑5609142
                                                                                                                                                 |                                                                                                                                                    
aaGGCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGACGATGGTCTGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCt                                                                                                                                                    <  1:277028/148‑1 (MQ=255)
 aGGCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCt                                                                                                                                                    >  3:264427/1‑147 (MQ=255)
   gCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCGGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAgg                                                                                                                                                <  3:152849/149‑1 (MQ=255)
     aTCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGc                                                                                                                                              <  5:345600/149‑1 (MQ=255)
          cACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGTCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGc                                                                                                                                          <  5:350658/148‑1 (MQ=255)
           aCCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCCCTGAGCGCGTGCaa                                                                                                                                        <  5:344444/149‑2 (MQ=255)
             cgtcgtGAGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCAccc                                                                                                                                      <  3:189636/149‑1 (MQ=255)
               tcgtGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTgacctggg                                                                                                                                    <  6:196158/149‑8 (MQ=255)
               tcgtGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAg                                                                                                                                    <  7:3050/149‑1 (MQ=255)
                    ggTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCAc                                                                                                                               >  3:156262/1‑149 (MQ=255)
                         tGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAg                                                                                                                          >  1:370656/1‑149 (MQ=255)
                           gcggcAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAgcg                                                                                                                        >  5:303933/1‑149 (MQ=255)
                                aaTCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGTCGGTGACTTGCACGCCCTGAGCGCGTGCaa                                                                                                                                        <  6:240842/128‑2 (MQ=255)
                                aaTCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGTCGGTGACTTGCACGCCCTGAGCGCGTGCaa                                                                                                                                        >  5:240842/1‑127 (MQ=255)
                                aaTCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACACAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTggg                                                                                                                   >  1:105807/1‑149 (MQ=255)
                                  tCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGcc                                                                                                                 >  2:371230/1‑149 (MQ=255)
                                    cacacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTc                                                                                                               <  4:264427/149‑1 (MQ=255)
                                     acacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCg                                                                                                              <  8:162650/149‑1 (MQ=255)
                                     acacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGGCGATGGGCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCg                                                                                                              <  5:16672/149‑1 (MQ=255)
                                      cacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTgg                                                                                                                    <  3:307973/142‑1 (MQ=255)
                                      cacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTgg                                                                                                                    >  4:307973/1‑142 (MQ=255)
                                       acaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGt                                                                                                            <  3:149739/149‑1 (MQ=255)
                                         aGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGttt                                                                                                          <  4:129257/149‑1 (MQ=255)
                                          gggCGCTGGCCGCCCCCCAACTGTTCACACAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTc                                                                                                         <  2:105807/149‑1 (MQ=255)
                                           ggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                        <  8:55520/149‑1 (MQ=255)
                                           ggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                        <  3:19717/149‑1 (MQ=255)
                                           ggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGGCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                        <  2:192825/149‑1 (MQ=255)
                                            gcgcTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCgg                                                                                                       <  4:156262/149‑1 (MQ=255)
                                             cgcTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGt                                                                                                      <  7:116353/149‑1 (MQ=255)
                                             cgcTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCAGTTTCGGt                                                                                                      <  2:370656/149‑1 (MQ=255)
                                               cTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCt                                                                                                    <  6:303933/149‑1 (MQ=255)
                                                     gccCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACgg                                                                                              <  2:72806/149‑1 (MQ=255)
                                                         cccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTc                                                                                                     >  1:84932/1‑138 (MQ=255)
                                                         cccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTc                                                                                                     <  2:84932/138‑1 (MQ=255)
                                                          ccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCAc                                                                                                >  1:53136/1‑142 (MQ=255)
                                                          ccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCAc                                                                                                <  2:53136/142‑1 (MQ=255)
                                                           cAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCa                                                                                        <  1:160178/149‑1 (MQ=255)
                                                                       cAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGAGCAGGCTCATCTCgg                                                                            <  5:330966/149‑1 (MQ=255)
                                                                        aaTGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCgg                                                                            >  4:265762/1‑148 (MQ=255)
                                                                          tGCCCTCGCGCGTGGCAGCAGCTGTGCGTCGATGGGCAGGTCGAACTTGATCGCGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTTGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCa                                                                         <  6:129552/149‑1 (MQ=255)
                                                                              cTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCAt                                                                     >  5:281647/1‑149 (MQ=255)
                                                                                 gcgcgTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCgg                                                                  >  8:49370/1‑149 (MQ=255)
                                                                                     gTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGTGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCgg                                                              <  2:78184/149‑1 (MQ=255)
                                                                                     gTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCgg                                                              <  1:99462/149‑1 (MQ=255)
                                                                                     gTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCgg                                                              <  1:371230/149‑1 (MQ=255)
                                                                                       ggCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGt                                                            >  1:36518/1‑149 (MQ=255)
                                                                                              gcTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCAATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTg                                                     <  2:270511/149‑1 (MQ=255)
                                                                                               cTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTgg                                                    >  7:219401/1‑149 (MQ=255)
                                                                                                                   cGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACgg                                <  7:31547/149‑1 (MQ=255)
                                                                                                                       cTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGcc                             >  3:306278/1‑148 (MQ=255)
                                                                                                                       cTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCg                            >  2:223874/1‑149 (MQ=255)
                                                                                                                            tCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACCGGCCGAgcg                        >  1:126230/1‑148 (MQ=255)
                                                                                                                              gAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCt                      >  7:208238/1‑148 (MQ=255)
                                                                                                                                gCCGGTGACTTTCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCt                      >  2:245600/1‑146 (MQ=255)
                                                                                                                                gCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCa                    <  5:322010/148‑1 (MQ=255)
                                                                                                                                         ttGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAgcgcg          >  6:62410/1‑149 (MQ=255)
                                                                                                                                          tGCACGCCCTGAGCGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGCCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAgcgcg          >  7:271525/1‑148 (MQ=255)
                                                                                                                                             aCGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCa                                                          >  2:71765/1‑97 (MQ=255)
                                                                                                                                             aCGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCa                                                          <  1:71765/97‑1 (MQ=255)
                                                                                                                                             aCGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCAtt      >  2:110656/1‑149 (MQ=255)
                                                                                                                                              cgcgCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTg     <  2:178581/149‑1 (MQ=255)
                                                                                                                                                 gcTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGcc  <  6:281647/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AAGGCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCC  >  NC_002947/5608849‑5609142

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: ‑