Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,945,329 T→C 5.3% intergenic (+644/+154) 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,3290TC 5.3% 272.8 / 2.0 89intergenic (+644/+154)PP_5673/flhBD‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base T (38/45);  new base C (3/2);  total (41/48)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAA  >  NC_002947/4945186‑4945465
                                                                                                                                               |                                                                                                                                        
gctgTGAATCGGGGTGTTCTCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCt                                                                                                                                    <  2:2066716/150‑1 (MQ=255)
 ctgTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTc                                                                                                                                   <  2:1713293/150‑1 (MQ=255)
   caataTCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTTCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGa                                                                                                                                 >  2:1126514/5‑150 (MQ=255)
     gAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGata                                                                                                                               <  1:1002587/150‑1 (MQ=255)
       aTCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATAtt                                                                                                                             >  2:1376089/1‑150 (MQ=255)
       aTCGGGGTGTTCGCCGTGAGGTTTTTGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATAtt                                                                                                                             <  1:788917/150‑1 (MQ=255)
        tCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTTCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTg                                                                                                                            <  1:1126514/150‑1 (MQ=255)
        tCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTg                                                                                                                            <  2:800714/150‑1 (MQ=255)
        tCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTg                                                                                                                            <  2:278860/150‑1 (MQ=255)
         cGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTgc                                                                                                                           <  2:194727/150‑1 (MQ=255)
           gggTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTgtgc                                                                                                                         <  1:2329282/150‑1 (MQ=255)
           gggTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTgcgc                                                                                                                         <  2:2171709/150‑1 (MQ=255)
           gggTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTgcgc                                                                                                                         <  1:825995/150‑1 (MQ=255)
            ggTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGcc                                                                                                                        >  1:1536727/1‑150 (MQ=255)
               gtTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTa                                                                                                                     <  1:1376089/150‑1 (MQ=255)
                ttCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTacac                                                                                                                                      <  2:1474104/132‑1 (MQ=255)
                ttCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTacac                                                                                                                                      >  1:1474104/1‑132 (MQ=255)
                     cGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTacaaaca                                                                                                               >  1:557023/1‑150 (MQ=255)
                      gTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTacaaacaa                                                                                                              <  1:2225798/150‑1 (MQ=255)
                         aGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGc                                                                                                           <  1:2370346/150‑1 (MQ=255)
                           gTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCgg                                                                                                         <  2:2566123/150‑1 (MQ=255)
                           gTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCgg                                                                                                         <  2:557023/150‑1 (MQ=255)
                            tttttGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTc                                                                                                                                   <  1:2102231/123‑1 (MQ=255)
                            tttttGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTc                                                                                                                                   >  2:2102231/1‑123 (MQ=255)
                                  tCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGat                                                                                                                                <  1:2469653/120‑1 (MQ=255)
                                  tCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGat                                                                                                                                >  2:2469653/1‑120 (MQ=255)
                                  tCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTcgcgcg                                                                                                  >  1:1378942/1‑150 (MQ=255)
                                         tgAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGcc                                                                                                                        <  1:2714391/121‑1 (MQ=255)
                                         tgAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGcc                                                                                                                        >  2:2714391/1‑121 (MQ=255)
                                         tgAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCggg                                                                                           >  2:284200/1‑150 (MQ=255)
                                         tgAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCggg                                                                                           <  2:1378942/150‑1 (MQ=255)
                                              tCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACag                                                                                      <  1:284200/150‑1 (MQ=255)
                                               cGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACagg                                                                                     <  2:487577/150‑1 (MQ=255)
                                                 aGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACaggag                                                                                   <  1:74954/150‑1 (MQ=255)
                                                  gcgcCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACaggagg                                                                                  >  2:2524721/1‑150 (MQ=255)
                                                        cccgcGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTc                                                                                                                                   >  1:1161190/1‑95 (MQ=255)
                                                        cccgcGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTc                                                                                                                                   <  2:1161190/95‑1 (MQ=255)
                                                         ccgcGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGc                                                                           <  1:2524721/150‑1 (MQ=255)
                                                           gcgcgGCGCATCGCGAGCAAGCTAGCTCCTACATTTGTTTCGGGTCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGccc                                                                         <  2:466720/150‑1 (MQ=255)
                                                              gggcgCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGGACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGaa                                                                      <  1:2628819/149‑1 (MQ=255)
                                                                gcgcATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGc                                                                           <  2:515969/143‑1 (MQ=255)
                                                                gcgcATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGc                                                                           >  1:515969/1‑143 (MQ=255)
                                                                gcgcATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAAc                                                                    <  1:2473768/150‑1 (MQ=255)
                                                                   cATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACCCCTCGATATTGCGCGGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACaaa                                                                 <  2:2847783/150‑1 (MQ=255)
                                                                    aTCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCTTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAAt                                                                <  1:1250937/150‑1 (MQ=255)
                                                                    aTCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAAt                                                                <  2:2776504/150‑1 (MQ=255)
                                                                    aTCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAAt                                                                <  2:1536727/150‑1 (MQ=255)
                                                                       gcgAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTa                                                                                                                                         >  2:1923968/1‑74 (MQ=255)
                                                                       gcgAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTa                                                                                                                                         <  1:1923968/74‑1 (MQ=255)
                                                                        cgAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAg                                                            <  2:394609/150‑1 (MQ=255)
                                                                         gAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAgg                                                           <  1:432969/150‑1 (MQ=255)
                                                                                     tcCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGcc                                                                          >  1:846475/1‑123 (MQ=255)
                                                                                     tcCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGcc                                                                          <  2:846475/123‑1 (MQ=255)
                                                                                       cTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGaaa                                                                     >  2:2009705/1‑126 (MQ=255)
                                                                                       cTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGaaa                                                                     <  1:2009705/126‑1 (MQ=255)
                                                                                         aCATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAg                                                            <  1:2248568/133‑1 (MQ=255)
                                                                                         aCATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAg                                                            >  2:2248568/1‑133 (MQ=255)
                                                                                         aCATTTCTTTCGGGCCAGTTATTCCTGCGGGATTTACTCGCATCCGCCTTCGCGCACACCTCGCTATAGCGCCGTACAAACAAGGCGGTCGCGCGCGCCCGGCACAGGAGGGACTCGCCCGAAACAAATGTAGGAGCGAGCTCGCCcgcc                                           >  2:469491/1‑149 (MQ=21)
                                                                                           aTTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGAt                                         >  2:1139736/1‑150 (MQ=255)
                                                                                            tttGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATg                                        >  1:802570/1‑150 (MQ=255)
                                                                                                tttCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACAGCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATCcgcc                                    >  1:1122053/1‑150 (MQ=255)
                                                                                                  tCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTg                                                                             <  2:1654991/107‑1 (MQ=255)
                                                                                                  tCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTg                                                                             >  1:1654991/1‑107 (MQ=255)
                                                                                                    gggCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCcgcgc                                >  1:306090/1‑150 (MQ=255)
                                                                                                     ggCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCcgcgcg                               >  1:196516/1‑150 (MQ=255)
                                                                                                      gCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGa                                                                                >  2:895470/1‑100 (MQ=2)
                                                                                                      gCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGa                                                                                <  1:895470/100‑1 (MQ=2)
                                                                                                          attATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACaaac                                                                >  2:665247/2‑110 (MQ=2)
                                                                                                          attATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACaaac                                                                <  1:665247/110‑2 (MQ=2)
                                                                                                           ttATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGggcggc                         >  2:1512412/1‑150 (MQ=255)
                                                                                                           ttATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGggcggc                         >  2:1073317/1‑150 (MQ=255)
                                                                                                             aTTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGgcgc                       >  1:1305869/1‑150 (MQ=255)
                                                                                                                 cTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGCGAGCGCTTGCTCGCTATGCGCCGCCCGCGCGgcgcatcg                   >  2:28183/1‑146 (MQ=37)
                                                                                                                  tGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGAt                  >  2:119955/1‑150 (MQ=255)
                                                                                                                  tGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCCATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCCCGCGCGgcgcatcgc                  <  1:28183/150‑6 (MQ=37)
                                                                                                                    cGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGTGCGGCGCTCGAtct                >  2:2935405/1‑150 (MQ=255)
                                                                                                                           tACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGc         >  1:2668093/1‑150 (MQ=255)
                                                                                                                              gcgcgACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAAc      >  1:136873/1‑150 (MQ=255)
                                                                                                                              gcgcgACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAAc      >  2:785484/1‑150 (MQ=255)
                                                                                                                                gcgACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTcaca            <  2:471596/142‑1 (MQ=255)
                                                                                                                                gcgACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTcaca            >  1:471596/1‑142 (MQ=255)
                                                                                                                                  gACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCaa  >  1:2406837/1‑150 (MQ=255)
                                                                                                                                  gACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCaa  >  2:2673658/1‑150 (MQ=255)
                                                                                                                                       ccTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGcgcc                                    >  1:2758483/1‑111 (MQ=255)
                                                                                                                                       ccTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGcgcc                                    <  2:2758483/111‑1 (MQ=255)
                                                                                                                                        cTTCGCGTACACCTCGATATTGCTCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTa                                                             <  1:144281/85‑1 (MQ=255)
                                                                                                                                        cTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTa                                                             >  2:144281/1‑85 (MQ=255)
                                                                                                                                          tCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAgg          <  1:2130389/134‑1 (MQ=255)
                                                                                                                                          tCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAgg          >  2:2130389/1‑134 (MQ=255)
                                                                                                                                               |                                                                                                                                        
GCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAA  >  NC_002947/4945186‑4945465

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