Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,945,332 A→G 7.3% intergenic (+647/+151) PP_5673 → / ← flhB D‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,945,3320AG 7.3% 188.2 / 2.3 66intergenic (+647/+151)PP_5673/flhBD‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base A (32/28);  new base G (2/3);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACC  >  NC_002947/4945194‑4945468
                                                                                                                                          |                                                                                                                                        
tCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGCTCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTg                                                                                                                               <  1:498575/150‑1 (MQ=255)
      tgtTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGt                                                                                                                         <  1:241841/150‑1 (MQ=255)
         tCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTaca                                                                                                                      <  1:177452/150‑1 (MQ=255)
            ctgTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTTTTACATTTGTTTCGTGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTacaaac                                                                                                                   <  2:1758383/148‑1 (MQ=255)
              gTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTacaaacaa                                                                                                                 <  2:1387714/150‑1 (MQ=255)
               tGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAg                                                                                                                <  2:736915/150‑1 (MQ=255)
                gAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGa                                                                                                                                    >  2:1729168/1‑129 (MQ=255)
                gAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGa                                                                                                                                    <  1:1729168/129‑1 (MQ=255)
                gAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAgg                                                                                                               <  1:458286/150‑1 (MQ=255)
                 aGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGc                                                                                                              <  2:1939068/150‑1 (MQ=255)
                  ggTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGCCCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCg                                                                                                             >  2:356559/1‑150 (MQ=255)
                  ggTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCg                                                                                                             <  1:2008859/150‑1 (MQ=255)
                                gtgAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCgg                                                                                               <  1:356559/150‑1 (MQ=255)
                                 tgAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCggg                                                                                              >  1:1150330/1‑150 (MQ=255)
                                 tgAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCggg                                                                                              >  1:116752/1‑150 (MQ=255)
                                     aTCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGc                                                                                                              <  2:2027717/130‑1 (MQ=255)
                                     aTCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGcaca                                                                                          <  2:2062347/150‑1 (MQ=255)
                                     aTCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACAGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGc                                                                                                              >  1:2027717/1‑130 (MQ=255)
                                           cgccgcCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAggg                                                                                    <  1:1197860/150‑1 (MQ=255)
                                               gcccgcGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTg                                                                                <  2:116752/150‑1 (MQ=255)
                                               gcccgcGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTg                                                                                <  2:1150330/150‑1 (MQ=255)
                                                 ccgcGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGc                                                                              <  2:272852/150‑1 (MQ=255)
                                                 ccgcGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGc                                                                              <  1:2254933/150‑1 (MQ=255)
                                                  cgcgcgGCAGATCGTGAGCATGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGAGACCGCCTTCGCGTACCCCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCCGGCACAGGAGGGACTGGct                                                                             <  2:598906/150‑2 (MQ=25)
                                                    cgcgGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCg                                                                           <  1:589101/150‑1 (MQ=255)
                                                      cggcgCATCGCGAGCAATCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGTTATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGaa                                                                         <  2:2183406/150‑1 (MQ=255)
                                                             tCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGata                                                                                                                                  <  2:1186859/86‑1 (MQ=255)
                                                             tCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGata                                                                                                                                  >  1:1186859/1‑86 (MQ=255)
                                                             tCGCAAGCAAGCTCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGGCCGCCTTCGCGTACGCCTCGATATTGCGACGTACAAACCAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATg                                                                  <  2:1332044/150‑1 (MQ=21)
                                                                         tcgctcCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCt                                                      >  1:36364/1‑150 (MQ=255)
                                                                          cgctcCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACGGGCCCGAAACAAATGTAGGAGCGAGCtt                                                     <  1:437875/150‑1 (MQ=255)
                                                                           gctcCTACATTTGTTTCGGGCCAGTTATCCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGACGTACAAACAAGGCGGCCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTg                                                    <  1:472518/150‑1 (MQ=255)
                                                                                    tttGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTAGCGCGCGCCGGGCACAGGAGGGACTGGCCCCAAACAAATGTAGGAGCGGGCTTGCTCGCGATg                                           >  1:1568693/1‑150 (MQ=255)
                                                                                       gTTTCGGGCCAGTTATTCCTGCGGGGTGTGCTCGCGACTGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATgcgc                                        >  1:643469/1‑150 (MQ=25)
                                                                                        tttCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGcgcc                                       >  2:2206822/1‑150 (MQ=255)
                                                                                        tttCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGACCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGcgcc                                       >  1:2239734/1‑150 (MQ=255)
                                                                                        tttCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGG‑CCGAAACAAACGTAGgagcggg                                                        >  1:388577/1‑130 (MQ=2)
                                                                                        tttCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGG‑CCGAAACAAACGTAGgagcggg                                                        <  2:388577/132‑3 (MQ=2)
                                                                                              gCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCgcgg                                 >  1:773956/1‑150 (MQ=255)
                                                                                              gCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGccgaaa                                                                         >  2:2106965/1‑106 (MQ=2)
                                                                                              gCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGccgaaa                                                                         <  1:2106965/110‑5 (MQ=2)
                                                                                               ccAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCgcggg                                >  2:1383224/1‑150 (MQ=255)
                                                                                               ccAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCgcggg                                >  2:2056370/1‑150 (MQ=255)
                                                                                               ccAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGT‑GGAGCGAGCTTGCTCGCGATGCGCCGCgcgggc                               >  1:1178852/1‑150 (MQ=255)
                                                                                                cAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGTAGCGAGCTTGCTAGCGATGCGCCGCgcgggc                               >  1:693975/1‑150 (MQ=255)
                                                                                                  gTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGc                                                                              >  2:2263949/1‑101 (MQ=255)
                                                                                                  gTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGc                                                                              <  1:2263949/101‑1 (MQ=255)
                                                                                                      ttCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCt                         >  1:1576657/1‑150 (MQ=255)
                                                                                                        ccTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCg                       >  2:1468265/1‑150 (MQ=255)
                                                                                                         cTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGa                      >  1:1167462/1‑150 (MQ=255)
                                                                                                         cTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGa                      >  2:58487/1‑150 (MQ=255)
                                                                                                          tGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGAt                     >  1:321371/1‑150 (MQ=255)
                                                                                                            cGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGAtct                   >  2:655224/1‑150 (MQ=255)
                                                                                                                 tttACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGT‑GGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAgg             <  2:1178852/150‑1 (MQ=255)
                                                                                                                 tttACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGAGCGCCGGGCACAGGAGGGACTGGCCCGaa                                                                         >  1:68841/1‑91 (MQ=255)
                                                                                                                 tttACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGAGCGCCGGGCACAGGAGGGACTGGCCCGaa                                                                         <  2:68841/91‑1 (MQ=255)
                                                                                                                   tACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGc            >  1:1352349/1‑150 (MQ=255)
                                                                                                                      gcgcgACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGga                                                             >  2:225141/1‑98 (MQ=255)
                                                                                                                      gcgcgACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGga                                                             <  1:225141/98‑1 (MQ=255)
                                                                                                                      gcgcgACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAAc         >  1:581699/1‑150 (MQ=255)
                                                                                                                        gcgACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGAATGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGc       >  2:113103/1‑150 (MQ=255)
                                                                                                                         cgACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGCGCTTCCTCCCGATGCGCCGCGCGGGCGGCGCGCTAGCTCACAGGCAACGCa      >  2:1423556/1‑150 (MQ=25)
                                                                                                                             cGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAAcc  >  2:431940/1‑150 (MQ=255)
                                                                                                                             cGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAAcc  >  2:73923/1‑150 (MQ=255)
                                                                                                                               ccTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACt                                                                                 <  2:2155754/69‑1 (MQ=255)
                                                                                                                               ccTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACt                                                                                 >  1:2155754/1‑69 (MQ=255)
                                                                                                                                          |                                                                                                                                        
TCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACC  >  NC_002947/4945194‑4945468

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