Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,960,967 G→C 7.2% intergenic (‑149/+115) fliE ← / ← atoC flagellar hook‑basal body complex protein FliE/two‑component system DNA‑binding transcriptional activator AtoC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,960,9670GC 7.2% 175.0 / 2.0 62intergenic (‑149/+115)fliE/atoCflagellar hook‑basal body complex protein FliE/two‑component system DNA‑binding transcriptional activator AtoC
Reads supporting (aligned to +/- strand):  ref base G (24/33);  new base C (3/2);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

CCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGAC  >  NC_002947/4960821‑4961105
                                                                                                                                                   |                                                                                                                                          
cccTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGAATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCaga                                                                                                                                           <  1:265845/149‑1 (MQ=255)
cccTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGTCGCagag                                                                                                                                          >  8:82281/1‑149 (MQ=255)
cccTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCagag                                                                                                                                          >  5:262335/1‑149 (MQ=255)
 ccTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGc                                                                                                                                         >  3:13207/1‑149 (MQ=255)
        tctcCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCg                                                                                                                                  >  8:229031/1‑149 (MQ=255)
        tctcCGGCGGCCCCATCTTTAGGACAGTCATGCGGCGAATGACCTACCTAAACCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGccc                                                                                                                                    <  6:252714/147‑1 (MQ=39)
          tcCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGACCGGCCCCGGc                                                                                                                                >  5:77377/1‑149 (MQ=255)
            cTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCNTCGCGATAGGGCCGCAGAGCGGCCCCGGCg                                                                                                                               <  1:33813/148‑1 (MQ=255)
            cTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTCTGCGGccctatc                                                                                                                                <  3:150722/147‑5 (MQ=255)
                ggCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCtt                                                                                                                          <  6:175994/149‑1 (MQ=255)
                  ccGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTTGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTg                                                                                                                        <  1:53463/149‑1 (MQ=255)
                  ccGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTg                                                                                                                        <  7:229031/149‑1 (MQ=255)
                      aTTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTCTGCGGcccta                                                                                                                                  >  8:167586/1‑133 (MQ=255)
                      aTTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTCTGCGGcccta                                                                                                                                  <  7:167586/135‑3 (MQ=255)
                       tttttttGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgcggcg                                                                                                                   <  1:79529/149‑1 (MQ=255)
                       tttttttGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgcggcg                                                                                                                   >  7:80075/1‑149 (MQ=255)
                       tttttttGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgcggcg                                                                                                                   >  1:243439/1‑149 (MQ=255)
                       tttttttGACAGTCATGCGGCGAATGACCTACCAAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgcggcg                                                                                                                   <  2:45918/149‑1 (MQ=255)
                          ttttGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGCTCTTTGCGGCGGAg                                                                                                                >  6:73183/1‑149 (MQ=255)
                              gACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTg                                                                                                             <  5:86077/148‑1 (MQ=255)
                                  gTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGAc                                                                                                        >  4:313583/1‑149 (MQ=255)
                                      tGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTgg                                                                                                    <  3:250254/149‑1 (MQ=255)
                                        cggcgAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGGCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTggg                                                                                                   <  7:82281/148‑1 (MQ=255)
                                            gAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGc                                                                                              <  5:122857/149‑1 (MQ=255)
                                                 acctacctAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTcgc                                                                                         <  3:262032/149‑1 (MQ=255)
                                                 acctacctAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTcg                                                                                          <  3:313583/148‑1 (MQ=255)
                                                 acctacctAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGAACTTTGCGGCGGTGGTGAAGACATGGGGCAGCTAcgc                                                                                         >  4:262032/1‑149 (MQ=38)
                                                  cctacctaAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTcgcg                                                                                        >  7:91255/1‑149 (MQ=255)
                                                  cctacctaAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTcgcg                                                                                        >  6:303058/1‑149 (MQ=255)
                                                   ctacctaAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCgcgc                                                                                       >  5:369974/1‑149 (MQ=255)
                                                        taAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCAt                                                                                  >  1:189816/1‑149 (MQ=255)
                                                         aaaGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGCCCGGGGGCCGCGTAGCGCCCCAtt                                                                                 >  2:9599/1‑149 (MQ=255)
                                                         aaaGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCGTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCAt                                                                                  >  7:302636/1‑148 (MQ=255)
                                                             caacaaAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGAGGTTTGCGGGGGCGCTGGAGACCTGGGGCCGCTTGGCGgccccgtggcg                                                                             >  8:303970/1‑142 (MQ=38)
                                                               acaaAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGa                                                                            >  4:303758/1‑148 (MQ=255)
                                                               acaaAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCATGCGTCGCGATAGGGCCGCACAGCGGCCCCGGCGAGCATTGCGGCGGAGCTGAAGACCTGGGGCCCCTTAGCGACCCATTcgacac                                                                           >  4:220312/1‑145 (MQ=38)
                                                                 aaaGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGGGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGacac                                                                         <  8:91255/149‑1 (MQ=255)
                                                                       gCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGcc                                                                   <  5:258016/149‑1 (MQ=255)
                                                                            ccTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTcc                                                              <  3:220312/149‑1 (MQ=255)
                                                                              ttCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCca                                                            >  8:246415/1‑149 (MQ=255)
                                                                                 gCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACag                                                         <  2:189816/149‑1 (MQ=255)
                                                                                  cAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACaga                                                        <  7:267961/149‑1 (MQ=255)
                                                                                         tttCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCcgc                                                 <  5:73183/149‑1 (MQ=255)
                                                                                              cAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAg                                            <  7:246415/149‑1 (MQ=255)
                                                                                               agagATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCag                                                                                  <  6:37388/110‑2 (MQ=255)
                                                                                               agagATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCag                                                                                  >  5:37388/1‑109 (MQ=255)
                                                                                                    aTCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCa                                       <  8:260383/149‑1 (MQ=255)
                                                                                                      cacacaATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACaga                                                        <  1:157793/130‑1 (MQ=255)
                                                                                                        cacaATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACaga                                                        >  7:230258/1‑128 (MQ=255)
                                                                                                        cacaATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACaga                                                        <  8:230258/128‑1 (MQ=255)
                                                                                                            aTTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGAttt                               <  3:303758/149‑1 (MQ=255)
                                                                                                            aTTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGAttt                               <  2:243439/149‑1 (MQ=255)
                                                                                                               ccTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGc                             <  5:303058/148‑1 (MQ=255)
                                                                                                                   cGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCCCGCCAGTTGCAGATGATTTGCCAGTc                        <  1:9599/149‑1 (MQ=255)
                                                                                                                      gAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGc                     >  7:35330/1‑149 (MQ=255)
                                                                                                                            ccTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGaa                 >  5:342672/1‑147 (MQ=255)
                                                                                                                              ttGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGAttt                               <  7:328749/131‑1 (MQ=255)
                                                                                                                              ttGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGAttt                               >  8:328749/1‑131 (MQ=255)
                                                                                                                              ttGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGt                         >  1:164445/1‑137 (MQ=255)
                                                                                                                                 cGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCt           <  8:35330/148‑1 (MQ=255)
                                                                                                                                 cGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTg          <  6:369974/149‑1 (MQ=255)
                                                                                                                                   tCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCgcgc                                               <  8:137548/110‑1 (MQ=255)
                                                                                                                                   tCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCgcgc                                               >  7:137548/1‑110 (MQ=255)
                                                                                                                                      cgATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTc     <  6:342672/149‑1 (MQ=255)
                                                                                                                                         cggggCCGCTCTGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCGGGCCTCGAc  >  8:326612/3‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
CCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGAC  >  NC_002947/4960821‑4961105

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