Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,949,005 T→C Q322R (CAG→CGG)  yifB ← magnesium chelatase family protein and putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,949,0050TC99.3% 520.8 / ‑4.1 140Q322R (CAG→CGG) yifBmagnesium chelatase family protein and putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (72/67);  total (73/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

Only 100 of 140 total aligned reads displayed.

AGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATG  >  NC_000913/3948858‑3949152
                                                                                                                                                   |                                                                                                                                                   
agCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGcc                                                                                                                                                  >  1:701068/1‑151 (MQ=255)
  ccGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCg                                                                                                                                                <  2:574725/151‑1 (MQ=255)
   cGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCgg                                                                                                                                               >  2:986775/1‑151 (MQ=255)
    ggTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGa                                                                                                                                              <  3:546197/151‑1 (MQ=255)
    ggTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGa                                                                                                                                              <  2:917929/151‑1 (MQ=255)
        gagaTAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCa                                                                                                                                          <  4:58439/151‑1 (MQ=255)
          gaTAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAAt                                                                                                                                        <  1:1684952/151‑1 (MQ=255)
           aTAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATc                                                                                                                                       <  1:166314/151‑1 (MQ=255)
           aTAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATc                                                                                                                                       >  2:663034/1‑151 (MQ=255)
             aaCGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCgg                                                                                                                                     <  2:1407562/151‑1 (MQ=255)
               cGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGct                                                                                                                                   <  1:103096/151‑1 (MQ=255)
               cGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGct                                                                                                                                   <  1:983000/151‑1 (MQ=255)
               cGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGct                                                                                                                                   <  2:215433/151‑1 (MQ=255)
                gTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGctc                                                                                                                                  <  1:869756/151‑1 (MQ=255)
                  aaTGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCtctc                                                                                                                                <  1:147289/151‑1 (MQ=255)
                   aTGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCg                                                                                                                               >  1:1328438/1‑151 (MQ=255)
                         gTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAgg                                                                                                                          <  4:6306/150‑1 (MQ=255)
                         gTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGc                                                                                                                         <  4:254863/151‑1 (MQ=255)
                         gTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGc                                                                                                                         <  2:454138/151‑1 (MQ=255)
                           tCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCa                                                                                                                        <  2:1537673/150‑1 (MQ=255)
                           tCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCa                                                                                                                        <  1:374041/150‑1 (MQ=255)
                           tCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCAt                                                                                                                       >  1:909198/1‑151 (MQ=255)
                             tGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATcc                                                                                                                     >  3:289744/1‑151 (MQ=255)
                              ggCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATcc                                                                                                                     >  4:394564/1‑150 (MQ=255)
                                 gTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAgtg                                                                                                                 <  1:1330001/151‑1 (MQ=255)
                                 gTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAgt                                                                                                                  <  2:136033/150‑1 (MQ=255)
                                   gCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTa                                                                                                               <  2:929199/151‑1 (MQ=255)
                                   gCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTa                                                                                                               <  3:314589/151‑1 (MQ=255)
                                     aGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACg                                                                                                             >  4:455800/1‑151 (MQ=255)
                                     aGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACg                                                                                                             >  2:802026/1‑151 (MQ=255)
                                     aGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACg                                                                                                             >  2:1493083/1‑151 (MQ=255)
                                        ggTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCg                                                                                                          >  1:118427/1‑151 (MQ=255)
                                        ggTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCg                                                                                                          <  1:482188/151‑1 (MQ=255)
                                            aTGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTc                                                                                                       >  2:1534248/1‑150 (MQ=255)
                                              ggTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCaa                                                                                                     >  1:1636185/1‑150 (MQ=255)
                                               gTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAAt                                                                                                   <  2:1244914/151‑1 (MQ=255)
                                                   cccTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAg                                                                                               <  2:658875/151‑1 (MQ=255)
                                                       gATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTa                                                                                            <  3:125083/150‑1 (MQ=255)
                                                           aTGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCtca                                                                                       <  2:1289670/151‑1 (MQ=255)
                                                           aTGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCtc                                                                                        >  1:1186143/1‑150 (MQ=255)
                                                             gTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCtcatc                                                                                     >  2:1284991/1‑151 (MQ=255)
                                                             gTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCtcat                                                                                      <  2:1385888/150‑1 (MQ=255)
                                                               ccGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCtcatca                                                                                    <  4:430476/150‑1 (MQ=255)
                                                                  gTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGaa                                                                                <  2:909312/151‑1 (MQ=255)
                                                                  gTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGaa                                                                                <  2:549909/151‑1 (MQ=255)
                                                                  gTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGa                                                                                 >  2:629826/1‑150 (MQ=255)
                                                                    aGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGaaa                                                                               >  1:1108786/1‑150 (MQ=255)
                                                                    aGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGaaa                                                                               <  2:53554/150‑1 (MQ=255)
                                                                         tGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCAc                                                                         >  3:572459/1‑151 (MQ=255)
                                                                          gggATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCa                                                                          <  1:1629712/149‑1 (MQ=255)
                                                                             aTTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGcct                                                                     <  2:1636349/151‑2 (MQ=255)
                                                                                 aTCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGtt                                                                   <  1:399477/149‑1 (MQ=255)
                                                                                 aTCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTAt                                                                 >  1:814020/1‑151 (MQ=255)
                                                                                 aTCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTAt                                                                 >  1:1738993/1‑151 (MQ=255)
                                                                                  tCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATg                                                                >  1:114759/1‑151 (MQ=255)
                                                                                   cgccgcAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATgc                                                               >  2:429657/1‑151 (MQ=255)
                                                                                       gcAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGcc                                                            >  4:218469/1‑150 (MQ=255)
                                                                                        caacaaGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCa                                                           <  1:454020/150‑1 (MQ=255)
                                                                                        caacaaGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCa                                                           <  2:745193/150‑1 (MQ=255)
                                                                                         aacaaGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGc                                                         <  3:571632/151‑1 (MQ=255)
                                                                                            aaGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGa                                                       <  2:713083/150‑1 (MQ=255)
                                                                                                tGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAAttt                                                  >  1:81093/1‑151 (MQ=255)
                                                                                                  gAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTc                                                 >  1:728799/1‑150 (MQ=255)
                                                                                                  gAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCa                                                <  2:1043048/151‑1 (MQ=255)
                                                                                                   aaaCGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCAc                                               >  1:1130349/1‑151 (MQ=255)
                                                                                                     aCGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCg                                             >  1:449480/1‑151 (MQ=255)
                                                                                                       gggCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCggg                                           <  1:251222/151‑1 (MQ=255)
                                                                                                       gggCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCggg                                           <  2:563444/151‑1 (MQ=255)
                                                                                                          cTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGccc                                        >  4:521162/1‑151 (MQ=255)
                                                                                                           tGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCt                                       >  1:1024339/1‑151 (MQ=255)
                                                                                                            ggATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCt                                       <  1:1685462/150‑1 (MQ=255)
                                                                                                             gATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTgg                                     >  2:528524/1‑151 (MQ=255)
                                                                                                                   ttATTTTTGCTCGTGTGGGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTg                                <  2:639645/150‑1 (MQ=255)
                                                                                                                    tATTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTgc                               <  1:1388598/150‑1 (MQ=255)
                                                                                                                     aTTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTgcgc                             >  1:1202705/1‑151 (MQ=255)
                                                                                                                     aTTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTgcgc                             >  1:245102/1‑151 (MQ=255)
                                                                                                                     aTTTTTGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTgcgc                             >  2:271071/1‑151 (MQ=255)
                                                                                                                       ttttGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGcc                            >  4:385201/1‑150 (MQ=255)
                                                                                                                         ttGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCAc                          <  2:285722/150‑1 (MQ=255)
                                                                                                                          tGCTCGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCg                        >  1:735075/1‑151 (MQ=255)
                                                                                                                              cGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCt                     >  1:1467343/1‑150 (MQ=255)
                                                                                                                              cGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCt                     >  2:1661046/1‑150 (MQ=255)
                                                                                                                              cGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCt                     <  2:118448/150‑1 (MQ=255)
                                                                                                                              cGTGTGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTa                    <  1:1384086/151‑1 (MQ=255)
                                                                                                                               gtgtGCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCt                     <  4:144523/149‑1 (MQ=255)
                                                                                                                                tgtgCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTAcc                  >  3:473551/1‑151 (MQ=255)
                                                                                                                                tgtgCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTAc                   <  2:598052/150‑1 (MQ=255)
                                                                                                                                 gtgCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTg                                >  3:414416/1‑136 (MQ=255)
                                                                                                                                   gCGTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATc               >  1:238938/1‑151 (MQ=255)
                                                                                                                                     gTGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGc             <  3:521178/151‑1 (MQ=255)
                                                                                                                                      tGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGc             >  2:1025403/1‑150 (MQ=255)
                                                                                                                                      tGAAAGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGc             >  3:40012/1‑150 (MQ=255)
                                                                                                                                        aaaGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGt          >  2:1154538/1‑151 (MQ=255)
                                                                                                                                        aaaGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGt          <  3:197242/151‑1 (MQ=255)
                                                                                                                                        aaaGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGt          >  2:1542098/1‑151 (MQ=255)
                                                                                                                                         aaGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGtt         <  2:345102/151‑1 (MQ=255)
                                                                                                                                          aGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGtt         >  2:1383098/1‑150 (MQ=255)
                                                                                                                                          aGATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGtt         >  2:454928/1‑150 (MQ=255)
                                                                                                                                           gATGGATCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTa        >  1:426083/1‑150 (MQ=255)
                                                                                                                                                aTCCGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATg  >  1:1397742/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATG  >  NC_000913/3948858‑3949152

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: