Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,684,506 Δ1 bp 5.5% coding (1113/1989 nt) yhjK ← cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,684,5060G. 5.5% 254.8 / 8.4 73coding (1113/1989 nt)yhjKcyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (32/37);  new base . (2/2);  total (34/39)
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.25e-01

CAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGA  >  NC_000913/3684359‑3684651
                                                                                                                                                   |                                                                                                                                                 
cAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTg                                                                                                                                                   <  5:90489/148‑1 (MQ=255)
  gTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGCTCGGCGGTgag                                                                                                                                                 <  7:253332/148‑1 (MQ=255)
      aTTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCg                                                                                                                                             >  5:71209/1‑148 (MQ=255)
       ttAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTCCGGCGGTGAGATCGc                                                                                                                                            <  6:218312/148‑1 (MQ=255)
               atCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGaa                                                                                                                                     <  4:287288/147‑1 (MQ=255)
               atCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGaa                                                                                                                                     >  2:126895/1‑147 (MQ=255)
                tCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAAc                                                                                                                                    <  6:7738/147‑1 (MQ=255)
                    tctctTCCGTCAACCGCTTCTGGGCGGCTTCCATATGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGATCGGCGGTGAGATCGCCGTAGAACATcgc                                                                                                                               >  2:90606/1‑148 (MQ=255)
                     ctctTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATcgc                                                                                                                               >  3:52281/1‑147 (MQ=255)
                        ttCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATcgccacg                                                                                                                           >  2:122719/1‑148 (MQ=255)
                        ttCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATcgcc                                                                                                                              <  4:62840/145‑1 (MQ=255)
                         tCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTgaga                                                                                                                                                >  8:48497/1‑126 (MQ=255)
                         tCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTgaga                                                                                                                                                <  7:48497/126‑1 (MQ=255)
                             tCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAAt                                                                                                                      <  6:6674/148‑1 (MQ=255)
                              cAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAAt                                                                                                                      >  1:157007/1‑147 (MQ=255)
                                 ccGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTa                                                                                                                  <  6:308619/148‑1 (MQ=255)
                                  cGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCaa                                                                                                                       <  7:270293/142‑1 (MQ=255)
                                  cGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTAc                                                                                                                 <  3:259747/148‑1 (MQ=255)
                                  cGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTAc                                                                                                                 <  2:157007/148‑1 (MQ=255)
                                  cGCTTCTGGGCGGCTTCCATCTGCTGCCGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCaa                                                                                                                       >  8:270293/1‑142 (MQ=255)
                                      tCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTg                                                                                                             >  6:10684/1‑148 (MQ=255)
                                          ggcggcTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGcgga                                                                                                         <  4:52281/148‑1 (MQ=255)
                                          ggcggcTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGc                                                                                                            <  5:10684/145‑1 (MQ=255)
                                          ggcggcTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGc                                                                                                            <  6:294696/145‑1 (MQ=255)
                                             ggcTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGcggacg                                                                                                       >  3:181745/1‑147 (MQ=255)
                                                 tCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATcgcc                                                                                                                              >  6:49844/1‑120 (MQ=255)
                                                 tCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATcgcc                                                                                                                              <  5:49844/120‑1 (MQ=255)
                                                  ccATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACCGTGCGGACGGAGtt                                                                                                  >  5:5001/1‑147 (MQ=255)
                                                  ccATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGt                                                                                                   <  6:23920/146‑1 (MQ=255)
                                                  ccATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGt                                                                                                   <  1:126895/146‑1 (MQ=255)
                                                   cATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTg                                                                                                 <  1:122719/147‑1 (MQ=255)
                                                   cATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGaa                                                                                               >  4:243320/1‑148 (MQ=255)
                                                   cATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGaa                                                                                               >  4:239311/1‑148 (MQ=255)
                                                    aTCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGaa                                                                                               <  8:32683/148‑1 (MQ=255)
                                                    aTCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGaa                                                                                               >  6:209184/1‑148 (MQ=255)
                                                        gctgcGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACg                                                                                           >  5:71537/1‑148 (MQ=255)
                                                           gcGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGtt                                                                                         >  6:201183/1‑147 (MQ=255)
                                                            cGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGtt                                                                                         <  5:160183/146‑1 (MQ=255)
                                                                tCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCg                                                                                   <  1:162999/148‑1 (MQ=255)
                                                                 cAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTgag                                                                                                                                                 >  8:201182/1‑85 (MQ=255)
                                                                 cAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTgag                                                                                                                                                 <  7:201182/85‑1 (MQ=255)
                                                                  aaaGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGc                                                                                 <  6:71209/148‑1 (MQ=255)
                                                                   aaGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCgg                                                                                  >  5:9421/1‑146 (MQ=255)
                                                                    aGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAg                                                                               >  2:80364/1‑148 (MQ=255)
                                                                     gAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAgg                                                                              >  8:298241/1‑148 (MQ=255)
                                                                      aaCTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGgc                                                                             <  7:74260/148‑1 (MQ=255)
                                                                       aCTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTAc                                                                                                                 <  6:177250/111‑1 (MQ=255)
                                                                       aCTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTAc                                                                                                                 >  5:177250/1‑111 (MQ=255)
                                                                            aTCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACg                                                                                           <  1:254136/128‑1 (MQ=255)
                                                                            aTCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACg                                                                                           >  2:254136/1‑128 (MQ=255)
                                                                                gATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCAt                                                                   >  7:81194/1‑148 (MQ=255)
                                                                                 aTTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCAt                                                                   >  6:109176/1‑147 (MQ=255)
                                                                                 aTTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATg                                                                  <  4:181745/148‑1 (MQ=255)
                                                                                       gCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGgcgc                                                                           >  7:183785/1‑133 (MQ=255)
                                                                                       gCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGgcgc                                                                           <  8:183785/133‑1 (MQ=255)
                                                                                       gCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATa                                                               <  1:80364/145‑1 (MQ=255)
                                                                                       gCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAg                                                              <  6:71537/146‑1 (MQ=255)
                                                                                         cTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAg                                                          >  3:152691/1‑148 (MQ=255)
                                                                                          tttATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAg                                                          >  2:210012/1‑147 (MQ=255)
                                                                                              tGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTg                                                    <  3:243320/148‑1 (MQ=255)
                                                                                              tGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTg                                                    <  3:239311/148‑1 (MQ=255)
                                                                                                          aaTGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAg                                          <  5:1380/147‑1 (MQ=255)
                                                                                                            tGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTg                                       <  8:81194/148‑1 (MQ=255)
                                                                                                              ccGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGAt                                     <  6:9421/148‑1 (MQ=255)
                                                                                                               cGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGAtt                                    >  3:229089/1‑148 (MQ=255)
                                                                                                                    aaaTAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTg                               >  4:188175/1‑148 (MQ=255)
                                                                                                                    aaaTAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTg                               <  5:109176/148‑1 (MQ=255)
                                                                                                                                 gTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTg                   >  8:95879/1‑147 (MQ=255)
                                                                                                                                   aaaGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTa                 >  4:284065/1‑147 (MQ=255)
                                                                                                                                       cTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTg                                                    <  6:78579/108‑1 (MQ=255)
                                                                                                                                       cTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTg                                                    >  5:78579/1‑108 (MQ=255)
                                                                                                                                         gTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGtt          <  1:90606/148‑1 (MQ=255)
                                                                                                                                          ttCGGCGGTGAGATCGCCGTAGCACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGCATACGTTCAATCGGCAGGCGCTCTCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTg         >  6:306109/1‑148 (MQ=255)
                                                                                                                                          ttCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGtt          <  7:298241/147‑1 (MQ=255)
                                                                                                                                          ttCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGt           <  1:210012/146‑1 (MQ=255)
                                                                                                                                                ggTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATg   <  3:188175/148‑1 (MQ=255)
                                                                                                                                                 gTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGa  <  4:152691/148‑1 (MQ=255)
                                                                                                                                                 gTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGa  >  4:192832/1‑148 (MQ=255)
                                                                                                                                                  tGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGa  >  8:20074/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGA  >  NC_000913/3684359‑3684651

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: