Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 856,695 A→G 5.7% T245A (ACC→GCC)  ybiT → ABC‑F family putative regulatory ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913856,6950AG 5.7% 205.1 / 2.2 71T245A (ACC→GCC) ybiTABC‑F family putative regulatory ATPase
Reads supporting (aligned to +/- strand):  ref base A (32/33);  new base G (2/2);  total (36/35)
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.75e-01

TGCTGAACGAGCGTGACAGCACCATGATCATCATCTCGCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACT  >  NC_000913/856552‑856843
                                                                                                                                               |                                                                                                                                                    
tGCTGAACGAGCGTGACAGCACCATGATCATCATCTCGCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAAGTACGTTGAGTACATGACGGCGGCGGCCCAg                                                                                                                                                 >  2:391062/1‑149 (MQ=255)
tGCTGAACGAGCGTGACAGCACCATGATCATCATCTCGCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGGCCCAg                                                                                                                                                 >  2:385717/1‑149 (MQ=255)
     aaCGAGCGTGACAGCACCATGATCATCATCTCGCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGcgcg                                                                                                                                            >  4:273062/1‑149 (MQ=255)
      aCGAGCGTGACAGCACCATGATCATCATCTCGCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGt                                                                                                                                           <  7:27652/149‑1 (MQ=255)
         aGCGTGACAGCACCATGATCATCATCTCGCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGaa                                                                                                                                        <  4:45335/149‑1 (MQ=255)
         aGCGAGACAGCACCATGATCATCATCTCGCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACAGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGaa                                                                                                                                        <  7:290414/149‑1 (MQ=255)
           cGTGACAGCACCATGATCATCATCTCGCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACg                                                                                                                                      >  7:127638/1‑149 (MQ=255)
             tGACAGCACCATGATCATCATCTCGCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTc                                                                                                                                    <  2:105477/149‑1 (MQ=255)
               aCAGCACCATGATCATCATCTCGCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTctg                                                                                                                                  <  5:311560/149‑1 (MQ=255)
                cAGCACCATGATCATCATCTCGCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTctgc                                                                                                                                 <  8:290205/149‑1 (MQ=255)
                     ccATGATCATCATCTCGCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGcc                                                                                                                            >  5:369588/1‑149 (MQ=255)
                      cATGATCATCATCTCGCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGcc                                                                                                                            >  8:359412/1‑148 (MQ=255)
                      cATGATCATCATCTCGCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGATTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGcc                                                                                                                            >  2:324494/1‑148 (MQ=255)
                       aTGATCATCATCTCGCACGACCGTCAGTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGa                                                                                                                          <  7:359412/149‑1 (MQ=255)
                           tcatcatcTCGCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAAc                                                                                                                      <  3:292092/149‑1 (MQ=255)
                           tcatcatcTCGCACGACCGTCACTTCCCTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAAc                                                                                                                      >  3:271832/1‑149 (MQ=255)
                             atcatcTCGCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGc                                                                                                                    >  3:272845/1‑149 (MQ=255)
                              tcatcTCGCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGc                                                                                                                    >  3:7344/1‑148 (MQ=255)
                                atcTCGCACGACCGTCACTTCCCTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCaa                                                                                                                 <  4:271832/149‑1 (MQ=255)
                                     gCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCaagaag                                                                                                             >  4:267710/1‑148 (MQ=255)
                                     gCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCaagaag                                                                                                             >  7:67147/1‑148 (MQ=255)
                                     gCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGACCGATAACGCCaag                                                                                                                <  3:273062/145‑1 (MQ=255)
                                       aCGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAagaaga                                                                                                            >  7:317267/1‑147 (MQ=255)
                                       aCGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGaaa                                                                                                          <  8:127638/149‑1 (MQ=255)
                                         gACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAgc                                                                                                        <  3:267710/149‑1 (MQ=255)
                                                 cTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTg                                                                                                <  6:369588/149‑1 (MQ=255)
                                                    ccTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCt                                                                                              <  1:154810/148‑1 (MQ=255)
                                                     cTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCt                                                                                              >  2:155888/1‑147 (MQ=255)
                                                       tAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGt                                                                                          >  1:226397/1‑149 (MQ=255)
                                                         aCATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTTCATGACGGCGGCGCCCCAGGCGCGTGAACGTCTGCTGGACGATAAAGCCAAGAAGAAAGCGCAGATTGCTGAGTTg                                                                                        >  6:23308/1‑149 (MQ=255)
                                                         aCATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGtt                                                                                         <  1:60025/148‑1 (MQ=255)
                                                           aTGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGc                                                                                       >  5:122242/1‑148 (MQ=255)
                                                            tGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCaa                                                                                     >  7:315808/1‑149 (MQ=255)
                                                              gTCTGTACCCACATGGCGGAGCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAAt                                                                                    <  7:298614/148‑1 (MQ=255)
                                                                    aCCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGt                                                                             >  2:134514/1‑149 (MQ=255)
                                                                     cccACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGt                                                                             >  2:365311/1‑148 (MQ=255)
                                                                      ccACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTctgct                                                                                                                                >  4:335182/1‑96 (MQ=255)
                                                                      ccACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTctgct                                                                                                                                <  3:335182/96‑1 (MQ=255)
                                                                       cacaTGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTc                                                                               <  4:52818/144‑1 (MQ=255)
                                                                       cacaTGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTc                                                                               >  3:52818/1‑144 (MQ=255)
                                                                       cacaTGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAg                                                                          <  5:363214/149‑1 (MQ=255)
                                                                        acaTGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGAGAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGc                                                                         >  1:139258/1‑149 (MQ=255)
                                                                           tGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGc                                                                      >  8:146462/1‑149 (MQ=255)
                                                                                 ttCTGGATTACGGGGAGCTGCGCGTTTATCCGGGTCACTCCGCTGAGTACATGACGGCGGCGACCCAGGCGCGTGACCGGCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAgc                                                                <  1:238748/148‑1 (MQ=255)
                                                                                      gATTACGTCGAGCTGCGCGTTTATCCTGGTAACTACGATGAGTACATGACGGCGGCGGCCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCaa                                                           <  1:391062/149‑1 (MQ=255)
                                                                                      gATTACGGCGCGCTGCGCGTTTATACGGGTAACTACGATGACTACATTACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTa                                                                  <  6:327045/142‑1 (MQ=255)
                                                                                      gATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGGCCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCaa                                                           <  1:385717/149‑1 (MQ=255)
                                                                                      gATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCCGGCGCGTGAACGTCTGATGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTa                                                                  >  5:327045/1‑142 (MQ=255)
                                                                                            ggCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTc                                                     >  5:197193/1‑149 (MQ=255)
                                                                                            ggCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTc                                                     >  7:19112/1‑149 (MQ=255)
                                                                                                   tGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAAtct                                              >  3:126077/1‑149 (MQ=255)
                                                                                                      gcgTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCgg                                           <  4:7344/149‑2 (MQ=255)
                                                                                                        gTTTATCCGTGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCa                                         <  2:76384/149‑1 (MQ=255)
                                                                                                        gTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCa                                         <  4:272845/149‑1 (MQ=255)
                                                                                                        gTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCa                                         <  7:146462/149‑1 (MQ=255)
                                                                                                         tttATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAg                                        <  8:317267/149‑1 (MQ=255)
                                                                                                             tCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCaa                                    >  7:10799/1‑149 (MQ=255)
                                                                                                                      cTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTcgcgcg                           >  2:178207/1‑149 (MQ=255)
                                                                                                                       tACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAagaaga                                                                                                            >  6:73646/1‑67 (MQ=255)
                                                                                                                       tACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAagaaga                                                                                                            <  5:73646/67‑1 (MQ=255)
                                                                                                                       tACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTcgcgcgc                          >  2:339488/1‑149 (MQ=255)
                                                                                                                            tGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCgcgcgcgcgc                      >  5:209401/1‑148 (MQ=255)
                                                                                                                               gTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGa                  <  1:155888/149‑1 (MQ=255)
                                                                                                                                 aCATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGAtt                >  7:14251/1‑149 (MQ=255)
                                                                                                                                  cATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGAtt                <  1:178207/148‑1 (MQ=255)
                                                                                                                                   aTGACGGCGGCGCCCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGa              >  2:107869/1‑149 (MQ=255)
                                                                                                                                    tGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGAt             <  8:19112/149‑1 (MQ=255)
                                                                                                                                     gACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATACCGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGAt             >  1:102399/1‑148 (MQ=255)
                                                                                                                                      aCGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATaa           <  6:209401/149‑1 (MQ=255)
                                                                                                                                        ggcggcGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATaaaa         <  2:102399/149‑1 (MQ=255)
                                                                                                                                         gcggcgACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAAt        <  2:226397/149‑1 (MQ=255)
                                                                                                                                               aCCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACt  >  4:352038/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TGCTGAACGAGCGTGACAGCACCATGATCATCATCTCGCACGACCGTCACTTCCTTAACATGGTCTGTACCCACATGGCGGATCTGGATTACGGCGAGCTGCGCGTTTATCCGGGTAACTACGATGAGTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACT  >  NC_000913/856552‑856843

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