Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,628,011 T→A 5.1% N310Y (AAT→TAT)  yfgF ← cyclic‑di‑GMP phosphodiesterase, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,628,0110TA 5.1% 252.3 / 3.5 78N310Y (AAT→TAT) yfgFcyclic‑di‑GMP phosphodiesterase, anaerobic
Reads supporting (aligned to +/- strand):  ref base T (32/42);  new base A (2/2);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

AATAACAAAGCGCAGACCAGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGATCATCAACACGACCGCCCAGAGCAGCGAAATCAGCTTATAA  >  NC_000913/2627864‑2628153
                                                                                                                                                   |                                                                                                                                              
aaTAACAAAGCGCAGACCAGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAAtt                                                                                                                                                <  7:287312/148‑1 (MQ=255)
 aTAACAAAGCGCAGACCAGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGa                                                                                                                                              <  4:119897/149‑1 (MQ=255)
  tAACAAAGCGCAGACCAGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGa                                                                                                                                              >  6:79233/1‑148 (MQ=255)
    aCAAAGCGCAGACCAGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACaa                                                                                                                                           <  8:195804/149‑1 (MQ=255)
     cAAAGCGCAGACCAGGGGGCATCACGCACCGCGCGATTCAGGGCGCCAACATTTGGCACATGAACCCCCGCGTCCACAAACGCAAGCCGCAGCAAGCCTCTGACTACCGCTCTCTGCCCGGTAGCCAGTACAGCCATAAAATTGACaa                                                                                                                                           >  1:185113/1‑148 (MQ=255)
        aGCGCAGACCAGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATaaa                                                                                                                                       <  4:55434/149‑1 (MQ=255)
           gcAGACCAGGGGGCATCACGCAACGCGCGTTTCAGGTCGCGATCATTTGGCAGATGAACGACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCGAGTACAGCGATATAATTGACAATtaaata                                                                                                                                    <  2:185113/149‑3 (MQ=255)
             aGACCAGGGGGCATCACGCAACGCGCGATTCCGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAgaga                                                                                                                                  <  5:40353/149‑1 (MQ=255)
              gACCAGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGaa                                                                                                                                 <  1:295680/149‑1 (MQ=255)
                  aGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATAcc                                                                                                                             <  5:177470/149‑1 (MQ=255)
                  aGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATAcc                                                                                                                             <  1:223991/149‑1 (MQ=255)
                   gggggCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCTCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCa                                                                                                                            <  7:280180/149‑1 (MQ=255)
                   gggggCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCa                                                                                                                            <  4:309402/149‑1 (MQ=255)
                   gggggCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCa                                                                                                                            <  4:171591/149‑1 (MQ=255)
                     gggCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGa                                                                                                                          <  5:79233/149‑1 (MQ=255)
                            aCGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTgg                                                                                                                   <  4:306255/149‑1 (MQ=255)
                            aCGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTg                                                                                                                    <  3:35374/148‑1 (MQ=255)
                                aaCGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTggag                                                                                                                 >  5:202295/1‑147 (MQ=255)
                                   gcgcgATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAAt                                                                                                                                          >  7:261529/1‑119 (MQ=255)
                                   gcgcgATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAAt                                                                                                                                          <  8:261529/119‑1 (MQ=255)
                                    cgcgATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGtt                                                                                                           <  3:215680/149‑1 (MQ=255)
                                           cAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTc                                                                                                     <  8:130118/148‑1 (MQ=255)
                                            aGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAg                                                                                                   >  6:169946/1‑149 (MQ=255)
                                                  cgcaCATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGCCGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACGGGAGCAGGTGATGGTCAGCTGCgg                                                                                             >  2:303418/4‑148 (MQ=255)
                                                   gAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTg                                                                                            <  5:257145/149‑1 (MQ=255)
                                                     aCATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTgg                                                                                           >  1:34368/1‑148 (MQ=255)
                                                          tGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAg                                                                                      <  2:211987/148‑1 (MQ=255)
                                                            gCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGcc                                                                                    <  2:32389/148‑1 (MQ=255)
                                                            gCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGcc                                                                                    <  8:205489/148‑1 (MQ=255)
                                                                   aaCCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATaaa                                                                            >  1:273982/1‑149 (MQ=255)
                                                                        ccGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGa                                                                       <  6:56289/149‑1 (MQ=255)
                                                                        ccGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCGGACGGGGTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGGTATGGGCAGCTGCGTGGTATAGCCAGGATAAATGGGa                                                                       <  1:303418/149‑1 (MQ=255)
                                                                          gggTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAAt                                                                     <  2:198637/149‑1 (MQ=255)
                                                                              ccACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTaa                                                                 >  7:258423/1‑149 (MQ=255)
                                                                                acaTACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAAtt                                                               >  4:315733/1‑149 (MQ=255)
                                                                                acaTACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAAtt                                                               >  7:224325/1‑149 (MQ=255)
                                                                                 caTACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAAttt                                                              >  4:160961/1‑149 (MQ=255)
                                                                                  aTACGCAAGCCGCTGAATGCGTCTGACTCCCGCTCTCTGACTGGTTGCCAGTACAGCCATATAATTGACAATAAACGAGCATACCCGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAAtttt                                                             <  2:289321/149‑1 (MQ=255)
                                                                                  aTACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAAtttt                                                             >  1:107672/1‑149 (MQ=255)
                                                                                    aCGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGTATGTAATTTTGa                                                           >  8:307972/1‑149 (MQ=255)
                                                                                      gCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGa                                                           >  8:125782/1‑147 (MQ=255)
                                                                                       cAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAg                                                        <  5:169946/149‑1 (MQ=255)
                                                                                       cAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAg                                                        >  1:187796/1‑149 (MQ=255)
                                                                                         agacGCTTAATGAGTCTTACTACAGATCTCTGAAGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAAATGGAGTAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCt                                                      <  3:174004/146‑1 (MQ=255)
                                                                                             gCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTgg                                                  <  8:32188/149‑1 (MQ=255)
                                                                                                 aaTGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGc                                              >  6:329866/1‑149 (MQ=255)
                                                                                                     cGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGa                                           <  8:258423/148‑1 (MQ=255)
                                                                                                     cGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGa                                           <  7:307972/148‑1 (MQ=255)
                                                                                                     cGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGa                                           <  3:160961/148‑1 (MQ=255)
                                                                                                     cGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGGTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGat                                          <  6:202295/149‑1 (MQ=255)
                                                                                                       tCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGTTATAGCCAGGATAAATGGGAATGTAATTTTGATATCTGTGGATGCTGatca                                        >  3:19061/1‑149 (MQ=255)
                                                                                                       tCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGatca                                        >  6:129958/1‑149 (MQ=255)
                                                                                                        cTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGatcat                                       >  4:130478/1‑149 (MQ=255)
                                                                                                         tGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAATTGGGAATGTAAt                                                                >  4:7309/1‑123 (MQ=255)
                                                                                                         tGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAAt                                                                <  3:7309/123‑1 (MQ=255)
                                                                                                         tGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGatcat                                       >  6:255161/1‑148 (MQ=255)
                                                                                                           aCTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGatcatc                                      <  5:255161/147‑1 (MQ=255)
                                                                                                            cTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGATCATCAac                                   <  2:34368/149‑1 (MQ=255)
                                                                                                             tACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATAGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGcc                                                                                    <  7:202216/99‑1 (MQ=255)
                                                                                                             tACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATAGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGcc                                                                                    >  8:202216/1‑99 (MQ=255)
                                                                                                             tACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATAGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGATCATCAaca                                  >  8:179923/1‑149 (MQ=255)
                                                                                                                 gCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGGGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGATCATCAACACGAc                              >  2:37773/1‑149 (MQ=255)
                                                                                                                       tGACGGGTTGCCAGTACAGCCATATAATTGACAATaaaa                                                                                                                                      >  6:81149/1‑39 (MQ=255)
                                                                                                                       tGACGGGTTGCCAGTACAGCCATATAATTGACAATaaaa                                                                                                                                      <  5:81149/39‑1 (MQ=255)
                                                                                                                       tGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGATCATCAACACGACCGcc                          >  1:253625/1‑147 (MQ=255)
                                                                                                                       tGACGGGTTGCCAGTACAGCCATATAATTGACAAAAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCATGTGGGTGGTATAGCCAGTATAAATGTGAATGTAATTTTGATAGCAGTGGATGGTGATCAACAACAGGCCCGcg                          >  5:308573/1‑146 (MQ=255)
                                                                                                                            ggTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGATCATCAACACGACCGCCCAGagca                   <  8:51055/149‑1 (MQ=255)
                                                                                                                            ggTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGATCATCAACACGACCGCCCAGagca                   <  5:329866/149‑1 (MQ=255)
                                                                                                                                gCCAGTACAGCCATATAATAGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGATCATCAACACGACCGCCCAGAGCAGCg                <  7:179923/148‑1 (MQ=255)
                                                                                                                                  cAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGATCATCAacac                                 <  6:300407/129‑1 (MQ=255)
                                                                                                                                  cAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGATCATCAacac                                 >  5:300407/1‑129 (MQ=255)
                                                                                                                                       cAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGATCATCAACACGACCGCCCAGAGCAGCGAAATCAGc        >  8:118989/1‑149 (MQ=255)
                                                                                                                                       cAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGACAGTATAAATGGGAATGTAATTTTGATATCTGTTGATGCTGa                                           >  6:83630/1‑114 (MQ=255)
                                                                                                                                       aatccATATAATTTACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATTGTCAGCTGCGTTGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATTCTGa                                           <  5:83630/111‑1 (MQ=255)
                                                                                                                                           cATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGATCATCAACACGACCGCCCAGAGCAGCGAAATCAGCTtat    <  2:253625/149‑1 (MQ=255)
                                                                                                                                            atatAATTGACAATAAAAGAGAATACCAGATAACTGGTGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGATGATGATCCACACGACCGCCCAGAGCAGCGAAATCAGATtat    >  5:250199/1‑148 (MQ=255)
                                                                                                                                              ataATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGATGTGGATGCTGATGATCAACACGAACGCCCAGAGAAGCGAAATCAGCTTATaa  >  3:47766/1‑148 (MQ=255)
                                                                                                                                                aaTTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAAttt                                                              >  2:42580/1‑86 (MQ=255)
                                                                                                                                                aaTTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAAttt                                                              <  1:42580/86‑1 (MQ=255)
                                                                                                                                                  ttGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGATCATCAACACGACCGCCCAGAGCAGCGAAATCa          >  8:15175/1‑136 (MQ=255)
                                                                                                                                                  ttGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGATCATCAACACGACCGCCCAGAGCAGCGAAATCa          <  7:15175/136‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
AATAACAAAGCGCAGACCAGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGAATACCAGATAACTGGAGGAGGTTATGGTCAGCTGCGTGGTATAGCCAGGATAAATGGGAATGTAATTTTGATAGCTGTGGATGCTGATCATCAACACGACCGCCCAGAGCAGCGAAATCAGCTTATAA  >  NC_000913/2627864‑2628153

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