Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 652,767 G→C 6.8% intergenic (‑104/+14) dsbD‑I ← / ← ctpL thiol/disulfide interchange protein 1/methyl‑accepting chemotaxis protein CtpL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947652,7670GC 6.8% 224.4 / 5.8 74intergenic (‑104/+14)dsbD‑I/ctpLthiol/disulfide interchange protein 1/methyl‑accepting chemotaxis protein CtpL
Reads supporting (aligned to +/- strand):  ref base G (33/36);  new base C (2/3);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

GGCTGGCCAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCTGTGGT  >  NC_002947/652621‑652911
                                                                                                                                                  |                                                                                                                                                
ggCTGGCCAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGC                                                                                                                                                <  8:246633‑M1/149‑1 (MQ=255)
 gCTGGCCAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCG                                                                                                                                               >  8:182303‑M1/1‑149 (MQ=255)
  cTGGCCAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGG                                                                                                                                              <  4:130052‑M1/149‑1 (MQ=255)
      ccAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCT                                                                                                                                          <  4:40537‑M1/149‑1 (MQ=255)
       cagcagTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGCAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTC                                                                                                                                         >  8:81898‑M1/1‑149 (MQ=255)
         gcagTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGGAGCGCCGCCCGCGCGGCGCTtcg                                                                                                                                       <  1:60685‑M1/149‑4 (MQ=255)
          cagTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGAT                                                                                                                                      <  5:126456‑M1/149‑1 (MQ=255)
          cagTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGAT                                                                                                                                      <  7:182303‑M1/149‑1 (MQ=255)
          cagTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCCCGCGCGGCGCTtcgc                                                                                                                                      <  5:303277‑M1/149‑5 (MQ=255)
             tAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCT                                                                                                                                    >  2:216782‑M1/1‑148 (MQ=255)
               ggAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCa                                                                                                                                  <  7:64559/148‑1 (MQ=255)
                    cgcttgtaaaaCAGGAGGCGAATGGCGGACTGTCGTGGGGGAACGTTGCCGGCTTTTTTAGGGATGCGGAGGAATAATTGGCCCTAAACATACGTAGGAGCGGTTTTACGCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCt                                                                                                                            <  6:234585/142‑1 (MQ=16)
                      cATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTg                                                                                                                           >  3:148079/1‑148 (MQ=255)
                      cATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGa                                                                                                                          <  2:284160/149‑1 (MQ=255)
                      cATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGAT                                                                                                                                      <  1:334521‑M1/137‑1 (MQ=255)
                      cATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGAT                                                                                                                                      >  2:334521‑M1/1‑137 (MQ=255)
                       aTGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGaa                                                                                                                         >  5:322258/1‑149 (MQ=255)
                          aaaaaCAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGc                                                                                                                      <  5:324178/149‑1 (MQ=255)
                           aaaaCAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCt                                                                                                                     <  6:223832/149‑1 (MQ=255)
                            aaaCAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTc                                                                                                                    <  6:322258/149‑1 (MQ=255)
                                                  tCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCAcc                                                                                              >  3:190887/1‑149 (MQ=255)
                                                  tCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCAcc                                                                                              >  2:324923/1‑149 (MQ=255)
                                                       aGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCa                                                                                         <  1:216782/149‑1 (MQ=255)
                                                            gAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGtt                                                                                                       >  2:213421/1‑130 (MQ=255)
                                                            gAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGtt                                                                                                       <  1:213421/130‑1 (MQ=255)
                                                            gAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTc                                                                                    <  2:217242/149‑1 (MQ=255)
                                                             aaCGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCAGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCaagatcg                                                                                   >  5:33347/1‑143 (MQ=255)
                                                             aaCGTGCCCGGCATGATAAAGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGccc                                                                                          <  6:33347/142‑1 (MQ=255)
                                                              aCGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTcag                                                                                  >  6:333976/1‑149 (MQ=255)
                                                                               aCGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCTCGCGGGCGGCGCTCTATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTATCCGCCCAGGCTCAGCAGCTGTTCACCCTGtt                                                                 <  5:134969/149‑1 (MQ=255)
                                                                                     ccgcaGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaaccccacccaggcgctgaaaaaaccccggcgaacccctcccacgcg                                                            >  1:267197‑M1/2‑77 (MQ=255)
                                                                                       gcAGGAAAAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGcccc                                                         >  5:112852/1‑149 (MQ=255)
                                                                                         aGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaataccccggagaacaccacccaggcgctga                                                        >  8:35772‑M1/1‑73 (MQ=255)
                                                                                           gAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaa                                                     >  4:164402‑M1/1‑71 (MQ=255)
                                                                                             ataatTGGCCCGAAACAAACGTAGGAGCTGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaa                                                     <  7:35772‑M1/147‑79 (MQ=255)
                                                                                             ataatTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGcc                                                    >  5:76962/1‑148 (MQ=255)
                                                                                                atTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATg                                                <  1:324923/149‑1 (MQ=255)
                                                                                                atTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATg                                                <  6:76962/149‑1 (MQ=255)
                                                                                                 ttGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATgc                                               >  2:248282/1‑149 (MQ=255)
                                                                                                 ttGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATgc                                               >  6:176987/1‑149 (MQ=255)
                                                                                                      ccGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTg                                           >  7:141163/1‑148 (MQ=255)
                                                                                                      ccGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTg                                           >  2:331248/1‑148 (MQ=255)
                                                                                                      ccGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTg                                           >  2:336409/1‑148 (MQ=255)
                                                                                                            caaccGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGAGCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTg                                    <  7:81898/149‑1 (MQ=255)
                                                                                                            caaacGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACg                                                             <  1:326837/124‑1 (MQ=255)
                                                                                                            caaacGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACg                                                             >  2:326837/1‑124 (MQ=255)
                                                                                                             aaacGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAg                                                                                                    <  4:289040/84‑1 (MQ=255)
                                                                                                             aaacGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAg                                                                                                    >  3:289040/1‑84 (MQ=255)
                                                                                                             aaacGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGt                                   >  1:135797/1‑149 (MQ=255)
                                                                                                              aacGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTc                                  >  5:75729/1‑149 (MQ=255)
                                                                                                              aacGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTc                                  <  4:148079/149‑1 (MQ=255)
                                                                                                               acGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCt                                 <  5:176987/149‑1 (MQ=255)
                                                                                                                cGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCtt                                <  8:141163/149‑1 (MQ=255)
                                                                                                                 gTAGGAGCGGGTTTACCCGCGAAGCGCCGCTCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCtt                                >  8:239248/1‑148 (MQ=255)
                                                                                                                  tAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCa                              <  6:112852/149‑1 (MQ=255)
                                                                                                                     gAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcacccgcccaacacctcccagcccctccacgcctaccgccaaactcctaacagcaa                            >  3:182499‑M1/1‑45 (MQ=255)
                                                                                                                         gggTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGtt                       >  7:83079/1‑149 (MQ=255)
                                                                                                                            tttACCCGCGAAGCGCCGCTCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGt                    <  7:239248/149‑1 (MQ=255)
                                                                                                                            tttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGt                    <  2:70747/149‑1 (MQ=255)
                                                                                                                            tttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGt                    >  1:70747/1‑149 (MQ=255)
                                                                                                                                 ccGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTc                      <  1:175314/142‑1 (MQ=255)
                                                                                                                                 ccGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTAACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTc                      >  2:175314/1‑142 (MQ=255)
                                                                                                                                   gcgAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCcg             >  6:76874/1‑149 (MQ=255)
                                                                                                                                   gcgAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCcg             >  6:81836/1‑149 (MQ=255)
                                                                                                                                    tcgaGCGCCGCCCGCGCGGCGCTTGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCcgct           >  8:127741/4‑149 (MQ=255)
                                                                                                                                     gAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCcgct           >  6:174651/1‑149 (MQ=255)
                                                                                                                                       aGCGCCGCCCGCGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAg                                        >  1:93306/1‑118 (MQ=255)
                                                                                                                                       aGCGCCGCCCGCGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAg                                        <  2:93306/118‑1 (MQ=255)
                                                                                                                                            cgcgcgGGGGGCGCTCGATGTCTGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCGCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCGTCGCGAATGGGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCtgtg    <  7:114064/149‑1 (MQ=255)
                                                                                                                                            cgcgcgGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCtgtg    <  1:248282/149‑1 (MQ=255)
                                                                                                                                            cgcgcgGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCtgtg    <  4:190887/149‑1 (MQ=255)
                                                                                                                                             gcgcgGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgccaaacacctacccgccaaagcgcgaattccgcccattgccct   <  3:164402‑M1/149‑129 (MQ=255)
                                                                                                                                              cgcgGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCTGTGGt  <  1:336409/149‑1 (MQ=255)
                                                                                                                                              cgcgGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCTGTGGt  <  1:331248/149‑1 (MQ=255)
                                                                                                                                               gcgggcggCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCTGTGGt  >  7:96149/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GGCTGGCCAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCTGTGGT  >  NC_002947/652621‑652911

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