Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,678,713 C→T 5.9% D159N (GAT→AAT)  yhjH ← cyclic‑di‑GMP phosphodiesterase, FlhDC‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,678,7130CT 5.9% 216.5 / 2.4 69D159N (GAT→AAT) yhjHcyclic‑di‑GMP phosphodiesterase, FlhDC‑regulated
Reads supporting (aligned to +/- strand):  ref base C (31/33);  new base T (2/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

AGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGACGCAGGGCGAT  >  NC_000913/3678566‑3678854
                                                                                                                                                   |                                                                                                                                             
aGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAatca                                                                                                                                              >  7:138054/1‑149 (MQ=255)
  cGATTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAatcatc                                                                                                                                            <  1:125645/149‑1 (MQ=255)
    aTTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCa                                                                                                                                          <  4:203331/149‑1 (MQ=255)
    aTTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCa                                                                                                                                          <  1:212887/149‑1 (MQ=255)
        aTTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCca                                                                                                                                      >  3:92912/1‑149 (MQ=255)
        aTTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCca                                                                                                                                      >  5:196173/1‑149 (MQ=255)
              tGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCgg                                                                                                                                <  7:186043/149‑1 (MQ=255)
                tAAAAGCTGATAAAAGCGTGTGCGGCCTTCCGGCGACTGACGCAGCATCACAAACCGTTCTCGCGCGATTTTGATGTAGTCATACCGCACTTCACTTACCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGc                                                                                                                              <  2:116088/149‑1 (MQ=255)
                tAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGc                                                                                                                              >  5:382502/1‑149 (MQ=255)
                    atcgtagaAAAGCGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACATTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCaaa                                                                                                                          <  8:275396/144‑1 (MQ=255)
                    aGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCaaa                                                                                                                          <  2:97540/149‑1 (MQ=255)
                      cTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTTTATCAGCATCACAAACAGTTCTAGCGCGATTTTGATGTAGTCATATCTCTCTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTAACAAAATGATCCAGCCACAGCGGGCCAAAtt                                                                                                                        <  6:271961/149‑1 (MQ=255)
                        gagaAAAGAGTGTGCGTCCTTCCGGCGACTGCCGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTca                                                                                                                      >  6:82546/1‑149 (MQ=255)
                         agaAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGTTGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTcac                                                                                                                     <  4:317716/149‑1 (MQ=255)
                         agaAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTca                                                                                                                      <  6:222902/148‑1 (MQ=255)
                          gaAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTcaca                                                                                                                    >  2:351348/1‑149 (MQ=255)
                           aaaaGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTcaca                                                                                                                    >  1:167193/1‑148 (MQ=255)
                                agTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGa                                                                                                              <  6:10132/149‑1 (MQ=255)
                                    tgCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCg                                                                                                               >  6:129173/1‑144 (MQ=255)
                                    tgCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCg                                                                                                               <  5:129173/144‑1 (MQ=255)
                                      cGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCaaa                                                                                                        <  7:173000/149‑1 (MQ=255)
                                      cGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCa                                                                                                          >  8:9370/1‑147 (MQ=255)
                                          cTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGtt                                                                                                    <  7:112305/149‑1 (MQ=255)
                                          cTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGtt                                                                                                    <  2:167193/149‑1 (MQ=255)
                                                gCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCg                                                                                                               <  8:89067/132‑1 (MQ=255)
                                                gCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCg                                                                                                               >  7:89067/1‑132 (MQ=255)
                                                  gACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAatca                                                                                                                                              <  6:334019/99‑1 (MQ=255)
                                                  gACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAatca                                                                                                                                              >  5:334019/1‑99 (MQ=255)
                                                   aCTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTc                                                                                           <  7:294848/149‑1 (MQ=255)
                                                    cTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCg                                                                                          <  3:374756/149‑1 (MQ=255)
                                                    cTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCg                                                                                          <  6:307770/149‑1 (MQ=255)
                                                    cTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCg                                                                                          <  5:315657/149‑1 (MQ=255)
                                                     tGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTca                                                                                                                      <  7:216703/120‑1 (MQ=255)
                                                     tGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTca                                                                                                                      >  8:216703/1‑120 (MQ=255)
                                                     tGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGATTCTTTCgg                                                                                         >  2:84586/1‑149 (MQ=255)
                                                       aCGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGc                                                                                        <  5:38449/148‑1 (MQ=255)
                                                          cagcaTCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGAc                                                                                    <  6:382502/149‑1 (MQ=255)
                                                           ggcaTCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTc                                                                                           >  4:176292/2‑141 (MQ=255)
                                                           ggcaTCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTc                                                                                           <  3:176292/140‑1 (MQ=255)
                                                                  caAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAAtt                                                                                                                        <  2:253833/105‑1 (MQ=255)
                                                                  caAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGAACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAAtt                                                                                                                        >  1:253833/1‑105 (MQ=255)
                                                                   aaaaAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCAATTAGCGCAGAGAAATTTTCCATCCCGGTACAAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGAAAAGTTTGAATCTTTCGGCATACGGATATGCt                                                                           <  7:347247/149‑1 (MQ=255)
                                                                   aaaCAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCt                                                                           >  7:164012/1‑149 (MQ=255)
                                                                   aaaCAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCt                                                                           >  8:323600/1‑149 (MQ=255)
                                                                      cAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTcca                                                                        <  2:193914/149‑1 (MQ=255)
                                                                       aGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATTATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTccac                                                                       >  8:14608/1‑149 (MQ=255)
                                                                       aGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATTATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTccac                                                                       >  8:20444/1‑149 (MQ=255)
                                                                            tcGCGCGATTTTGATGTAGTCATAACGCCCTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGtt                                                                  >  2:97071/1‑149 (MQ=255)
                                                                             cgcgcgATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTc                                                                 >  3:326444/1‑149 (MQ=255)
                                                                              gcgcgATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGtt                                                                  <  7:323600/147‑1 (MQ=255)
                                                                              gcgcgATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTc                                                                 >  6:376014/1‑148 (MQ=255)
                                                                              gcgcgATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTc                                                                 >  8:337618/1‑148 (MQ=255)
                                                                              gcgcgATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCg                                                                <  1:351348/149‑1 (MQ=255)
                                                                                gcgATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGaa                                                              <  5:376014/149‑1 (MQ=255)
                                                                                 cgATTTTGATGTAGTCAAAACGCACTGCCCTTACCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGaaa                                                             <  6:183124/149‑1 (MQ=255)
                                                                                              gTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGa                                                 <  5:82546/148‑1 (MQ=255)
                                                                                                 aTAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGa                                              <  6:196173/148‑1 (MQ=255)
                                                                                                 aTAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGAc                                             >  8:268899/1‑149 (MQ=255)
                                                                                                 aTAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGAc                                             >  4:139580/1‑149 (MQ=255)
                                                                                                         cttcacttAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATc                                     >  5:213089/1‑149 (MQ=255)
                                                                                                           tcacttAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTg                                   >  3:373031/1‑149 (MQ=255)
                                                                                                                 aGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCa                              <  7:9370/148‑1 (MQ=255)
                                                                                                                        agaAATTTGCCATCCCGGTACCAAAATTATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGAttttt                      <  7:20444/149‑1 (MQ=255)
                                                                                                                         gaAATTTGCCATCCCGGTACCAAAATTATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGAttttt                      <  7:14608/148‑1 (MQ=255)
                                                                                                                            aTTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGtt                  >  7:257379/1‑149 (MQ=255)
                                                                                                                                  cATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGAc            >  3:269130/1‑149 (MQ=255)
                                                                                                                                   aTCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGgcg                                <  4:376686/128‑1 (MQ=255)
                                                                                                                                   aTCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGgcg                                >  3:376686/1‑128 (MQ=255)
                                                                                                                                    tCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGACg           >  5:223578/1‑148 (MQ=255)
                                                                                                                                        ggTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGACGCAggg      >  7:393132/1‑149 (MQ=255)
                                                                                                                                             cAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGAAGGATTTTTGGTTGCTGACGCAGGGCGAt  >  2:3539/1‑148 (MQ=255)
                                                                                                                                              aaaaTCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGCCGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGACGCAGGGCGAt  >  5:286007/1‑147 (MQ=255)
                                                                                                                                              aaaaTCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGACGCAGGGCGAt  >  8:17167/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
AGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGACGCAGGGCGAT  >  NC_000913/3678566‑3678854

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