Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,586,450 C→A 6.1% L240L (CTC→CTA narQ → sensory histidine kinase in two‑component regulatory system with NarP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,586,4500CA 6.1% 308.8 / 2.8 98L240L (CTC→CTAnarQsensory histidine kinase in two‑component regulatory system with NarP
Reads supporting (aligned to +/- strand):  ref base C (50/41);  new base A (2/5);  total (52/46)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

AACACGGGCAGTTCGACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGTTAAATGTCGGTGAA  >  NC_000913/2586302‑2586597
                                                                                                                                                    |                                                                                                                                                   
aaCACGGGCAGTTCGACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTAGGTCTGCTTTCAAAAACCTTTAACCAGCTGTCGCGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTCGAAGAAAAGACCAGCGAtctc                                                                                                                                                     >  6:240416/1‑149 (MQ=255)
      ggCAGTTCGACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAg                                                                                                                                               <  2:364440/149‑1 (MQ=255)
       gCAGTTCGACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGTATTTTTAAGTATAAGTATAGACCCGAGATCTACACGAgg                                                                                                                                              <  8:231104/149‑1 (MQ=255)
         aGTTCGACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGCTGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGcc                                                                                                                                            <  5:171823/149‑1 (MQ=255)
         aGTTCGACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGc                                                                                                                                             >  5:222064/1‑148 (MQ=255)
           ttCGACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAgg                                                                                                                                              >  8:203751/1‑145 (MQ=255)
           tgcGACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAgg                                                                                                                                              <  7:203751/143‑1 (MQ=255)
              gACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAg                                                                                                                                               >  4:195239/1‑140 (MQ=255)
              gACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAg                                                                                                                                               <  3:195239/140‑1 (MQ=255)
                 tCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGcgtcg                                                                                                                                    >  1:86329/1‑149 (MQ=255)
                 tCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGcgtcg                                                                                                                                    >  3:82764/1‑149 (MQ=255)
                   gccgccgcTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGcgtcgt                                                                                                                                   <  1:4014/148‑1 (MQ=255)
                      gccgcTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGACCAAGCGTCGTCTg                                                                                                                                >  6:1166/1‑148 (MQ=255)
                        cggtGGCTACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGCGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAg                                                                                                                             <  7:48871/146‑1 (MQ=255)
                            ggATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCATTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGgtgt                                                                                                                         <  8:343000/149‑1 (MQ=255)
                            ggATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCCCGAGGCCAAGCGTCGTCTGGAGgtgt                                                                                                                         >  2:250272/1‑149 (MQ=255)
                                 ccAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGa                                                                                                                              >  1:107080/1‑139 (MQ=255)
                                 ccAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGa                                                                                                                              <  2:107080/139‑1 (MQ=255)
                                  cAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTAt                                                                                                                    >  8:75323/1‑148 (MQ=255)
                                  cAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTAt                                                                                                                    >  4:170416/1‑148 (MQ=255)
                                   aaCCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGAGATCTACACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCa                                                                                                                  <  4:206929/149‑1 (MQ=255)
                                     ccTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAgt                                                                                                                <  7:75323/149‑1 (MQ=255)
                                          cGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCg                                                                                                           >  2:40478/1‑149 (MQ=255)
                                          cGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCg                                                                                                           >  4:39944/1‑149 (MQ=255)
                                              tgcgcgtGGTCTGCTTGCAAAAACCTTTACCCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGTATTTTTAAGTATAAGAAAAGACCCGAGATCTACACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAgg                                                                                                       <  2:86329/143‑1 (MQ=255)
                                                aGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGgcg                                                                                                     >  6:153341/1‑149 (MQ=255)
                                                 gCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGCAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGgcgc                                                                                                    >  2:163443/1‑149 (MQ=255)
                                                   ttGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCt                                                                                                   >  7:371573/1‑148 (MQ=255)
                                                   ttGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCt                                                                                                   >  7:235553/1‑148 (MQ=255)
                                                     ggTCGGCTTGCAACAACCTTTAACCAGAGGTCGCGCGAGCTGCATACATTGTACCGGTCGCTGGAAGCGTCAGGAGAAGAACAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGaa                                                                                                <  5:240416/149‑1 (MQ=255)
                                                       tCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAaca                                                                                              <  3:64724/149‑1 (MQ=255)
                                                        cTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAacac                                                                                             >  8:309819/1‑149 (MQ=255)
                                                         tGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAacac                                                                                             >  3:96803/1‑148 (MQ=255)
                                                         tGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGACGCCAAGCGTCGTCTGGAGGTCTTGTATAAGTGTTCTCAGGCGCTGAacac                                                                                             >  2:15189/1‑148 (MQ=255)
                                                          gCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAacac                                                                                             >  4:212278/1‑147 (MQ=255)
                                                           ctagCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCCCGATCTCCACGAGGCCAAGCGTCGTCTGGATGTGTTGTATCAGTCTTCGCAGGCGCTGAACACt                                                                                            <  5:169252/144‑1 (MQ=255)
                                                           cTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGAACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGAGCTCCACGAGGCCAAGCGTCGTCTGGAGGTGATGTATCAGTGTTCGCAGGCGCTGAACACTAg                                                                                          <  5:187786/149‑1 (MQ=255)
                                                              gCAAACACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCTTTTTTTTTTTAATTATAAGTATAGAACCGAGATCTACACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGc                                                                                         <  2:65847/147‑1 (MQ=255)
                                                              gCAAAAACCTTTAACCCGATGTCGAGCGAGCTGCATAAATTGTACCGGTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGc                                                                                         <  1:117361/147‑1 (MQ=255)
                                                               cAAAAACCTGTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGGCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGcgtc                                                                                                                                     >  7:135296/1‑102 (MQ=255)
                                                                aaaaaCCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGGCGTCTGGAGGTGTTGTATCAGTGTTCGCATGCTCTGCACACTatc                                                                                         >  2:296593/1‑143 (MQ=255)
                                                                  aaaCCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGAtt                                                                                   >  3:67003/1‑149 (MQ=255)
                                                                  aaaCCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGAtt                                                                                   >  7:326680/1‑149 (MQ=255)
                                                                  aaaCCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGAtt                                                                                   >  6:88412/1‑149 (MQ=255)
                                                                   aaCCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGAtt                                                                                   >  6:3982/1‑148 (MQ=255)
                                                                    aCCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGAGTAAGTAGAAGAAAAGAACCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGa                                                                                 <  8:171272/149‑1 (MQ=255)
                                                                     ccTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGAt                                                                                <  3:170416/149‑1 (MQ=255)
                                                                     ccTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGAt                                                                                <  3:302765/149‑1 (MQ=255)
                                                                            ccAGATGCCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCAtt                                                                         >  6:40442/1‑149 (MQ=255)
                                                                              aGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGt                                                                       >  2:39341/1‑149 (MQ=255)
                                                                                      agcgagcTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCa                                                               >  4:64375/1‑149 (MQ=255)
                                                                                       gcgagcTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCa                                                               >  4:74782/1‑148 (MQ=255)
                                                                                            cTGCATAAATTGTACCGTTCGCTGGAAGCGTAAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATAtttt                                                         <  8:326680/149‑1 (MQ=255)
                                                                                            cTGCATAAATTGTACCGTTCGCTGGAAGAGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATAtttt                                                         <  3:39944/149‑1 (MQ=255)
                                                                                             tGCATAAATTGTACCGTTCGCTTTAATTGTAATTATAAGTATAGAACCGAGATCTACACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGAt                                                                                <  2:32600/125‑1 (MQ=255)
                                                                                             tGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGAt                                                                                >  1:32600/1‑125 (MQ=255)
                                                                                              gCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCAtt                                                                         >  4:283895/1‑131 (MQ=255)
                                                                                              gCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCAtt                                                                         <  3:283895/131‑1 (MQ=255)
                                                                                              gCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTg                                                        >  4:133991/1‑148 (MQ=255)
                                                                                               cATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCa                                                      <  1:39341/149‑1 (MQ=255)
                                                                                                aTAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTg                                                        >  5:79085/1‑146 (MQ=255)
                                                                                                 tAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGTAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGa                                                    <  8:235553/149‑1 (MQ=255)
                                                                                                 tAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGa                                                    <  5:88412/149‑1 (MQ=255)
                                                                                                     ttGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCa                                                      <  5:71249/143‑1 (MQ=255)
                                                                                                     ttGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCa                                                      >  6:71249/1‑143 (MQ=255)
                                                                                                     ttGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGt                                                <  7:70850/149‑1 (MQ=255)
                                                                                                     ttGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGt                                                >  7:327052/1‑149 (MQ=255)
                                                                                                        tACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTcg                                             <  8:371573/149‑1 (MQ=255)
                                                                                                        tACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTcg                                             >  4:110945/1‑149 (MQ=255)
                                                                                                          ccGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCa                                                                                                                  <  5:312138/78‑1 (MQ=255)
                                                                                                          ccGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTACACGAGTCAAATCGTCTTCTGGAGGTGTTTTATCa                                                                                                                  >  6:312138/1‑78 (MQ=255)
                                                                                                           cGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGa                                          >  2:239829/1‑149 (MQ=255)
                                                                                                                 gtGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATg                                     >  1:115096/2‑148 (MQ=255)
                                                                                                                  tGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAg                                                                                                                                         <  6:182037/47‑1 (MQ=255)
                                                                                                                  tGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAg                                                                                                                                         >  5:182037/1‑47 (MQ=255)
                                                                                                                    gAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAg                                  <  1:40478/148‑1 (MQ=255)
                                                                                                                    gAAGAGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGc                                 <  4:230648/149‑1 (MQ=255)
                                                                                                                         gTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTc                                              <  1:148180/131‑1 (MQ=255)
                                                                                                                         gTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTc                                              >  2:148180/1‑131 (MQ=255)
                                                                                                                         gTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCt                             >  5:199806/1‑148 (MQ=255)
                                                                                                                          tCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGa                           <  3:74782/149‑1 (MQ=255)
                                                                                                                           cAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTg                                                 <  7:88420/126‑1 (MQ=255)
                                                                                                                            aGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGaa                          >  3:298130/1‑148 (MQ=255)
                                                                                                                             gTAGAAGAAAAGCCCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAata                        >  5:328483/1‑149 (MQ=255)
                                                                                                                              tAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAatat                       >  5:45941/1‑149 (MQ=255)
                                                                                                                              tAGAAGAAAAGACCAGCGATCTACACGAGGCCAAGCGTCGTCTTGAGGTGTTTTATCAGAGTTCGCAGGAGCTCAACACTACCAAGAATGATGTGAATTGTTTCCGCAATATTTTGCAGATTGTTAGCGACAATGAACAGGCTGAatat                       >  7:142537/1‑149 (MQ=255)
                                                                                                                                gaagaaAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCt                     <  4:96803/149‑1 (MQ=255)
                                                                                                                                gaagaaAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGAATGATGTGCATTGTTTCCGCCATATTTTGAAGATTGTTcgc                                            >  7:24932/1‑126 (MQ=255)
                                                                                                                                   gaaAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTg                    >  5:133378/1‑147 (MQ=255)
                                                                                                                                   gaaAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGa                  <  3:212278/149‑1 (MQ=255)
                                                                                                                                    aaaaGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAg                 <  6:222064/149‑1 (MQ=255)
                                                                                                                                        gcccggagATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTTTCAGTTTTCGCAGGCGCTGAACACTAGCCAGATTGATGTTGAGTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGTTaa             <  2:100402/142‑1 (MQ=255)
                                                                                                                                         aCCCGCGAGCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTAGCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGAGGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGTTaaa            <  1:250272/149‑1 (MQ=255)
                                                                                                                                          cccGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGTTAAAt           <  6:328483/149‑1 (MQ=255)
                                                                                                                                           ccGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACTCTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTg                    <  7:242336/139‑1 (MQ=255)
                                                                                                                                           ccGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTg                    >  8:242336/1‑139 (MQ=255)
                                                                                                                                           ccGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCAGTGTTTCCGCGATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGTTAAATg          >  1:64293/1‑149 (MQ=255)
                                                                                                                                              cgATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGTTAAATGTCg       <  5:40442/149‑1 (MQ=255)
                                                                                                                                                   tcCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGAGACAATGAAGAGGCTGAATATCTGGAGTTAAATGTCGGTGa   >  4:101483/1‑148 (MQ=255)
                                                                                                                                                    ccACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGTTAAATGTCGGTGaa  >  3:244755/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AACACGGGCAGTTCGACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGTTAAATGTCGGTGAA  >  NC_000913/2586302‑2586597

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: