Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 18,905 A→C 5.1% E64A (GAG→GCG)  nhaR → transcriptional activator of nhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091318,9050AC 5.1% 232.1 / 3.0 79E64A (GAG→GCG) nhaRtranscriptional activator of nhaA
Reads supporting (aligned to +/- strand):  ref base A (37/38);  new base C (2/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.74e-01

TGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGC  >  NC_000913/18757‑19050
                                                                                                                                                    |                                                                                                                                                 
tGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAg                                                                                                                                                  <  2:329932/150‑1 (MQ=255)
tGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAg                                                                                                                                                  <  2:565757/150‑1 (MQ=255)
     tGTCTATACAGAAGGTTCCGTGGGGGGCGCAGCGGAGGCGCTTTATTTCACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAACCGCAAGGGACGTGGTCTCGAACCCAGCGAGCTggg                                                                                                                                             <  2:149852/150‑1 (MQ=255)
     tGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTggg                                                                                                                                             <  2:1088400/150‑1 (MQ=255)
       tCTATAAAGCAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGgag                                                                                                                                           <  2:828468/150‑1 (MQ=255)
       tCTATAAAGAATGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAATGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGgag                                                                                                                                           <  2:774408/150‑1 (MQ=255)
          ataAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAAc                                                                                                                                        >  2:1253125/1‑150 (MQ=255)
           taAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCTCGCTgggttcgag                                                                                                                                       <  1:883724/150‑7 (MQ=255)
            aaaGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCTCGCTgggttcgaga                                                                                                                                      <  2:1182827/150‑8 (MQ=255)
            aaaGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTg                                                                                                                                      >  1:29995/1‑150 (MQ=255)
              agcagGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAACCCATTACCGCACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGt                                                                                                                                    <  2:634747/147‑1 (MQ=255)
              agatgtttCCGTGGTTGCAGCTGAGTAGGCGCTTTATTTTACTCATAAAACCATTACCTGACAGATTCGAGCGCTGTTAGATTGACTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGt                                                                                                                                    <  2:1017923/144‑1 (MQ=255)
                    ttCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCg                                                                                                                              <  2:1270423/150‑1 (MQ=255)
                       cGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTa                                                                                                                           >  2:780640/1‑150 (MQ=255)
                        gTGGTTGGCGCAGCGGAGGCGCTTTATTTATCTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAATCGCAAGGGTCGTGGTCTCGAACCCAGCGAGCTGGGTGATCTGGTGTTTCGGTAt                                                                                                                          >  1:1224463/1‑150 (MQ=255)
                        gTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTAt                                                                                                                          >  2:234012/1‑150 (MQ=255)
                        gTGGTTGGCGCAGCGGAGGCGCGTTATTTAACTCCACAAACCATTACCGGACAGCTTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAAGTATTTTAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGTGATCTGGTCTTTGGCTAt                                                                                                                          >  1:1219593/1‑150 (MQ=255)
                          ggTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGc                                                                                                                        <  1:1228934/150‑1 (MQ=255)
                              ggCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGAt                                                                                                                    >  1:83614/1‑150 (MQ=255)
                               gcgcAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATa                                                                                                                   >  2:464880/1‑150 (MQ=255)
                                         ggCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCAc                                                                                                         <  2:685338/150‑1 (MQ=255)
                                         ggCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCAc                                                                                                         <  1:870653/150‑1 (MQ=255)
                                            gcTTTATTTAACTCCACACACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAATTTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGACTATAAAATGTTCACAtt                                                                                                      >  2:202108/1‑150 (MQ=255)
                                            gcTTTATTTAACTCCACAAACCATTACCGGACAGATTCGCGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCtt                                                                                                      >  2:121203/1‑150 (MQ=255)
                                            gcTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCtt                                                                                                      >  1:208241/1‑150 (MQ=255)
                                             cTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTa                                                                                                     <  1:1227305/150‑1 (MQ=255)
                                               ttatttaACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAg                                                                                                   >  1:1080213/1‑150 (MQ=255)
                                               ttatataACTCCACAATCCAATACCGGACTTATGCGAGCGCTGGAAGAGCGCCTGCATGGAAAATTATTTAAACGCTAGGGACGTGGTATCGAACCCAGCGAGGTGGGAGAACTGGTCTATCGCTATGCCGATAAAATTTTCACCTTAAg                                                                                                   <  1:1242075/150‑1 (MQ=255)
                                                 atttaACTACACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAATGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGcc                                                                                                 >  2:656659/1‑150 (MQ=255)
                                                   ttaGCTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCACGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACACAGCGAGCTGGGAGATCTGGTCTATCGGTATGCCGATAAAATGTTCACCTTAAGCCAg                                                                                               >  1:1008506/1‑150 (MQ=255)
                                                   ttaACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAATTTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAg                                                                                               >  1:1013697/1‑150 (MQ=255)
                                                   ttaACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAg                                                                                               <  2:753190/150‑1 (MQ=255)
                                                    taACTCCACAAACCATTACCGGACAGATTCGAGCGCGGGCAGAGCGCCTGCAAGGCAAATTATTTATACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTAACGCTATGCCGATAAAATGTTCACCTTAAGCCAgg                                                                                              <  1:149009/150‑1 (MQ=255)
                                                         ccACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAATTg                                                                                         >  2:954775/1‑150 (MQ=255)
                                                         ccACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATg                                                                                         >  2:101606/1‑150 (MQ=255)
                                                             aaaCCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTgg                                                                                     >  1:674089/1‑150 (MQ=255)
                                                             aaaCCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTgg                                                                                     <  2:873381/150‑1 (MQ=255)
                                                              aaCCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGa                                                                                    <  2:29995/150‑1 (MQ=255)
                                                                      ccGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGa                                                                            <  2:1159568/150‑1 (MQ=255)
                                                                        ggACAGATTCTAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGCACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAAc                                                                          <  1:1140786/150‑1 (MQ=255)
                                                                             gATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTAAAACGCAAGGGACGTGGTCTCGAACCCTGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCg                                                                     <  1:1165854/150‑1 (MQ=255)
                                                                              aTTCGAGCGCTGTAAGAGCGCCTGCAATGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGc                                                                    <  1:1108363/150‑1 (MQ=255)
                                                                                 cGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCaaa                                                                 >  1:92497/1‑150 (MQ=255)
                                                                                       cTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATACCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATcc                                                           <  1:11361/150‑1 (MQ=255)
                                                                                       cTGGAAGAGCGCCTGCAAGGAAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAATCCAGGTAATGCTGGATTTTGTGAACTATCGCAAAGAATcc                                                           <  2:563090/150‑1 (MQ=255)
                                                                                        tGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCa                                                          >  2:910261/1‑150 (MQ=255)
                                                                                         ggAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATTGCAAAGAATCCaa                                                         >  2:1134649/1‑150 (MQ=255)
                                                                                          gAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCGGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAAt                                                        <  2:1055668/150‑1 (MQ=255)
                                                                                              agCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTAt                                                    <  2:263061/150‑1 (MQ=255)
                                                                                              agCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTAt                                                    <  2:724580/150‑1 (MQ=255)
                                                                                                cgcCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCTATAATATGTTCACCTTAATCCAGGAAATGCTGGAGATTGTGAACTATGGCAATGAATCCAattgattg                                                  >  2:642160/1‑150 (MQ=255)
                                                                                                  ccTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGtt                                                <  2:1113261/150‑1 (MQ=255)
                                                                                                      cAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGAc                                            <  1:234012/150‑1 (MQ=255)
                                                                                                        aGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGt                                          <  2:41096/150‑1 (MQ=255)
                                                                                                          gCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTg                                        >  2:139806/1‑150 (MQ=255)
                                                                                                                tattTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTgg                                  >  1:1025500/1‑150 (MQ=255)
                                                                                                                 attTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGc                                 <  1:671316/150‑1 (MQ=255)
                                                                                                                    tAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGa                              >  1:1023506/1‑150 (MQ=255)
                                                                                                                      aaCGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATg                            <  1:1063829/150‑1 (MQ=255)
                                                                                                                      aaCGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATg                            <  1:940158/150‑1 (MQ=255)
                                                                                                                        cGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGCATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCa                          >  1:390643/1‑150 (MQ=255)
                                                                                                                            aGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACttt                      <  1:1253125/150‑1 (MQ=255)
                                                                                                                             gggACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTc                     <  1:1240806/150‑1 (MQ=255)
                                                                                                                             gggACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTc                     >  1:298095/1‑150 (MQ=255)
                                                                                                                             gggACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTAACGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTc                     <  2:83614/150‑1 (MQ=255)
                                                                                                                              gcaCGTGGTCTCGAACCCAGCAAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCATAAGCCAGGAAAAGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGTCTGATGCACTTTcc                    >  1:190376/3‑150 (MQ=255)
                                                                                                                                   tGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGAAAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACg               <  2:135371/150‑1 (MQ=255)
                                                                                                                                    ggTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGc              >  1:864205/1‑150 (MQ=255)
                                                                                                                                    ggTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGc              >  1:699473/1‑150 (MQ=255)
                                                                                                                                      tctcGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCt            >  2:1145714/1‑150 (MQ=255)
                                                                                                                                       ctcGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGAACTTTCCAAACGCCTg           >  2:338231/1‑150 (MQ=255)
                                                                                                                                        tcGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTgg          <  1:973536/150‑1 (MQ=255)
                                                                                                                                          tctcccAGCTCGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTc        >  1:1060055/4‑150 (MQ=255)
                                                                                                                                          gAACCCAGCGAGCTGGGAGAACTGGTTTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTc        >  2:1169392/1‑150 (MQ=255)
                                                                                                                                           aaCCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCa       >  2:3804/1‑150 (MQ=255)
                                                                                                                                              ccAGCTCGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTa    >  1:175356/1‑150 (MQ=255)
                                                                                                                                              ccAGCGAGCTGGGAGAACTGGTCTATCGCTTTGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCATTTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGtt    >  1:815856/1‑149 (MQ=255)
                                                                                                                                               cAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAg   >  2:449783/1‑150 (MQ=255)
                                                                                                                                                aGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGc  >  1:328255/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGC  >  NC_000913/18757‑19050

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: