Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,538,334 T→G 15.3% T174P (ACC→CCC)  narY ← nitrate reductase 2 (NRZ), beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,538,3340TG15.3% 143.8 / 2.0 73T174P (ACC→CCC) narYnitrate reductase 2 (NRZ), beta subunit
Reads supporting (aligned to +/- strand):  ref base T (19/38);  new base G (12/3);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 2.46e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCGTAGTAATCATCAATCTG  >  NC_000913/1538104‑1538568
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
atttggTGGGAAACAGTAATTTCTTTTTTCTGACTTGCCGTTTTTCAAGTTGAAGTGGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCCGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCaa                                                                                                                                                                                                                                      <  2:508643/234‑1 (MQ=255)
aaTTTGGGGGTAATAGAAGATGCTTTTTTCTGACTTGCCGCTTTTCCAGTGGAAGTAGATTTTGTTGGACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCaa                                                                                                                                                                                                                                      <  2:500288/239‑1 (MQ=255)
 tttCGTGGGTAACAGAATGTCCATTTTTCTGACTTGCCGGTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCCCAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCCCAGGCGCGGCAGGTACATCATGAAGGTGTTTTCaaa                                                                                                                                                                                                                                     <  2:891118/238‑1 (MQ=255)
     gTGGGTAACGAGATATGCATTTTTCTGACTTGCCTTTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGc                                                                                                                                                                                                                                 <  1:421018/239‑1 (MQ=255)
          tAACAGAATGGGCATTTTTTCGGGTTACCTCTTTTCCAGTTGAAGTGTATTTTTTTTGCAGGAAACCCGCTTATGCACAAACGCCAGCCGCGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCGCGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGt                                                                                                                                                                                                                             <  2:902779/238‑1 (MQ=255)
          tAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACCATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGGGGCAGGAGCATCCTGAAGGGGTTTTCAACCTGccccgc                                                                                                                                                                                                                            >  1:441407/1‑236 (MQ=255)
                 agtttctttttttgggcTTGCCGGTTTTCCAGTTAAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGGTCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACAt                                                                                                                                                                                                                         <  2:932282/219‑1 (MQ=255)
                 agaTGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACAt                                                                                                                                                                                                                         <  1:350382/235‑1 (MQ=255)
                  gaTGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGGGTTTTCAACATACCCGTACATCTCCt                                                                                                                                                                                                                    >  1:38930/1‑239 (MQ=255)
                        ttttttTGATTTTCCGTTTTTCCAGTTTAAGTAGATTTTTTTGTCAGGACAACCGCTTAGGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGGCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCa                                                                                                                                                                                                              <  2:854452/239‑1 (MQ=255)
                        tttttGTGACTGGCCGCTTTTCCAGTTGAAGTAGATTTTTTTTGACGGACCACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCa                                                                                                                                                                                                              <  1:442489/239‑1 (MQ=255)
                        tttttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCa                                                                                                                                                                                                              <  2:660605/239‑1 (MQ=255)
                        ttttcttgtcttgcgttttttCCATTTGAAGTAGATTTTTTTTTCCGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCGTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCa                                                                                                                                                                                                              <  2:472639/222‑1 (MQ=255)
                          tttCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGGACATACGTAAGGGGTTTTCCAACCACCCGGCAATCCCCTTTTGCGTg                                                                                                                                                                                                            >  2:811812/1‑239 (MQ=255)
                                aCTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGGGTTCGCACAGGGCGGGCAGGGACATCATGAAGGGGTTTTCAAACTGCCCGGCAAACCCCTTTTGAATgggccgg                                                                                                                                                                                                      >  2:314050/1‑234 (MQ=21)
                                          tttCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATAATGAAGGGGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCgg                                                                                                                                                                                             >  1:924550/1‑238 (MQ=255)
                                                              ttttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGCTTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACACCCTGACGGGGTTTTTCAACTGCCCGTACATCTCCTTTTTCTAGGCCCCGAAGTGCGGGCGCgggggccggttttgtaac                                                                                                                                                                      >  2:259565/1‑224 (MQ=255)
                                                                         cAACCGCTTATGCACAAACGCCAGCCACGGCATTTAACCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCGTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCcag                                                                                                                                                            <  1:254989/239‑1 (MQ=255)
                                                                                    gCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGAACATCATGAAGGTGTTTTCAAACTGCCCCTACATCTCCTTTTTCATGGCCTCGAAGTTTGGGGCGCGGGCCCTGTTTTCGAAATCACCGCCCACCAAGTCTTccc                                                                                                                                                 >  2:476612/1‑239 (MQ=255)
                                                                                      acaAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGGTTTCAAACTGCCCGTCCATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCCTGGGCACGTTTTTCGCACCCACCGCCCCGCAGTTTTTCCCAg                                                                                                                                               >  1:950144/1‑239 (MQ=255)
                                                                                          aCGCCCTCCACGGCCTTTATCCATATCAATCAGGAAAATGCGGGCTTCTGCGGGTTTGTGGATCGTGCCGCTTGGGCCGGGCGCCACGCAGCTGGAATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTg                                                                                                                                           <  2:738621/239‑1 (MQ=255)
                                                                                          aCGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGGGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCCCGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCCGGTTg                                                                                                                                           >  1:356018/1‑239 (MQ=255)
                                                                                            gccagcctCGGCATTTATCCTGATCAATCAGCACAAGGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGc                                                                                                                                         <  1:314050/239‑1 (MQ=255)
                                                                                              cagccaCGGCATTTATCCTGATCAATCAGCAAAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGccc                                                                                                                                       <  1:189276/239‑1 (MQ=255)
                                                                                              cagccaCGGCATTTAACCTGATCAATCAGCCCAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGccc                                                                                                                                       <  1:259565/239‑1 (MQ=255)
                                                                                                        tttgatccttgataagcggCACAAGTCCGTCTTCGTTGCGTTTTGAGATGGCGCCGGTTGGGCAGGTCGCCGCGCAGCTGGGATGGAGGCAGGGTTCGCACAGGCGAGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACt                                                                                                                             <  2:214693/221‑1 (MQ=255)
                                                                                                          ttATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGGACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGGACTCACCGCCCCGCAGTTCTTCCCAGTTTGGCCCCCAGATCACttt                                                                                                                           >  1:80880/1‑239 (MQ=255)
                                                                                                               cTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCCTTTGCATGTCCTCGAAGTTGCGGGCGCGGGCACGGTTTTCGAACTCACAGCCCGGCAGTTCTTCCCCGTTTGGCCCCCAGATCCCttttgcc                                                                                                                       >  1:994074/1‑234 (MQ=255)
                                                                                                                 gATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGCAGGGGTTTTCAAACTGCCCTGACATCTGCTTTTTCATGGCCTCGCAGTTGCGGTCGCCGGCCCGGTTTTCGGACTCCCCCCCCAGCAGTTCTTCCCCGTTTTGCCCCCAGAACaacttggccctc                                                                                                                    >  2:587817/1‑228 (MQ=255)
                                                                                                                 gATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGGGTTTTCAAACTGCCCGTACCTCCCCTGTTGCCGGGCCCCGAAGTGTCGTCGTCGGGCAGGTTTTTTGACCTAACCCCCCACCAGTTTTTCCCCTTTTGGccccacgacacctttgtcatc                                                                                                                     >  2:1078017/1‑221 (MQ=255)
                                                                                                                 gATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGGGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCCCGGGCCCGTTTTTCGAACTTCCCCCCCCACAGGTTTTTCCCGTTTGGCCCCCAGAACACTTTGGCCATc                                                                                                                    >  2:592412/1‑239 (MQ=255)
                                                                                                                  atcaatcaGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGGGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGGTTGGCCCCCAGATCACTTTGTcactcc                                                                                                                   >  1:430720/1‑234 (MQ=255)
                                                                                                                       tcaGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGGGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTAGCGGTTGCGGGCACGTTTTTCGAACTTACCACCCACCAGGTCTTCCCAGTTTTGCACCCAGATCACTTTTGCCCTCCGttt                                                                                                               >  1:687262/1‑238 (MQ=255)
                                                                                                                       tcaGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGGACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGGGGGCACGTTTTTCGAACTCACCGCCCGGCAGTTCTTCCCATTTTGGCCCCCGGGTAACTTTTTCCATCCGttt                                                                                                               >  1:1077896/1‑238 (MQ=255)
                                                                                                                             ccatcccGTCTTCTTCGCGTTTGTAGATAGCGCCGCTGGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCa                                                                                                        <  2:924550/234‑1 (MQ=255)
                                                                                                                                   cggtttcttcGCGTTTGTAGATAGTGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCa                                                                                                    <  1:637326/233‑1 (MQ=255)
                                                                                                                                     tcttctTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCGGCAGTTCTTCCCAGTTTGGCCCCCCGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGa                                                                                                >  1:756293/1‑239 (MQ=255)
                                                                                                                                     tcttctTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGa                                                                                                >  1:419189/1‑239 (MQ=255)
                                                                                                                                         ttacgcgTTTGTAGATTGCCCCGGTTTGGCGGGGCACCCGGCAGGTGGGATGGAGGCAGTGGTCGCACAGGGCCGGCGGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAATTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAAcgc                                                                                            <  2:160296/236‑1 (MQ=255)
                                                                                                                                          ttcGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGCTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATTGCCTCGCAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTTACAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAATGAActcg                                                                                           >  1:572760/1‑239 (MQ=255)
                                                                                                                                              cggttGTAGTTGGCGGCGGTTGGGAAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCg                                                                                       <  1:582333/236‑1 (MQ=255)
                                                                                                                                               gTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCgg                                                                                      >  1:633601/1‑238 (MQ=255)
                                                                                                                                                 ttgtgtaaggcgcCGGTGGGGCGGGGTGCACCGCAGCGGGGTTTGGGGCAGTGGTCTCACAGGCGTGGCAGGTACATCATGAAGGTGTTTTCAAGCTGCCCGTACATCTCCTTTTGCATGTCCTCGAAGTTGCGGTCGCGGGCCCGTTTTTCGAACTCACCGCCCAGCAGTTCTTCACAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTAAATCAGTGAACGCGGGCGAg                                                                                     <  2:572760/229‑1 (MQ=255)
                                                                                                                                                 ttGGAGGTAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAg                                                                                     <  1:592412/238‑1 (MQ=255)
                                                                                                                                                  tttgaattcccccggTTGGGCCGGGTGCCCCGCAGCTGGGATTGAGGCAGTTTTCGCACAGGGGGGGCAGGTACATCATGAAGGTTTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCa                                                                                   <  2:873921/228‑1 (MQ=255)
                                                                                                                                                  ggTAGAGAGCGCCGCTTGGGCCGGTCGCCACGTCGCTGGGGTTGAGTCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTTCCCGTACATCTCCTTGTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCa                                                                                   <  2:687262/238‑1 (MQ=255)
                                                                                                                                                    tgaatggcgcCGCTTGGGCCGGTCGCCACGCAGGTGGGGTGGAGGCAGTGGTCGCAGAGGCGGGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGt                                                                                 <  2:80880/233‑1 (MQ=255)
                                                                                                                                                    taaaatgggccggTTGGGCGGGTCGCCACGCAGCGGGGATGGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGTGGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGt                                                                                 <  2:860237/231‑1 (MQ=255)
                                                                                                                                                           cgccgcTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCCGTTCTTCCCCGTTTGGCCCCCAGATCACTTTGTCCCTCCGTTTGCCGTCAATCAAGGAAccgggggggcgcatggggaaa                                                                          >  2:276303/1‑219 (MQ=255)
                                                                                                                                                           cgccgcTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATAATGAAGGTGTTTTCAAACTGCCCGGAACTCTACTTTTGCCTGGCCTCGAAGTTGCGGTCCTGGGCCCGCTTTTCGAACTTACCGCCCAGCCGTTCTTCCCAAGTTGGGCCCCCGTTCCCCTTGTACCTCCGTTTCCCGGCAACCAGTGTACGCAGGGGacgcggtggaata                                                                          >  2:120382/1‑227 (MQ=255)
                                                                                                                                                              cggttGGGCAGGTGGCCCCGCAGGGGGGATTGAGGCAGGGTTCGCACAGGCGCGGCAGGTACATCATGAAGGGGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGt                                                                       <  2:756293/236‑1 (MQ=255)
                                                                                                                                                              cgcTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCCCGGGCACGGTTTTCGAACTTACCGCCACGCAGGTTTTCCCCGTTTGGCCCCCCGATCACCTTTTCCATCCGTTTTCCGTCAAACAGGGAACCGGGGCGGGCAGTCGGaaaagtt                                                                       >  2:698533/1‑234 (MQ=255)
                                                                                                                                                                        ggTCGCCCCGGAGCGGGGATTGGGGGAGTGTTCGCGCGGGCGCGGGAGGGACATCATGAAGGTGTTTTCAAACCGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCg                                                             <  2:586435/239‑1 (MQ=255)
                                                                                                                                                                        ggTCGCCACGCAGCTGGGGTTGAGGCAGTGTTCGCCCAGGGGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCg                                                             <  2:356018/239‑1 (MQ=255)
                                                                                                                                                                           cGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACCCGGGCGAGCATTAGGAATATGTTTGCCTTCCggg                                                          >  1:294113/1‑237 (MQ=255)
                                                                                                                                                                              cccgcgGCGGGGGTTGGGGCAGTTTTCGCACAGGCGGGGCAGGTACATCATGAAGGGGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTgcg                                                       <  2:950144/236‑1 (MQ=255)
                                                                                                                                                                              cACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTgcg                                                       <  1:697373/238‑1 (MQ=255)
                                                                                                                                                                                  ccggtGGGATTGGGGGAGTGTTCGCACGGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGCCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATg                                                  <  2:419189/235‑1 (MQ=255)
                                                                                                                                                                                                        cGCACAGGCGCGGCAGGTACATCATGAAGGAGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCCCGGGCACGGTTTTCGAACTCACCGCCCAGCAGTTCTTCCCCGTTTGGCCCCCAGATCACTTTTTCCATCCCTTTTCCGTCAATCCGTTAACCCGGGCGAGCAGTAGGGATATGTTTTCTTTCCGTGTGGCTATGACAATTTTCGTAGGCGGAgggg                            >  2:629651/1‑237 (MQ=255)
                                                                                                                                                                                                         gCACAGGCGCGGCAGGTACCACATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGa                           <  1:208816/239‑1 (MQ=255)
                                                                                                                                                                                                         ccccaGGGGCGGCAGGGTCATCATGAAGGTGTTTTCAAACTGCCCGGTCATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGa                           <  2:633601/236‑1 (MQ=255)
                                                                                                                                                                                                              ggCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGACACCTTTGTCCATCCGTTTTCCGTCCATCCGTGAACCCGGCGGCGCGACAGGAACATGTTTTCCCTTCGGGGGGCTCGGCAAAGGTTCTAAGTCTCAGGGGAAGGGt                      >  2:206099/1‑239 (MQ=255)
                                                                                                                                                                                                                cgcgGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAAATCACCGCCCCGCAGTTCTTCCCAGTTTGGCCCCCCGATAACTTTGTCCATCCGTTTGCCGTCAATCCGTGAACGCGGGCTGGCCGTGGGAATATGTTTCCCTTCCGGTGCGCTTTCCAAATGTTGCAATTCGAGGGGGAAAGGtac                    >  2:180185/1‑237 (MQ=255)
                                                                                                                                                                                                                  gggCAGGTAAATCATGAAGGTGTTTTCAAACTGCCCGTACTTCGCCTTTTGCATGGCCTAGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCg                   <  1:576968/237‑1 (MQ=255)
                                                                                                                                                                                                                  cgggtggTCCATGAGGAAGGTGTTTTCAAACTGCCCGTACATTCCCTTTGGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCg                   <  1:476612/233‑1 (MQ=255)
                                                                                                                                                                                                                  cgGCAGGTACATCATGAAGGTTTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCg                   <  1:698533/238‑1 (MQ=255)
                                                                                                                                                                                                                  cgGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGGGCAGTAGGAATATGTTTGCCTTCCGCTGCGCTATGCAAATGTTCGTAGTCGAAGGGGAACGGTTCgt                  >  1:813942/1‑239 (MQ=255)
                                                                                                                                                                                                                   ggCAGGTACACTAGGAAGGTTTTTCCAAATGGCCCGTACATCTCCTTTGGCATGGCCCCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCgta                 <  2:620800/239‑1 (MQ=255)
                                                                                                                                                                                                                    gCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGGTTGACTTCCGGGTCGCTATGCAAATGTTCGTAGTCGAAGGGGAAAGGTTCggag                >  1:620800/1‑239 (MQ=255)
                                                                                                                                                                                                                      aGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCgtagt               <  1:276303/238‑1 (MQ=255)
                                                                                                                                                                                                                        tttcatcatGTAGGTTTTTGTAAACTTCCCGTACATCTCCTTTTGCATGGCCCCTAAGTGGCGGTGGCGGGCACGTTTTTCGAACTCCCCGCCCTGGAGTGCTTCCCAGTTTGGCGCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCGTAGTAat            <  2:317628/236‑1 (MQ=255)
                                                                                                                                                                                                                         tACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCCCTTTGTCCATCCGTTTGCCGTCAATCAGTGAACCCGGGCGGCCAGTAGGAATATGTTTGCCTTCCGGGGCGCCATGCAAAGGTTCGTAGCCGAAGGGGAAAGGTTCggatgaaat           >  2:63984/1‑231 (MQ=255)
                                                                                                                                                                                                                                  aaGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATTTGTTTGCCTTCCGGTTCGCTATGAAAATGTTCGTAGTCGAAGGTTAAAGGGTCGAGGCAATCATAAATCtc  >  2:767297/1‑238 (MQ=255)
                                                                                                                                                                                                                                      tgtTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTgcg                                                       <  2:746056/182‑1 (MQ=255)
                                                                                                                                                                                                                                      tgtTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTgcg                                                       >  1:746056/1‑182 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
AATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCGTAGTAATCATCAATCTG  >  NC_000913/1538104‑1538568

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