Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,536,735 A→C 13.9% V192G (GTG→GGG)  narW ← nitrate reductase 2 (NRZ), delta subunit (assembly subunit)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,536,7350AC13.9% 187.9 / 2.3 74V192G (GTG→GGG) narWnitrate reductase 2 (NRZ), delta subunit (assembly subunit)
Reads supporting (aligned to +/- strand):  ref base A (41/19);  new base C (3/11);  total (44/30)
Fisher's exact test for biased strand distribution p-value = 2.03e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGAGGCGCGCCAGGTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATA  >  NC_000913/1536499‑1536973
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
gaggagGCGCGCCAGGTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGCCACACGCCGTGGCGTTATCTTCATTAAACTTCAcc                                                                                                                                                                                                                                              >  1:300985/1‑239 (MQ=255)
     ggtccgcCAGGTGTACTGCCCGTAGTCGTAGAGCAGCCAGCTGCCGAGGAAAAACCCCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGtt                                                                                                                                                                                                                                          <  1:482283/234‑1 (MQ=255)
     ggCGCGCCAAGTGTCTCGCCCGTAGGCGTAGCCCCGACAGCCGGCGAGGAAAAACACCAGCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCCCCTCCCGCCCTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGtt                                                                                                                                                                                                                                          <  2:355738/238‑1 (MQ=255)
     ggCGCCACAGGTGTACTGCCCGTAGTTGTAGCGGCGCCAGGTGCCGAGGAAAAACCCCGTCGCAAAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCACCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGtt                                                                                                                                                                                                                                          <  2:625843/238‑1 (MQ=255)
           ccAGGTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGAAATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCCGCGAGCTGTCACACCCCGTGGGGTTATCTTCAAACAACTTCCCCTGTTCCTCtt                                                                                                                                                                                                                                    >  2:585875/1‑238 (MQ=255)
                      gCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTACCCTGTTCCTCTTCCCAAACTGCa                                                                                                                                                                                                                         >  1:798223/1‑238 (MQ=255)
                      gCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACCTCACCTGTTCCTCTTCCCACACTGCAt                                                                                                                                                                                                                        >  1:114106/1‑239 (MQ=255)
                        cccgttgtgTAGCCCAGCCAGCTGGCAAGGAAAACCCCCCTCGCACAAATGTACGGGTATATATCGGAAAAAAAGCCGGTCAGATCCTGAATCATTTCCCTCCTCCCGCACTGATGTCGACATATCGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTACGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCa                                                                                                                                                                                                                      <  2:444387/231‑1 (MQ=255)
                                 tAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGCCATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACCCTGCATCAAGCGCCTgg                                                                                                                                                                                                             >  1:297144/1‑239 (MQ=255)
                                 tAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGCTATTGATTAAGCGGCGAGCTGTCACACGCCCTGGCGTTATCTTCAATAAACTTCACCTTTTCCTCTTCCCACCCCCCCTCAAGCCCCTgg                                                                                                                                                                                                             >  2:56443/1‑239 (MQ=255)
                                 tAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACCCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAgccgcggg                                                                                                                                                                                                             >  2:433688/1‑233 (MQ=255)
                                        gccagcTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGGGCTGTCACACGCCGGGGGGTTATTTTCAATAAACTTAACCTGTTCCTCTTCCCACCCTGCATCAAGACCCCGGCggggg                                                                                                                                                                                                       >  2:768700/1‑236 (MQ=255)
                                        gccagcTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGGGCCGTTATCTTAAATAAACTTAACCTGTTCCTCTTCCCCCACGCAATCAAGCGCCTgcgggggt                                                                                                                                                                                                       >  2:863045/1‑231 (MQ=255)
                                              tGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGCCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTGCCCTGTTCCCTCTCCCCCACTGCATCAACCCCCTGCCGGGGGGCAAcgc                                                                                                                                                                                                >  2:649360/1‑239 (MQ=255)
                                                      aaaaaCACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATTCTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCtt                                                                                                                                                                                          <  1:357845/237‑1 (MQ=255)
                                                       aaacccccGTTGCAAAAATGTGCGGGGAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCCCCTGTGCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATAGCGCTCTTCg                                                                                                                                                                                        <  2:96499/232‑1 (MQ=255)
                                                               gTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCCCCCCCTTCCCTGttt                                                                                                                                                                                  >  2:146848/1‑236 (MQ=255)
                                                               gTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGGGTCATCGGCATCTTCCCTGTTTACt                                                                                                                                                                               >  1:747577/1‑239 (MQ=255)
                                                                       aaTGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGGCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTtgt                                                                                                                                                                       <  1:146848/239‑1 (MQ=255)
                                                                            aCGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGacac                                                                                                                                                                  <  1:650605/239‑1 (MQ=255)
                                                                             cGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTTCATCAAGCGCCTGCCGGGTGTCATCCCCCCCTTCCCGGTTTCCTTTTTTtgggca                                                                                                                                                                    >  2:143335/1‑232 (MQ=255)
                                                                                 tAGATATCGTAAAAAAATGCGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGGCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTc                                                                                                                                                              <  1:863045/238‑1 (MQ=255)
                                                                                 tAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTc                                                                                                                                                              <  1:56443/238‑1 (MQ=255)
                                                                                 tAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTc                                                                                                                                                              <  1:768700/238‑1 (MQ=255)
                                                                                    atatCGTAAAAAAAGCGCTTCAGATACTGAATCATTTCCCTCCCCCCGCATGGATGTCGAAATATTGCGGCGCGACATCCGGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACt                                                                                                                                                           <  2:482517/238‑1 (MQ=255)
                                                                                     tatCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCCCGCTCTTCGCTGTTTACTTGTTTTTTGAACCTGTCCCTTg                                                                                                                                                         >  1:248314/1‑239 (MQ=255)
                                                                                          tAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTTGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTACCTTGAAAgg                                                                                                                                                    >  1:637698/1‑239 (MQ=255)
                                                                                          tAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCATCTTCCCACACTGCATCAAGCGCCTGCCGGGGGGCATCGCGCTCTTCGCGGTTTCCTTGTTTGTTGGTACTGTCACTTGCAAgg                                                                                                                                                    >  1:121509/1‑239 (MQ=255)
                                                                                                              cTGAATCATTTCCCGCCTCCCGACCTGAGGTCGACATTTTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCCCCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTgcag                                                                                                                                <  1:711108/239‑1 (MQ=255)
                                                                                                               tGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCAAACTGCATCAATCGCCTGGCGGGTGTCATCGCGCACTTCACTGTTTACTTGTTTTTTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTGgcccc                                                                                                                               >  1:526802/1‑236 (MQ=255)
                                                                                                               cgAATCTTTTCCCCCCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGGGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTgcagca                                                                                                                              <  1:612538/238‑1 (MQ=255)
                                                                                                                gAACCATTTCCCCCCCCCCGCCCTGATGCCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCTCTGTTATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTgcagca                                                                                                                              <  1:756382/239‑1 (MQ=255)
                                                                                                                 aaTCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGcagcag                                                                                                                             >  1:53205/1‑239 (MQ=255)
                                                                                                                 aaTCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTAACTTGAAAGGGAGCTTCCCGCCAGTTGcagcag                                                                                                                             >  1:167487/1‑239 (MQ=255)
                                                                                                                   tCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGGGCTTCCCGCCAGTTGcagcagg                                                                                                                            >  1:217432/1‑237 (MQ=255)
                                                                                                                   ccttttCCCCCCCCCCGCTCTGTTTTCTATATTTTGCGGCCTACCATCCTGGCTAACGGGTGTGTGATATTGATTAAGCGGCGAGCTGTCCACCGCCGGGGCTTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCaga                                                                                                                            <  2:297144/235‑1 (MQ=255)
                                                                                                                     aTTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGGGTCATCGCGCTCTTCGCTGTTTACTTGTTTTTTGACCCCGCACCTTGAAAGGGGACTTCCCCCCCGtttccacaaggat                                                                                                                         >  2:711108/1‑228 (MQ=255)
                                                                                                                      tttCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTTTGACCCTTTCACTTGAAAGGGAGCTTCCCGCCAGTTgcaccaaaacac                                                                                                                        >  2:144688/1‑233 (MQ=255)
                                                                                                                                         aTGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTTTGCCACTGTCACTTGAAAGAGAGCTTCCCGCCAGTTGCAGCAGACAGTCAAACACAGCGTACCACgg                                                                                                     >  1:97498/1‑239 (MQ=255)
                                                                                                                                         aTGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTTACCCTTTCCCTTGCAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTCCCACgg                                                                                                     >  1:433053/1‑239 (MQ=255)
                                                                                                                                           gTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTTCAGCAGAGCATCAAACAACGCGTACCACGGt                                                                                                    >  1:310036/1‑238 (MQ=255)
                                                                                                                                                   atgtgggccCGAAATCCGGGCTAAAGCGTGGCTGATTTTGATTAGGCGGCCAGCTTTCCCAGTCCGGGGCGTTATCTTCAATAAACTTCCCGTGTTCCTCTACCCTCCCGTCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGGGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgt                                                                                           <  2:637698/234‑1 (MQ=255)
                                                                                                                                                   atTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgt                                                                                           <  1:610299/239‑1 (MQ=255)
                                                                                                                                                   atTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGGGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGTGGCCTCGCgt                                                                                           >  1:493599/1‑239 (MQ=255)
                                                                                                                                                     tGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCCGTTGCAGCAGAGCATCAAACAAAGCGTACCACGGGGCCTCGCgttg                                                                                         >  1:44138/1‑239 (MQ=255)
                                                                                                                                                                                cGCTGATATTGCTTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCCCTCTTAGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTGCCCGCCAGTTGCACCAGAGCAGCAAAAAAACCGTCACAAGGGGCCCCCCGGTGTTTTAAGCGCCCGCCAAAAAGGGCggg                                                              >  2:403962/1‑237 (MQ=255)
                                                                                                                                                                                cGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGCAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACCCGTAACACGGTGCCTCGCGTTGTTTTAAGCGACCCCCAACCAGGGCggg                                                              >  2:725150/1‑237 (MQ=255)
                                                                                                                                                                                cGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGGGCCTCGCGTTGTTTTAAGGGACGGCCAAGCAGGGCGAg                                                              >  1:203447/1‑239 (MQ=255)
                                                                                                                                                                                    gATTTTGATTAAGCGGCGAGCTGTCCCGCGCCGGGGCGTTAACCTCAATAAACTTCCCCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGAt                                                           <  1:433688/238‑1 (MQ=255)
                                                                                                                                                                                       atTGATTAAGCGGCGAGCTGTCCCACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTCCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGgc                                                       <  1:97106/239‑1 (MQ=255)
                                                                                                                                                                                        ttGATTAAGCGGCGAGCTGTCACACGCCGTGGCTTTATCTTCAATAAACTTCCCCTGTTCCTCTTCCCACCCTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGgcg                                                      <  1:725150/239‑1 (MQ=255)
                                                                                                                                                                                        gtgtttAACCGGGGTGGTGTCCCACGCCGTGGGGTTATTTTCAATAAACTTCCCCTGTTCCGCTTCCCACCCTGCTTCAAGCGCCTGGGGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGGGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGgcg                                                      <  2:53205/235‑1 (MQ=255)
                                                                                                                                                                                               aGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAAAGCGTACCCCGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGCCGCAACCtt                                               >  1:514087/1‑239 (MQ=255)
                                                                                                                                                                                                   ggcgccTGTCCCACCCCGTGGGGTTTCCTCAAAAAACCTTCACCTGTTCCCTTTCCCCCACTGCAGCAAGCGCAGGGCGGGTGTCATCGCGGTCTTCGTTGTTTTCTTGTTTTGTGACACTGTCGCTTGAAAGGGAGCGTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGc                                           <  2:805193/234‑1 (MQ=255)
                                                                                                                                                                                                    cgggCTGCCCCACGACGTGGCGTTATCTTAAAAAAACTTCCCCTTTTCCTCTCCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACCCTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGc                                           <  1:138537/235‑1 (MQ=255)
                                                                                                                                                                                                    cgAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCCGTTGCAGCAGAGCATCAAACAACGCGTCCCACGGGGCCTCGCGTTGTTTTAAGCGACCGCCAACAAGGGCGGGGGTCGGCGGAAAATTGTGCa                                          >  2:745577/1‑239 (MQ=255)
                                                                                                                                                                                                     gAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACAGTGCCTCGCGTTGTTTTAAGCGACCGCCACGCAGGGCGAGGGGCGGCGCAAAATGGGGCaa                                         >  2:479343/1‑239 (MQ=255)
                                                                                                                                                                                                     gAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACCCGTACCACGGTGCCTCGCCTTGTTTTAACCGACCGCCAAACCGGGCGAAGAACGGCGCAAAATTGAGCa                                          >  2:351182/1‑238 (MQ=255)
                                                                                                                                                                                                                   cggtgtgTTTTCTGCAATAACCTCCCCCTGTCCCTCTTCCTCCACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTTTGTTGTGACACTGTCCCTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTGCCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCt                           <  2:273419/233‑1 (MQ=255)
                                                                                                                                                                                                                   cGTGGCTTTTTCTTGAAATACTTTACCCTGCTCCTCTTCCCCCCCTGCATCAAGCCACTGAGGGGTGCGCTGGGTCTTTTCGCGGTTTACTTGTTTTGTGACACTGACACTTGAAAGGGAGTGTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCt                           <  2:741210/239‑1 (MQ=255)
                                                                                                                                                                                                                   cGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCCGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCg                              >  1:464219/1‑236 (MQ=255)
                                                                                                                                                                                                                   cGTGGCGTTATCCTCAATAAACTTCCCCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTTTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCt                           <  1:21462/239‑1 (MQ=255)
                                                                                                                                                                                                                     ttgCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTAATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGGGAGGATCGGGGCAACATTGAGCAATCCTTCTTTCCCCTgg                         >  2:756366/3‑238 (MQ=255)
                                                                                                                                                                                                                     tGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTTCAGCAGAGCATCAAACAACGCGTCCCACGGTTCCTCGCGTTGTTTTAAGCGACCGCCCAACAGGGCTGGGATCGGCGCAACCTTGGGCAATCCTTTTTTcccccac                         >  2:806246/1‑236 (MQ=255)
                                                                                                                                                                                                                     tGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCCCGTTGTTTTAAGCGACCGCCAAGCAGGGCGGGGATCGGGGCAACATTTGGCAAACCTTCTTTTGCCTg                          >  2:378232/1‑238 (MQ=255)
                                                                                                                                                                                                                             tCTTCAATAAACTTACCCTGTTCCTCTTCCCACACTTCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCTGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAgg                 <  1:585875/239‑1 (MQ=255)
                                                                                                                                                                                                                               ttCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGAATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGGCGGGATCGGCCCAACATTGAGCAATCCTTCTTTCGCCTGATCggcaggca               >  2:115121/1‑239 (MQ=255)
                                                                                                                                                                                                                                  aaTAACCTTCCCCTGTTCCCTTTCCCACCCTGCATCAAGGGCCTGGCGGGTGGCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTa            <  2:479776/239‑1 (MQ=255)
                                                                                                                                                                                                                                  aaTAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGCGCTTCCCGCCAGTTGCAGCAGAGCATCACACAACGCGTACCACGGTGCCTCGCGTTGTTTTAATCGACCGTCAAGCAGGGAGAGGCTCGGCGCAACATTGAGCAATCCTCCTTCCGCCTGATCGCCCGGCAgca            >  2:602283/1‑237 (MQ=255)
                                                                                                                                                                                                                                  aaTAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTTAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGt             >  1:54737/1‑238 (MQ=255)
                                                                                                                                                                                                                                       aCTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCCCGTTGTTTTAAGCGCCCGCCAAGCAGGGCGAGGATCGGCGCAACATTGGGCAATCCTTCTTTCGCCTTATCGTCAGGCAGGACGCtt       >  2:387922/1‑239 (MQ=255)
                                                                                                                                                                                                                                           tCACCTTTTCCCCTTCCCACACTGAATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTaaa    <  1:144688/239‑1 (MQ=255)
                                                                                                                                                                                                                                            cACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAAt   >  1:479776/1‑239 (MQ=255)
                                                                                                                                                                                                                                            cACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGAACGCTTaaaa   >  2:494366/1‑238 (MQ=255)
                                                                                                                                                                                                                                             aCCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCCACATTGAGCAATCCTTCTTTCCCCTGGTCGTCAGGCAGTAAGCTTaaaaa  >  2:746035/1‑237 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
GAGGAGGCGCGCCAGGTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCCACCTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATA  >  NC_000913/1536499‑1536973

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