Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,346,451 C→T 6.0% D98N (GAC→AAC)  gmr ← cyclic‑di‑GMP phosphodiesterase; csgD regulator; modulator of RNase II stability

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,346,4510CT 6.0% 212.8 / 3.2 67D98N (GAC→AAC) gmrcyclic‑di‑GMP phosphodiesterase; csgD regulator; modulator of RNase II stability
Reads supporting (aligned to +/- strand):  ref base C (27/36);  new base T (2/2);  total (29/38)
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.00e-01

TTCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTC  >  NC_000913/1346303‑1346598
                                                                                                                                                    |                                                                                                                                                   
ttCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTc                                                                                                                                                     <  4:144747/149‑1 (MQ=255)
   aaCCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCAGTGTGTcatg                                                                                                                                                  <  4:156409/149‑1 (MQ=255)
       cTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTcatgccat                                                                                                                                              <  2:410814/149‑1 (MQ=255)
            tAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATGAATCACGAGCACCCGCCGGTGTGTCATGCCATGCAGa                                                                                                                                         >  1:551373/1‑149 (MQ=255)
                cTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGAcgcg                                                                                                                                     >  4:518937/1‑149 (MQ=255)
                       tAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAg                                                                                                                              >  4:562782/1‑149 (MQ=255)
                         aCCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGcc                                                                                                                            >  2:487005/1‑149 (MQ=255)
                                 aaGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCg                                                                                                                    >  1:110452/1‑149 (MQ=255)
                                  aGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGt                                                                                                                   <  1:538004/149‑1 (MQ=255)
                                      ttGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTTATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTAt                                                                                                               >  4:200341/1‑149 (MQ=255)
                                      ttGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTTATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTAt                                                                                                               >  2:39614/1‑149 (MQ=255)
                                         gATATTCCCCCGGCTGTCGAGAATCACTATTTCGGAATGTGCTTCAGAAACTACCTGCTCGGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGTGTTGCCAGCCCATTCCCGTTTATTGa                                                                                                            >  2:516099/1‑149 (MQ=255)
                                         gATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTTATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGa                                                                                                            <  3:200341/149‑1 (MQ=255)
                                         gATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTTATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGa                                                                                                            <  1:39614/149‑1 (MQ=255)
                                         gATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGa                                                                                                            <  2:82400/149‑1 (MQ=255)
                                            atTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGAttt                                                                                                         >  1:452818/1‑149 (MQ=255)
                                             ttCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGc                                                                                                                             <  1:517857/128‑1 (MQ=255)
                                             ttCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGc                                                                                                                             >  2:517857/1‑128 (MQ=255)
                                             ttCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGAtttt                                                                                                        <  4:28308/149‑1 (MQ=255)
                                              tCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGAttttt                                                                                                       >  1:144437/1‑149 (MQ=255)
                                               cccccGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTc                                                                                                      <  1:532802/149‑1 (MQ=255)
                                                ccccGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCg                                                                                                     >  2:58488/1‑149 (MQ=255)
                                                 cccGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGt                                                                                                    >  4:264531/1‑149 (MQ=255)
                                                  ccGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTc                                                                                                   >  4:179248/1‑149 (MQ=255)
                                                     gCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTa                                                                                                <  4:289984/149‑1 (MQ=255)
                                                               aTCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGcac                                                                                      <  1:439328/149‑1 (MQ=255)
                                                                   cTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGaa                                                                                  <  2:144437/149‑1 (MQ=255)
                                                                    tATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGaaa                                                                                 >  2:445793/1‑149 (MQ=255)
                                                                    tATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGaaa                                                                                 >  4:207458/1‑149 (MQ=255)
                                                                    tATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGAACAGaaa                                                                                 <  2:313960/149‑1 (MQ=255)
                                                                      ttGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGa                                                                               <  3:414702/149‑1 (MQ=255)
                                                                         cGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCg                                                                            <  3:469004/149‑1 (MQ=255)
                                                                         cGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCg                                                                            <  3:273383/149‑1 (MQ=255)
                                                                           gAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGa                                                                                           >  3:289685/1‑132 (MQ=255)
                                                                           gAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGa                                                                                           <  4:411317/132‑1 (MQ=255)
                                                                           gAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGa                                                                                           >  3:411317/1‑132 (MQ=255)
                                                                           gAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGa                                                                                           <  4:289685/132‑1 (MQ=255)
                                                                            aaTGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCtt                                                                         <  2:315748/149‑1 (MQ=255)
                                                                             aTGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCttt                                                                        <  1:357835/149‑1 (MQ=255)
                                                                                   ttCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGa                                                                  >  4:263753/1‑149 (MQ=255)
                                                                                   ttCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGa                                                                  >  2:315692/1‑149 (MQ=255)
                                                                                        aaaCTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTc                                                             <  4:310442/149‑1 (MQ=255)
                                                                                           cTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGa                                                          <  2:110452/149‑1 (MQ=255)
                                                                                                tGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTc                                                                                                   <  1:97913/103‑1 (MQ=255)
                                                                                                tGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTc                                                                                                   >  2:97913/1‑103 (MQ=255)
                                                                                                 gCTCAGCAAAGGATACCCCGTGTGATAAATCACGAGCAACCGCTGGTGTGTCCTGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGAGCGCTTTCATCGACGGTCAGACtcatc                                                    <  2:385912/149‑1 (MQ=255)
                                                                                                    cAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACtcatcatc                                                 <  4:564437/149‑1 (MQ=255)
                                                                                                       cAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAgg                                              >  4:364622/1‑149 (MQ=255)
                                                                                                        aaaGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGc                                             >  4:30609/1‑149 (MQ=255)
                                                                                                               aGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGt                                      <  1:491057/149‑1 (MQ=255)
                                                                                                               aGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGt                                      <  4:562195/149‑1 (MQ=255)
                                                                                                                 cccGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGtga                                    <  3:263753/149‑1 (MQ=255)
                                                                                                                    gtgtGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGTCGGTCAGACTCATCATCAGGCTGGAGGttatga                                 <  4:777/149‑1 (MQ=255)
                                                                                                                         aTAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACtcat                                                     <  3:310672/124‑1 (MQ=255)
                                                                                                                         aTAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACtcat                                                     >  4:310672/1‑124 (MQ=255)
                                                                                                                             aTCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGTGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCAttt                        <  1:516099/149‑1 (MQ=255)
                                                                                                                              tCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTc                       >  3:62213/1‑149 (MQ=255)
                                                                                                                              tCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTc                       >  1:211850/1‑149 (MQ=255)
                                                                                                                                aCGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGTCTGGCGGTGATGACCGTCATTTccc                     <  4:405911/149‑1 (MQ=255)
                                                                                                                                aCGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTccc                     <  4:373738/149‑1 (MQ=255)
                                                                                                                                aCGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTccc                     <  3:30609/149‑1 (MQ=255)
                                                                                                                                   aGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGa                  <  1:445793/149‑1 (MQ=255)
                                                                                                                                    gCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGATGGTc                                                             >  3:530978/1‑105 (MQ=255)
                                                                                                                                    gCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGATGGTc                                                             <  4:530978/105‑1 (MQ=255)
                                                                                                                                        ccgccgGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGa             <  1:103248/149‑1 (MQ=255)
                                                                                                                                        ccgccgGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCATCATTTCCCTGATGCGa             <  3:207458/149‑1 (MQ=255)
                                                                                                                                           ccgGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGagta          <  2:240858/149‑3 (MQ=255)
                                                                                                                                           ccgGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGagca          <  3:364622/149‑1 (MQ=255)
                                                                                                                                              gtgtgtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGcc       >  4:33153/1‑149 (MQ=255)
                                                                                                                                                   tCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTc  >  3:489537/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTC  >  NC_000913/1346303‑1346598

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑