Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,473,680 C→T G125E (GGG→GAG)  bdcA ← c‑di‑GMP‑binding biofilm dispersal mediator protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,473,6800CT98.6% 232.9 / ‑4.2 69G125E (GGG→GAG) bdcAc‑di‑GMP‑binding biofilm dispersal mediator protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (29/39);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

TGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCA  >  NC_000913/4473546‑4473815
                                                                                                                                      |                                                                                                                                       
tGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCt                                                                                                                                    <  1:396289/140‑1 (MQ=255)
 gATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCtt                                                                                                                                   >  2:352814/1‑140 (MQ=255)
 gATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCtt                                                                                                                                   >  2:456958/1‑140 (MQ=255)
 gATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCtt                                                                                                                                   >  2:547847/1‑140 (MQ=255)
  aTCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTc                                                                                                                                  >  2:141827/1‑140 (MQ=255)
  aTCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTc                                                                                                                                  >  1:304527/1‑140 (MQ=255)
             tCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGcc                                                                                                                       >  1:697051/1‑140 (MQ=255)
             tCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGcc                                                                                                                       >  2:57328/1‑140 (MQ=255)
               cagaaTCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCgg                                                                                                                     <  1:271147/137‑1 (MQ=255)
               cAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCATCAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCgg                                                                                                                     <  1:806527/140‑1 (MQ=255)
               cAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCgg                                                                                                                     >  1:380839/1‑140 (MQ=255)
               cAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCgg                                                                                                                     <  2:630103/140‑1 (MQ=255)
               cAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCgg                                                                                                                     <  2:399680/140‑1 (MQ=255)
                 aaaTCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGggc                                                                                                                   <  2:443976/140‑1 (MQ=255)
                 aaaTCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGggc                                                                                                                   <  1:136319/140‑1 (MQ=255)
                 aaaTCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGAAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGggc                                                                                                                   <  2:305214/140‑1 (MQ=255)
                    tCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGggcggc                                                                                                                <  1:141827/140‑1 (MQ=255)
                       cGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTc                                                                                                             >  2:833510/1‑140 (MQ=255)
                       cGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTc                                                                                                             >  1:811024/1‑140 (MQ=255)
                            cAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACag                                                                                                        <  1:263943/140‑1 (MQ=255)
                            cAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACag                                                                                                        <  1:546531/140‑1 (MQ=255)
                             aGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACaga                                                                                                       >  2:674847/1‑140 (MQ=255)
                              gCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCTGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACagag                                                                                                      <  2:761019/140‑1 (MQ=255)
                              gCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACagag                                                                                                      <  2:811024/140‑1 (MQ=255)
                              gCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACagag                                                                                                      <  2:380839/140‑1 (MQ=255)
                              gCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACagag                                                                                                      <  2:102393/140‑1 (MQ=255)
                              gCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACagag                                                                                                      <  1:784332/140‑1 (MQ=255)
                              gCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACagag                                                                                                      <  1:167990/140‑1 (MQ=255)
                                ccGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACTGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGc                                                                                                    <  1:453689/140‑1 (MQ=255)
                                ccGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGc                                                                                                    <  2:687324/140‑1 (MQ=255)
                                      gcCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATa                                                                                              <  1:883020/140‑1 (MQ=255)
                                      gcCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATa                                                                                              <  2:894644/140‑1 (MQ=255)
                                      gcCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATa                                                                                              <  2:697051/140‑1 (MQ=255)
                                         aTGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAgg                                                                                           <  1:813335/140‑1 (MQ=255)
                                           gCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAg                                                                                         <  1:352814/140‑1 (MQ=255)
                                             cTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCa                                                                                       <  1:890155/140‑1 (MQ=255)
                                                  aGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAat                                                                                  >  2:418022/1‑140 (MQ=255)
                                                  aGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAat                                                                                  >  1:885282/1‑140 (MQ=255)
                                                      cAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTg                                                                              >  1:772348/1‑140 (MQ=255)
                                                      cAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTg                                                                              >  2:743595/1‑140 (MQ=255)
                                                            tGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTg                                                                        <  1:141357/140‑1 (MQ=255)
                                                                 gcggcATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAg                                                                   <  1:491185/140‑1 (MQ=255)
                                                                 gcggcATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAg                                                                   >  1:161852/1‑140 (MQ=255)
                                                                 gcggcATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAg                                                                   <  1:57328/140‑1 (MQ=255)
                                                                   ggcATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGc                                                                 <  2:706696/140‑1 (MQ=255)
                                                                    gcATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCg                                                                <  1:418022/140‑1 (MQ=255)
                                                                             gcCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATAtc                                                       >  2:232134/1‑140 (MQ=255)
                                                                               cacgCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATAtcgt                                                     <  2:52725/137‑1 (MQ=255)
                                                                               cATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATAtcgt                                                     <  2:449164/140‑1 (MQ=255)
                                                                                     tGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCAt                                               >  2:103809/1‑140 (MQ=255)
                                                                                       cAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCAttt                                             <  1:547847/140‑1 (MQ=255)
                                                                                       cAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCAttt                                             <  2:304527/140‑1 (MQ=255)
                                                                                       cAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCAttt                                             >  1:74371/1‑140 (MQ=255)
                                                                                                tACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAgg                                    <  2:94861/140‑1 (MQ=255)
                                                                                                 aCGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAggg                                   >  1:772555/1‑140 (MQ=255)
                                                                                                      cGCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCg                              >  2:552431/1‑140 (MQ=255)
                                                                                                       gCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGc                             >  1:380986/1‑140 (MQ=255)
                                                                                                       gCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGc                             >  1:618514/1‑140 (MQ=255)
                                                                                                       gCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGc                             >  1:885229/1‑140 (MQ=255)
                                                                                                       gCCATTCACGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGc                             >  1:835038/1‑140 (MQ=255)
                                                                                                               cGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACg                     >  1:514869/1‑140 (MQ=255)
                                                                                                               cGGAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACg                     >  2:29690/1‑140 (MQ=255)
                                                                                                                 gAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGcc                   >  2:672296/1‑140 (MQ=255)
                                                                                                                 gAGCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGcc                   <  1:103809/140‑1 (MQ=255)
                                                                                                                   gCCGATGATTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCaa                 >  2:571146/1‑140 (MQ=255)
                                                                                                                       atgatTAAGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATAcc             >  2:332090/1‑140 (MQ=255)
                                                                                                                             aaGATGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCAtt       <  2:813293/140‑1 (MQ=255)
                                                                                                                                aTGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATTCCTGCAtttgc    <  2:772348/140‑3 (MQ=255)
                                                                                                                                 tGCGCTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTacc   <  2:393778/140‑1 (MQ=255)
                                                                                                                                  gcgcTCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTacca  <  1:674847/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
TGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCA  >  NC_000913/4473546‑4473815

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