Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,863,159 A→G 5.0% P275P (CCA→CCG)  ygbK → FliA‑regulated DUF1537 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,863,1590AG 5.0% 248.5 / 2.5 80P275P (CCA→CCG) ygbKFliA‑regulated DUF1537 family protein
Reads supporting (aligned to +/- strand):  ref base A (40/35);  new base G (2/2);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

ATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGTTTTTGCCACCGCCAGCACTGACGCATTGGCAGCAATT  >  NC_000913/2863013‑2863306
                                                                                                                                                  |                                                                                                                                                   
aTTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCCTGCTCTCAAATGACCACCCGCCAGGTAGCACATGACCGTCAAATTGCACCAg                                                                                                                                                    >  5:292895/1‑148 (MQ=255)
aTTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAg                                                                                                                                                    <  6:306225/148‑1 (MQ=255)
aTTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCCGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAg                                                                                                                                                    >  7:224371/1‑148 (MQ=255)
      ctggccCGGCAGTGGGCGCAAGAAAACGGTAAACAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGt                                                                                                                                               >  6:213595/1‑147 (MQ=255)
        ggccCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTg                                                                                                                                              >  1:80110/1‑146 (MQ=255)
         gccCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCTAATTGCACCAGCCCGTGaa                                                                                                                                            <  3:382237/147‑1 (MQ=255)
            cGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTg                                                                                                                                        <  3:317879/148‑1 (MQ=255)
               cAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGAtg                                                                                                                                     <  1:372002/148‑1 (MQ=255)
                   gggCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACAATACCGTCATATTGCACCAGCCCGTGAAGTTGATGTgg                                                                                                                                  >  1:107203/1‑147 (MQ=255)
                       gcAAGAAAACGGTAACCAGGCGCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGc                                                                                                                             <  6:110850/148‑1 (MQ=255)
                         aaGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTg                                                                                                                           <  7:137002/148‑1 (MQ=255)
                         aaGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTg                                                                                                                           >  3:74233/1‑148 (MQ=255)
                          aGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGc                                                                                                                          >  1:464562/1‑148 (MQ=255)
                          aGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGc                                                                                                                          >  1:467358/1‑148 (MQ=255)
                           gAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGcc                                                                                                                         >  3:155418/1‑148 (MQ=255)
                                cGGTAACCAGGCTGGCAAAGCAGGGCTTCCGCTCTATGGGCTCGGCGTAGTGCTCTCCTGTTCATGCTCTCAAATGACAAAACGCCAGGTAGCACATTACCGTCAAATTGCACCAGACCGTGAAGTTTATTTGGCACGCTGCCTCTCa                                                                                                                    <  4:56955/148‑1 (MQ=255)
                                cGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCa                                                                                                                    <  5:213595/148‑1 (MQ=255)
                                       cAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGaaa                                                                                                             <  1:191256/148‑1 (MQ=255)
                                       cAGGCACGCAAAGCAGGGCTGCCGCTCGCTGGGCGCGGCGTAGTGCGCTCCGGGTCATGCGCTCAAATGACCAACCGCCAGGTAGCACATAACCGTCAAATTGCACCAGCCCGTGAAGGTGATGTGGCACGCTGCCTCTCAATTGaaa                                                                                                             <  8:407775/148‑1 (MQ=255)
                                           cTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGTAGTTGATGTGGCACGCTGCCTCTCAATTGAAActct                                                                                                         >  3:9083/1‑148 (MQ=255)
                                            tCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTg                                                                                                        >  1:123713/1‑148 (MQ=255)
                                              gCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCGGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATGGAAACTCTGGc                                                                                                      >  8:478106/1‑148 (MQ=255)
                                              gCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCGGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAAGGGAAACTCTGGc                                                                                                      >  8:471694/1‑148 (MQ=255)
                                               cAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGcc                                                                                                     >  2:374263/1‑148 (MQ=255)
                                               cAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTGCATGCTCTCAAATGACCGCCCGCCCGGGAGCACAGTACCGTCACATGGCACCCGCCCGGGAAGGTGATGTGGCCCGCGGCCGCTCAATTGAAACTCTGGcc                                                                                                     >  1:36612/1‑148 (MQ=255)
                                                aaaGCTCGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCCTGATCTCAAAAGACCAACCGCCAGGTAGCACATAACCGTCAAATTGCACCAGCCCGTGAAGTTGATGGGGCACGCTGCCTCTCAATTGAAACTCTGGCCg                                                                                                    <  6:366027/148‑1 (MQ=255)
                                                aaaGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCg                                                                                                    <  8:54660/148‑1 (MQ=255)
                                                     aGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTAt                                                                                               <  1:374263/148‑1 (MQ=255)
                                                     aGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGCCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTAt                                                                                               <  1:243039/148‑1 (MQ=255)
                                                            ccGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCacac                                                                                         >  6:60804/1‑147 (MQ=255)
                                                                   ctGGGCGCGGCGTAGTGCTCTCAGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTTGCCGCTTATGCACACGAACTGGc                                                                                 <  8:224371/148‑1 (MQ=255)
                                                                     gggCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTgg                                                                                  >  4:57127/1‑145 (MQ=255)
                                                                     gggCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCag                                                                               <  7:193027/148‑1 (MQ=255)
                                                                              gTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTc                                                                      <  2:80110/148‑1 (MQ=255)
                                                                                aGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTg                                                                    >  8:450247/1‑148 (MQ=255)
                                                                                aGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTg                                                                    >  7:222342/1‑148 (MQ=255)
                                                                                   gCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTTCACCAGCCCGTGAAGTTGa                                                                                                                                       >  6:205080/1‑78 (MQ=255)
                                                                                   gCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGa                                                                                                                                       <  5:205080/78‑1 (MQ=255)
                                                                                    ctatcCGGTTCAGGCGCTCAACTGCCCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGAGGTGGCCAGCTGCCTATCAATTGAAACTCTGGCCGCTTATGCACAAGAACTGGCAGAGTGGGTTCTGGGcc                                                                <  2:36612/145‑1 (MQ=255)
                                                                                        ccGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAgg                                                             >  4:278543/1‑147 (MQ=255)
                                                                                          ggTTCATGCTCTCAAATGACCAACCGCCATGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGaaa                                                          <  8:222342/148‑1 (MQ=255)
                                                                                           gTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAg                                                         <  3:1142/148‑1 (MQ=255)
                                                                                               aTGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTa                                                     <  7:192969/148‑1 (MQ=255)
                                                                                               aTGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTa                                                     <  5:60804/148‑1 (MQ=255)
                                                                                               aTGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTa                                                     <  3:278543/148‑1 (MQ=255)
                                                                                                  ctctCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACt                                                   <  3:432359/147‑1 (MQ=255)
                                                                                                    ctcAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTg                                                                    <  4:185859/128‑1 (MQ=255)
                                                                                                    ctcAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTg                                                                    >  3:185859/1‑128 (MQ=255)
                                                                                                     tcAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTg                                                 <  2:406109/146‑1 (MQ=255)
                                                                                                     tcAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCt                                               >  6:281362/1‑148 (MQ=255)
                                                                                                      cAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTc                                              >  5:50809/1‑148 (MQ=255)
                                                                                                        aaTGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCa                                            >  3:157756/1‑148 (MQ=255)
                                                                                                         aTGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCAc                                           >  4:187647/1‑148 (MQ=255)
                                                                                                            aCCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTg                                         <  3:180885/147‑1 (MQ=255)
                                                                                                                 ccGCCAGGTAGCACATTACCGTCAAATTGCACCGGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGttttt                                   <  7:478106/148‑1 (MQ=255)
                                                                                                                 ccGCCAGGTAGCACATTACCGTCAAATTCCACCGGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAACCTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGttttt                                   <  7:471694/148‑1 (MQ=255)
                                                                                                                  cGCCAGGTAGCACATTACAGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCATCTCAATTGAAACTCTGGCCGCTTATGCAAACGAACTGGCAGAGTTGGTTCt                                                                     >  3:186633/1‑113 (MQ=255)
                                                                                                                   gCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTGCTGGGCCAGGATAGGGTACTTGCTCCACTGGTTTTTGc                                 >  8:184919/1‑148 (MQ=255)
                                                                                                                     cAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGTTTTTGCCa                               >  3:198119/1‑148 (MQ=255)
                                                                                                                       ggTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGTTTTTGCCAcc                             >  5:290137/1‑148 (MQ=255)
                                                                                                                          aGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGTTTTTGCCAcctc                           >  4:156276/1‑145 (MQ=255)
                                                                                                                             acaTTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGTTTTTGCCACCGCCAGc                       <  6:50809/148‑1 (MQ=255)
                                                                                                                               aTTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGTTTTTGCCACCGCCAGCAc                     >  8:394443/1‑148 (MQ=255)
                                                                                                                                ttACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGTTTTTGCCACCGCCAGCACt                    >  7:466463/1‑148 (MQ=255)
                                                                                                                                 tACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGtt                                      <  8:181366/129‑1 (MQ=255)
                                                                                                                                 tACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGtt                                      >  7:181366/1‑129 (MQ=255)
                                                                                                                                 aaCCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGATAGTGTACTTGCTCCACTGGTTTTTGCCACCGCCAGCACt                    >  3:398277/2‑147 (MQ=255)
                                                                                                                                  aCCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGTTTTTGCCACCGCCAg                        <  5:42937/142‑1 (MQ=255)
                                                                                                                                  aCCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGTTTTTGCCACCGCCAg                        >  6:42937/1‑142 (MQ=255)
                                                                                                                                    cGTCAAATTGCCCCAGCCCGTGAAGTTCACGTGGCACGCTGCCCCTCAATTGAAACTCTGGGCGCTAATGCACACGAACTGGCATAGTGGGTTCTGGGCCAGGAAAGTTTACTTGCTCCACTGGTTTTTGCCACCGCCAGCACTGACg                <  2:107203/148‑1 (MQ=255)
                                                                                                                                    cGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGaa                                                           >  7:462108/1‑105 (MQ=255)
                                                                                                                                    cGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGaa                                                           <  8:462108/105‑1 (MQ=255)
                                                                                                                                     gTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGTTTTTGCCACCGCCAGCACTGACg                >  3:216786/1‑147 (MQ=255)
                                                                                                                                     gTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGTTTTTGCCACCGCCAGCACTGACg                >  8:470585/1‑147 (MQ=255)
                                                                                                                                     gTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCAATGAAAGTTTACTTGCTCCACTGGTTTTTGCCACCGCCAGCACTGACg                >  6:252394/1‑147 (MQ=255)
                                                                                                                                      tCAAATTGCACCATCCCGTGAAGTTGATGTGGCACGATGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACT‑GCAGAGTTGGTTCTGGGCCAGGAAAGTGTACATGCTACACAGGTTTTTGCCACCGCCACCACTGACGCAt             >  2:440782/1‑148 (MQ=255)
                                                                                                                                      tCAAATTGCACCAGCCCGTGTAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGATAGTGTACTTGCTCCACTGGTTTTTGCCACCGCCAGCACTGACGCa              <  4:9083/148‑1 (MQ=255)
                                                                                                                                      tCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGTTTTTGCCACCGCCAGCACTGACGCa              >  8:382992/1‑148 (MQ=255)
                                                                                                                                      tCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGTTTTTGCCACCGCCAGCACTGACGCa              >  6:270459/1‑148 (MQ=255)
                                                                                                                                        aaaTTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGTTTTTGCCACCGCCAGCACTGACGCAt             <  2:123713/147‑1 (MQ=255)
                                                                                                                                                  aGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGTTTTTGCCACCGCCAGCACTGACGCATTGGCAGCAAtt  <  7:382992/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTGGTTTTTGCCACCGCCAGCACTGACGCATTGGCAGCAATT  >  NC_000913/2863013‑2863306

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