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

AATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCG‑GTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCGTAGTAATCATCAATCTG  >  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)
                                                              ttttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGCTTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACACCCTGACGGGGTTTTTCAACTGCCCGTACATCTCCTTTTTCTAGGCCCCGAAG‑TGCGGGCG‑Cgggggccggttttgtaac                                                                                                                                                                      >  2:259565/1‑224 (MQ=255)
                                                                         cAACCGCTTATGCACAAACGCCAGCCACGGCATTTAACCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCGTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCcag                                                                                                                                                            <  1:254989/239‑1 (MQ=255)
                                                                                    gCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGAACATCATGAAGGTGTTTTCAAACTGCCCCTACATCTCCTTTTTCATGGCCTCGAAGTTTGGGGCG‑CGGGCCCTGTTTTCGAAATCACCGCCCACCAAGTCTTccc                                                                                                                                                 >  2:476612/1‑239 (MQ=255)
                                                                                      acaAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGGTTTCAAACTGCCCGTCCATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCC‑TGGGCACGTTTTTCGCACCCACCGCCCCGCAGTTTTTCCCAg                                                                                                                                               >  1:950144/1‑239 (MQ=255)
                                                                                          aCGCCCTCCACGGCCTTTATCCATATCAATCAGGAAAATGCGGGCTTCTGCGGGTTTGTGGATCGTGCCGCTTGGGCCGGGCGCCACGCAGCTGGAATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTg                                                                                                                                           <  2:738621/239‑1 (MQ=255)
                                                                                          aCGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGGGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCCCGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCCGGTTg                                                                                                                                           >  1:356018/1‑239 (MQ=255)
                                                                                            gccagcctCGGCATTTATCCTGATCAATCAGCACAAGGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGc                                                                                                                                         <  1:314050/239‑1 (MQ=255)
                                                                                              cagccaCGGCATTTATCCTGATCAATCAGCAAAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGccc                                                                                                                                       <  1:189276/239‑1 (MQ=255)
                                                                                              cagccaCGGCATTTAACCTGATCAATCAGCCCAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGccc                                                                                                                                       <  1:259565/239‑1 (MQ=255)
                                                                                                        tttgatccttgataagcggCACAAGTCCGTCTTCGTTGCGTTTTGAGATGGCGCCGGTTGGGCAGGTCGCCGCGCAGCTGGGATGGAGGCAGGGTTCGCACAGGCGAGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACt                                                                                                                             <  2:214693/221‑1 (MQ=255)
                                                                                                          ttATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGGACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGGACTCACCGCCCCGCAGTTCTTCCCAGTTTGGCCCCCAGATCACttt                                                                                                                           >  1:80880/1‑239 (MQ=255)
                                                                                                               cTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCCTTTGCATGTCCTCGAAGTTGCGGGCG‑CGGGCACGGTTTTCGAACTCACAGCCCGGCAGTTCTTCCCCGTTTGGCCCCCAGATCCCttttgcc                                                                                                                       >  1:994074/1‑234 (MQ=255)
                                                                                                                 gATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGCAGGGGTTTTCAAACTGCCCTGACATCTGCTTTTTCATGGCCTCGCAGTTGCGGTCG‑CCGGCCCGGTTTTCGGACTCCCCCCCCAGCAGTTCTTCCCCGTTTTGCCCCCAGAACaacttggccctc                                                                                                                    >  2:587817/1‑228 (MQ=255)
                                                                                                                 gATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGGGTTTTCAAACTGCCCGTACCTCCCCTGTTGCCGGGCCCCGAAG‑TGTCGTCGTCGGGCAGGTTTTTTGACCTAACCCCCCACCAGTTTTTCCCCTTTTGGccccacgacacctttgtcatc                                                                                                                     >  2:1078017/1‑221 (MQ=255)
                                                                                                                 gATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGGGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCC‑CGGGCCCGTTTTTCGAACTTCCCCCCCCACAGGTTTTTCCCGTTTGGCCCCCAGAACACTTTGGCCATc                                                                                                                    >  2:592412/1‑239 (MQ=255)
                                                                                                                  atcaatcaGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGGGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGGTTGGCCCCCAGATCACTTTGTcactcc                                                                                                                   >  1:430720/1‑234 (MQ=255)
                                                                                                                       tcaGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGGGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTAGCGGTTG‑CGGGCACGTTTTTCGAACTTACCACCCACCAGGTCTTCCCAGTTTTGCACCCAGATCACTTTTGCCCTCCGttt                                                                                                               >  1:687262/1‑238 (MQ=255)
                                                                                                                       tcaGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGGACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑GGGGCACGTTTTTCGAACTCACCGCCCGGCAGTTCTTCCCATTTTGGCCCCCGGGTAACTTTTTCCATCCGttt                                                                                                               >  1:1077896/1‑238 (MQ=255)
                                                                                                                             ccatcccGTCTTCTTCGCGTTTGTAGATAGCGCCGCTGGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCa                                                                                                        <  2:924550/234‑1 (MQ=255)
                                                                                                                                   cggtttcttcGCGTTTGTAGATAGTGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCa                                                                                                    <  1:637326/233‑1 (MQ=255)
                                                                                                                                     tcttctTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCGGCAGTTCTTCCCAGTTTGGCCCCCCGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGa                                                                                                >  1:756293/1‑239 (MQ=255)
                                                                                                                                     tcttctTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGa                                                                                                >  1:419189/1‑239 (MQ=255)
                                                                                                                                         ttacgcgTTTGTAGATTGCCCCGGTTTGGCGGGGCACCCGGCAGGTGGGATGGAGGCAGTGGTCGCACAGGGCCGGCGGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAATTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAAcgc                                                                                            <  2:160296/236‑1 (MQ=255)
                                                                                                                                          ttcGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGCTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATTGCCTCGCAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTTACAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAATGAActcg                                                                                           >  1:572760/1‑239 (MQ=255)
                                                                                                                                              cggttGTAGTTGGCGGCGGTTGGGAAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCg                                                                                       <  1:582333/236‑1 (MQ=255)
                                                                                                                                               gTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCgg                                                                                      >  1:633601/1‑238 (MQ=255)
                                                                                                                                                 ttgtgtaaggcgcCGGTGGGGCGGGGTGCACCGCAGCGGGGTTTGGGGCAGTGGTCTCACAGGCGTGGCAGGTACATCATGAAGGTGTTTTCAAGCTGCCCGTACATCTCCTTTTGCATGTCCTCGAAGTTGCGGTCG‑CGGGCCCGTTTTTCGAACTCACCGCCCAGCAGTTCTTCACAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTAAATCAGTGAACGCGGGCGAg                                                                                     <  2:572760/229‑1 (MQ=255)
                                                                                                                                                 ttGGAGGTAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAg                                                                                     <  1:592412/238‑1 (MQ=255)
                                                                                                                                                  tttgaattcccccggTTGGGCCGGGTGCCCCGCAGCTGGGATTGAGGCAGTTTTCGCACAGGGGGGGCAGGTACATCATGAAGGTTTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCa                                                                                   <  2:873921/228‑1 (MQ=255)
                                                                                                                                                  ggTAGAGAGCGCCGCTTGGGCCGGTCGCCACGTCGCTGGGGTTGAGTCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTTCCCGTACATCTCCTTGTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCa                                                                                   <  2:687262/238‑1 (MQ=255)
                                                                                                                                                    tgaatggcgcCGCTTGGGCCGGTCGCCACGCAGGTGGGGTGGAGGCAGTGGTCGCAGAGGCGGGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGt                                                                                 <  2:80880/233‑1 (MQ=255)
                                                                                                                                                    taaaatgggccggTTGGGCGGGTCGCCACGCAGCGGGGATGGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGTGGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGt                                                                                 <  2:860237/231‑1 (MQ=255)
                                                                                                                                                           cgccgcTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCCGTTCTTCCCCGTTTGGCCCCCAGATCACTTTGTCCCTCCGTTTGCCGTCAATCAAGGAAccgggggggcgcatggggaaa                                                                          >  2:276303/1‑219 (MQ=255)
                                                                                                                                                           cgccgcTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATAATGAAGGTGTTTTCAAACTGCCCGGAACTCTACTTTTGCCTGGCCTCGAAGTTGCGGTCC‑TGGGCCCGCTTTTCGAACTTACCGCCCAGCCGTTCTTCCCAAGTTGGGCCCCCGTTCCCCTTGTACCTCCGTTTCCCGGCAACCAGTGTACGCAGGGGacgcggtggaata                                                                          >  2:120382/1‑227 (MQ=255)
                                                                                                                                                              cggttGGGCAGGTGGCCCCGCAGGGGGGATTGAGGCAGGGTTCGCACAGGCGCGGCAGGTACATCATGAAGGGGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGt                                                                       <  2:756293/236‑1 (MQ=255)
                                                                                                                                                              cgcTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCC‑CGGGCACGGTTTTCGAACTTACCGCCACGCAGGTTTTCCCCGTTTGGCCCCCCGATCACCTTTTCCATCCGTTTTCCGTCAAACAGGGAACCGGGGCGGGCAGTCGGaaaagtt                                                                       >  2:698533/1‑234 (MQ=255)
                                                                                                                                                                        ggTCGCCCCGGAGCGGGGATTGGGGGAGTGTTCGCGCGGGCGCGGGAGGGACATCATGAAGGTGTTTTCAAACCGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCg                                                             <  2:586435/239‑1 (MQ=255)
                                                                                                                                                                        ggTCGCCACGCAGCTGGGGTTGAGGCAGTGTTCGCCCAGGGGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCg                                                             <  2:356018/239‑1 (MQ=255)
                                                                                                                                                                           cGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACCCGGGCGAGCATTAGGAATATGTTTGCCTTCCg‑gg                                                          >  1:294113/1‑237 (MQ=255)
                                                                                                                                                                              cccgcgGCGGGGGTTGGGGCAGTTTTCGCACAGGCGGGGCAGGTACATCATGAAGGGGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCG‑GTgcg                                                       <  2:950144/236‑1 (MQ=255)
                                                                                                                                                                              cACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCG‑GTgcg                                                       <  1:697373/238‑1 (MQ=255)
                                                                                                                                                                                  ccggtGGGATTGGGGGAGTGTTCGCACGGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGCCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCG‑GTGCGCTATg                                                  <  2:419189/235‑1 (MQ=255)
                                                                                                                                                                                                        cGCACAGGCGCGGCAGGTACATCATGAAGGAGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCC‑CGGGCACGGTTTTCGAACTCACCGCCCAGCAGTTCTTCCCCGTTTGGCCCCCAGATCACTTTTTCCATCCCTTTTCCGTCAATCCGTTAACCCGGGCGAGCAGTAGGGATATGTTTTCTTTCCGTGTG‑GCTATGACAATTTTCGTAGGCGGAgggg                            >  2:629651/1‑237 (MQ=255)
                                                                                                                                                                                                         gCACAGGCGCGGCAGGTACCACATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCG‑GTGCGCTATGCAAATGTTCGTAGTCGAAGGTGa                           <  1:208816/239‑1 (MQ=255)
                                                                                                                                                                                                         ccccaGGGGCGGCAGGGTCATCATGAAGGTGTTTTCAAACTGCCCGGTCATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCG‑GTGCGCTATGCAAATGTTCGTAGTCGAAGGTGa                           <  2:633601/236‑1 (MQ=255)
                                                                                                                                                                                                              ggCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGACACCTTTGTCCATCCGTTTTCCGTCCATCCGTGAACCCGGCGGCGCGACAGGAACATGTTTTCCCTTCG‑GGGGGCTCGGCAAAGGTTCTAAGTCTCAGGGGAAGGGt                      >  2:206099/1‑239 (MQ=255)
                                                                                                                                                                                                                cgcgGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAAATCACCGCCCCGCAGTTCTTCCCAGTTTGGCCCCCCGATAACTTTGTCCATCCGTTTGCCGTCAATCCGTGAACGCGGGCTGGCCGTGGGAATATGTTTCCCTTCCG‑GTGCGCTTTCCAAATGTTGCAATTCGAGGGGGAAAGGtac                    >  2:180185/1‑237 (MQ=255)
                                                                                                                                                                                                                  gggCAGGTAAATCATGAAGGTGTTTTCAAACTGCCCGTACTTCGCCTTTTGCATGGCCTAGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCG‑GTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCg                   <  1:576968/237‑1 (MQ=255)
                                                                                                                                                                                                                  cgggtggTCCATGAGGAAGGTGTTTTCAAACTGCCCGTACATTCCCTTTGGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCG‑GTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCg                   <  1:476612/233‑1 (MQ=255)
                                                                                                                                                                                                                  cgGCAGGTACATCATGAAGGTTTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCG‑GTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCg                   <  1:698533/238‑1 (MQ=255)
                                                                                                                                                                                                                  cgGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGGGCAGTAGGAATATGTTTGCCTTCCG‑CTGCGCTATGCAAATGTTCGTAGTCGAAGGGGAACGGTTCgt                  >  1:813942/1‑239 (MQ=255)
                                                                                                                                                                                                                   ggCAGGTACACTAGGAAGGTTTTTCCAAATGGCCCGTACATCTCCTTTGGCATGGCCCCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCG‑GTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCgta                 <  2:620800/239‑1 (MQ=255)
                                                                                                                                                                                                                    gCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGGTTGACTTCCG‑GGTCGCTATGCAAATGTTCGTAGTCGAAGGGGAAAGGTTCggag                >  1:620800/1‑239 (MQ=255)
                                                                                                                                                                                                                      aGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCG‑GTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCgtagt               <  1:276303/238‑1 (MQ=255)
                                                                                                                                                                                                                        tttcatcatGTAGGTTTTTGTAAACTTCCCGTACATCTCCTTTTGCATGGCCCCTAAGTGGCGGTGG‑CGGGCACGTTTTTCGAACTCCCCGCCCTGGAGTGCTTCCCAGTTTGGCGCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCG‑GTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCGTAGTAat            <  2:317628/236‑1 (MQ=255)
                                                                                                                                                                                                                         tACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCCCTTTGTCCATCCGTTTGCCGTCAATCAGTGAACCCGGGCGGCCAGTAGGAATATGTTTGCCTTCCG‑GGGCGCCATGCAAAGGTTCGTAGCCGAAGGGGAAAGGTTCggatgaaat           >  2:63984/1‑231 (MQ=255)
                                                                                                                                                                                                                                  aaGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATTTGTTTGCCTTCCG‑GTTCGCTATGAAAATGTTCGTAGTCGAAGGTTAAAGGGTCGAGGCAATCATAAATCtc  >  2:767297/1‑238 (MQ=255)
                                                                                                                                                                                                                                      tgtTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCG‑GTgcg                                                       <  2:746056/182‑1 (MQ=255)
                                                                                                                                                                                                                                      tgtTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCG‑GTgcg                                                       >  1:746056/1‑182 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
AATTCGTGGGTAACAGAAGATGCATTTTTCTGACTTGCCGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCAATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCG‑CGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCG‑GTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCGTAGTAATCATCAATCTG  >  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: ‑