Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,664,594 G→A L9F (CTC→TTC)  gadW ← DNA‑binding transcriptional dual regulator GadW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,664,5940GA95.7% 133.0 / ‑1.5 47L9F (CTC→TTC) gadWDNA‑binding transcriptional dual regulator GadW
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (17/27);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.16e-01

CGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATGCAGCACGAATATATTTTCGCACAGCGTATAGCTTATGTTTATAAAAAAATGGCTGATCT  >  NC_000913/3664352‑3664835
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
cGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCCGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGACTGAACGACGAATGAAGATCACCg                                                                                                                                                                                                                                           >  2:122375/1‑251 (MQ=255)
            tGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATg                                                                                                                                                                                                                                 >  1:27113/1‑249 (MQ=255)
                 ggTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATg                                                                                                                                                                                                                          <  1:46882/251‑1 (MQ=255)
                  gTTAATTAAACAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGa                                                                                                                                                                                                                         <  2:236429/251‑1 (MQ=255)
                          aCAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCt                                                                                                                                                                                                                 <  2:301550/251‑1 (MQ=255)
                              aCGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCACAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTtat                                                                                                                                                                                                              >  1:191426/1‑250 (MQ=255)
                                                      aTTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCACCCTGTCAGCTCGTCGATATCCAGTAGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGACGATCACCGAGTAGACATGAGTACTGATTATCCCTTATATTTCTTACTGCGATTctttcagtt                                                                                                                                                                                     >  2:242071/1‑248 (MQ=255)
                                                      aTTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTatttcaatt                                                                                                                                                                                     <  1:362383/251‑1 (MQ=255)
                                                        tGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTatttcaattt                                                                                                                                                                                    <  1:477653/250‑1 (MQ=255)
                                                             ggTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCAGCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGtat                                                                                                                                                                              <  1:227019/251‑1 (MQ=255)
                                                              gTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATGTCCAGTCGTCGCGTAACCGGTGAACAAAACGCAAAATCGGCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGtata                                                                                                                                                                             <  2:119874/251‑1 (MQ=255)
                                                                         atatTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCATCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTaatgaaatg                                                                                                                                                                  <  2:29299/251‑1 (MQ=255)
                                                                                                       tCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCa                                                                                                                                                                                                                            >  1:270702/1‑163 (MQ=255)
                                                                                                       tCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCa                                                                                                                                                                                                                            <  2:270702/163‑1 (MQ=255)
                                                                                                        cATCGATATCCAGTCGTCGCGTATCCAGTGAACAAAACGAAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCtt                                                                                                                                                 <  1:25836/237‑1 (MQ=255)
                                                                                                        cATCGATATCCAGTCGTCGCGTATCCAGTGAACAAAACGAAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCtt                                                                                                                                                 >  2:25836/1‑237 (MQ=255)
                                                                                                                 ccAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTgg                                                                                                                          <  2:27113/251‑1 (MQ=255)
                                                                                                                       gtcgCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAACTTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTa                                                                                                                                                                                                                      >  2:144802/1‑153 (MQ=255)
                                                                                                                       gtcgCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTa                                                                                                                                                                                                                      <  1:144802/153‑1 (MQ=255)
                                                                                                                         cgcgTATCCGGTGAACAAAACGCAAAAGCCTCTGCCAAATTTTTCGCCAGCAGCAGAATACGCTCGTTATGCAGCGATATTTTTTGCTGGTCATGATAAATATCGAATGAACTACGAATGAAGATCACCGAGCAGACATGATTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTTTACTAATGAAATGATTCCAGACTGTTTCTTAGCGCAAAAACTGccgtatggtggaggacc                                                                                                                   <  2:480281/250‑16 (MQ=255)
                                                                                                                          gcgTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCaa                                                                                                                 <  1:268751/251‑1 (MQ=255)
                                                                                                                              aTCCGGGGCACAGAACGCACAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTAATACTGCGATTATTTCACTTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCCAAAACTGACAGCTTTGGTAGGACCAAAAtt                                                                                                             >  2:93073/1‑251 (MQ=255)
                                                                                                                                       aCAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCcaca                                                                                                     >  1:350045/1‑250 (MQ=255)
                                                                                                                                                 aaaTCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTg                                                                                          >  1:115807/1‑251 (MQ=255)
                                                                                                                                                     cGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTACCCACAAAGTAGTTTTGCAt                                                                                       >  1:88889/1‑250 (MQ=255)
                                                                                                                                                      gTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGccta                                                                                      >  1:426818/1‑247 (MQ=255)
                                                                                                                                                               aTTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAAtt                                                                            <  1:313241/251‑1 (MQ=255)
                                                                                                                                                                 ttttCGCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGGTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAACGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTg                                                                           <  2:110376/250‑1 (MQ=255)
                                                                                                                                                                       ccAGCAGCAGAATAATCTCGTTTTGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAAAGACTAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTAGTTCAATTTCAGTATACTAATGAAATGATGCCATACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGAACAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGAtaataa                                                                    <  2:397839/251‑1 (MQ=255)
                                                                                                                                                                            agcagAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATGTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCTCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTTCATAACAGTTAATTGATAATAAGTaaa                                                               <  2:310338/251‑1 (MQ=255)
                                                                                                                                                                             gcagAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGTGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAAt                                                              <  1:415163/251‑1 (MQ=255)
                                                                                                                                                                                   tACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAACTGATGCCAGACTGTTTCTTAGCGCAAAACCTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTTATAATAAGTAAATgcagc                                                         >  2:56591/1‑250 (MQ=255)
                                                                                                                                                                                    aCTCTCGTTATGCATCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGTACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATGCAGCAc                                                       <  1:372201/251‑1 (MQ=255)
                                                                                                                                                                                       ctcGTTATGCAGCGATATTTGTTGCTGTTCATGATAAATAGCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATgcagcatg                                                      <  1:300011/249‑3 (MQ=255)
                                                                                                                                                                                        tcGTTATGCAGCGATATTTTTTGCTGTTAATGATAAATATCGAATGAACGACGAATGAAGATCACCGTGCAGCCATGAGTCATGATTATACCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTATTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATgcagaaagaat                                                   >  2:302411/1‑244 (MQ=255)
                                                                                                                                                                                           ttATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATGCAACACGAatat                                                 >  2:110584/1‑250 (MQ=255)
                                                                                                                                                                                                          ttttGCTGTTCATGATAAATTTCGAATGTACGACGAATTAAGATCAACGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATGCCGCACGAATATATTTTCGCACAGCGt                                  <  2:341547/250‑1 (MQ=255)
                                                                                                                                                                                                          ttttGCTGTTCATGATAAATATCGAATGAACGTCGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATGCAGCACGAATATATTTTCGCACAGCGt                                  <  2:115807/250‑1 (MQ=255)
                                                                                                                                                                                                              gCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCaaaa                                                                                                                                        <  2:470798/144‑1 (MQ=255)
                                                                                                                                                                                                              gCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCaaaa                                                                                                                                        >  1:470798/1‑144 (MQ=255)
                                                                                                                                                                                                               cTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTa                                                                                                                                                >  1:364764/1‑135 (MQ=255)
                                                                                                                                                                                                               cTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTa                                                                                                                                                <  2:364764/135‑1 (MQ=255)
                                                                                                                                                                                                                   tCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATGCAGCACGAATATATTTTCGCACAGCGTATAGCTTAt                         <  1:98302/250‑1 (MQ=255)
                                                                                                                                                                                                                          aaaaatCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATGCAGCACGAATATATTTTCGCACAGCGTATAGCTTATGTTtata                  >  1:415018/5‑250 (MQ=255)
                                                                                                                                                                                                                                         gacgaAGGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAAGTTATCCACAAAGTAGTTTTGCATAACAGTTAAGTGATAATAAGTAAATGCAGCACGAATTTATTTTCGCACAGCGTATAGCTAATGTTTATAAAAAAATGGCTGATCt  <  2:495950/251‑1 (MQ=255)
                                                                                                                                                                                                                                          acgcaTGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAACCTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATGCAGCACGAATATATTTTCGCACAGCGTATAGCTTAt                         >  2:168757/5‑227 (MQ=255)
                                                                                                                                                                                                                                          acgaATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATGCAGCACGAATATATTTTCGCACAGCGTATAGCTTAt                         <  1:168757/227‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
CGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATGCAGCACGAATATATTTTCGCACAGCGTATAGCTTATGTTTATAAAAAAATGGCTGATCT  >  NC_000913/3664352‑3664835

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