Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,030,269 C→T 14.9% S69F (TCT→TTT)  yccX → weak acylphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,030,2690CT14.9% 86.4 / 2.1 63S69F (TCT→TTT) yccXweak acylphosphatase
Reads supporting (aligned to +/- strand):  ref base C (19/33);  new base T (7/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 3.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCCATTAAGGATAAGGAGAGCAGCAATGTCGAAAGTCTGCATAATTGCCTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAG  >  NC_000913/1030037‑1030490
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                  
cgTCCATTAAGGATAAGGAGAGCAGCAATGTCGAAAGTCTGCATAATTGCCTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGGGGTTAAAGGAGGGCGGGCCCCGCTTTgcggt                                                                                                                                                                                                                              >  2:333084/1‑237 (MQ=255)
       taaggataaggaGAGCAGCAATGTCGAAAGTCTGCATAATTGCCTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGaa                                                                                                                                                                                                                      <  1:47399/239‑1 (MQ=255)
           gataaggaGAGCAGCAATGTCGAAAGTCTGCATAATTGCCTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGGTCTGCGCCTGGAGAACcgg                                                                                                                                                                                                                  >  1:894266/1‑239 (MQ=255)
              aaggaGAGCAGCAATGTCGAAAGTCTGCATAATTGCCTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCCCGGTTTGCGGCTGTGGAACggggg                                                                                                                                                                                                                >  1:1123119/1‑236 (MQ=255)
               aggaGAGCAGCAATGTCGAAAGTCTGCATAATTGCCTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGGAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGGTCTGGGGGTGTGGAACggggggt                                                                                                                                                                                                              >  1:127419/1‑235 (MQ=255)
               aggaGAGCAGCAATGTCGAAAGTCTGCATAATTGCCTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGGGAAGAAGGGCAGGGGGAAAAATTAATGCAGTGGTTAAAGGGGGGGGGGCCCCGGTCTTCGGgtgtgacacgggggga                                                                                                                                                                                                               >  2:1195444/1‑227 (MQ=255)
                    agcagcAATGTCGAAAGTCTGCATAATTGCCTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGCCGATGGCAGCGTGGAAGTGGTTGCTTGCGGGTGAGAAGGGCAGGTGGAAAAATTAATGCATGGGTTAAAGCGTGGCAGGCGGCGGTCTgccggggggaaaaaggggcgctcag                                                                                                                                                                                                          >  2:266424/1‑216 (MQ=255)
                    agcagcAATGTCGAAAGTCTGCATAATTGCCTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGGGGTTAAAAAGTGGCGGCCCCCGGTCTCGGGGTGTTGAACgggggggcgggg                                                                                                                                                                                                         >  2:859797/1‑230 (MQ=255)
                     gcagcaATGTCGAAAGTCTGCATAATTGCCTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGGTCTGcgcggggtaaacggggggtcacagg                                                                                                                                                                                                         >  2:749315/1‑218 (MQ=255)
                     gcagcaATGTCGAAAGTCTGCATAATTGCCTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTCACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATCCAGTGGTTAAAGAGGGGCGGACAGGGGGATGCGCGCGGGGAAAGGGGGCTAAGCGa                                                                                                                                                                                                        >  2:1183921/1‑239 (MQ=255)
                       ataaaTGTCGAAAGTTGGCATAATTCCCGGGGTTTACGGGGGATGTCAGGGCGTAGGATTGCGCTACCGCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAAGGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGAg                                                                                                                                                                                                       <  2:563506/235‑1 (MQ=255)
                         caATGTCGAAAGTCTGCATAATTGCCTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAACGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGGGGTTAAAGAGTGGCGGCCCCGGTGCTGGGCGCGTGGaaagaggggcgcacgggcg                                                                                                                                                                                                     >  2:76248/1‑221 (MQ=255)
                           aTGTCGAAAGTCTGCATAATTGCCTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGGGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCCCGATTTGCGCGCGTAGAACGGGGGCGACGCGAGCCGc                                                                                                                                                                                                  >  1:845685/1‑238 (MQ=255)
                                 aaaGTCTGCATAATTGCCTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGGGGCGGCCCGCGGTGTGGGCGGGTAGAACAGGGGGCAAGCGAGACGcgcaacc                                                                                                                                                                                            >  2:421773/1‑232 (MQ=255)
                                     tCTGCATAATTGCCTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCATGGGTTAAAGGGTGGCGGGCCCGCGTTTTGGGCGTGTAGACGGGGGGCTAAGAGACGCGCATCAccccgc                                                                                                                                                                                        >  2:340488/1‑237 (MQ=255)
                                       ggaaaaaatGCGTGGGTTTAGGGGGGTGTGCGGGGCGTAGGTTTCCGCAACCCCACACAGGACGAAGCGAAAAGACTGGGTTTCACCGGGTACGCCAAAAATTTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGGGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCggg                                                                                                                                                                                     <  2:845685/231‑1 (MQ=255)
                                        gCATAATTGCCTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCCCGTTCTGCCGGGGTGGAACGGGGGTTCCGCGGGCCGCGTCCTCCCCCgggg                                                                                                                                                                                    >  1:475914/1‑239 (MQ=255)
                                              ttGCCTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGGCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTa                                                                                                                                                                              <  1:1183921/239‑1 (MQ=255)
                                                 ccTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTGGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGGTTTGCGCGTGTAGAACGGGGGGTCAGCGAGCCCACTCACTCCCCGGGGGGATTTAc                                                                                                                                                                            >  1:906926/1‑238 (MQ=255)
                                                   tGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGGTCTGGCCGTGTAGAACGGGGGCTGAGCGGGCCGGCGCCGCCCACGGGGGGATTAAACgg                                                                                                                                                                         >  1:338548/1‑238 (MQ=255)
                                                      gTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGGTTTGCGCGTGTAGAACGGGGGCTCAGCGAGCCCCACCATCCCTCGGGGGGATTAACCGAttt                                                                                                                                                                      >  1:335284/1‑239 (MQ=255)
                                                       tttACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAAc                                                                                                                                                                                                                     <  1:1220194/192‑1 (MQ=255)
                                                            gggCGTGTTCGGGGCGAAGGATTTGGCAACCCCCACCAGGCAGAGTCAAAAAGCATGGGTTTCAACGGGTACGCCAAAAATCTTGGCGGTGGAAGCGTGGACGTGGTTGCTTGGGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTAt                                                                                                                                                                <  2:846193/239‑1 (MQ=255)
                                                            gggCGAGGTAAGGGCGTAGGATTTCGCTTCACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGGGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGGGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTAt                                                                                                                                                                <  1:266424/239‑1 (MQ=255)
                                                                ggtttCGGGGCTTGGGGTTTCGCTTCCCCCCAAAGTAAGGAGCGAAAAGACTGGGTTTAACGGGGTCCGCAAAAAATGTTGCAGGTGGCAGGGTGGAAGGGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGGGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCg                                                                                                                                                             <  2:245560/235‑1 (MQ=255)
                                                                                       tacaccgaccaGTACGGAGCAAAAAGGCTGGGTTTAACCGGGTACGCCAAAAATCTTGACGTTGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGGGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTgg                                                                                                                                   <  1:1135153/230‑1 (MQ=255)
                                                                                       tACACCAAAAAGGACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAATTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTgg                                                                                                                                   <  2:338548/239‑1 (MQ=255)
                                                                                        acacCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGc                                                                                                                                  <  1:421773/239‑1 (MQ=255)
                                                                                              cacaGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGGGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCa                                                                                                                           <  1:1276321/239‑1 (MQ=255)
                                                                                              cacaGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGGGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTTACCGGTTTTCGTATCCGTTAAATACATTTTAACCGTTTTTGCGGGcgc                                                                                                                           >  1:1317053/1‑237 (MQ=255)
                                                                                              cacaGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGGTCTGCGCGCGTAGAACGGGTGCTCAGCGAGCCGACTCATCCCTCGGGGGGATTAACCGATTTTTGTATCCGTTAAATACATTTTACCCGTTTTTGCAGGccc                                                                                                                           >  1:512785/1‑238 (MQ=255)
                                                                                              cacaGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGCGCTGCGCGGGGAGAAAAGGGGCTCAGCGAGGCCCCGCCTCCCCCGGGGGACGTTACAGATTTTTGTAGCCCGTAAATAAAAATTACCGGTTTTGGGAgggcg                                                                                                                           >  2:624984/1‑236 (MQ=255)
                                                                                                   tACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCCCGGTCTTCGCGTGTAGAACGGGGGCTCCGCGAGCCGCCACATCCCTCGGGGGGATTAACCGGTTTTTGTATCCCTTAAATAAATTTTTACTGTTTTTGCAGGccacca                                                                                                                        >  1:368106/1‑237 (MQ=255)
                                                                                                   tACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGGGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCCCGCGCTCGGCGGGCAGGACGGGTGCTCAGCGGCGCCCACCCTACCTCCGGGGATTTAACCGGTGTTTTGATCCCTTACAATAATTTTACCGGTGTTTggtgggctacaca                                                                                                                        >  2:963989/1‑228 (MQ=255)
                                                                                                             aaaGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGGCCGCCGTTTTCGCCGGTAGAACCGGTGGTCAGCGGGCCGCATACTCCCCCCGGGGGATTAAACGGTTTTTTTATCCGTTAAATACATTTTAACCGTTTTTGGAAGCCAGCAAatttttggggt                                                                                                            >  2:591246/1‑232 (MQ=255)
                                                                                                             aaaGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGGGTTGAAAGGGTGGTCCGCGGGCGCCCTCCACCCTCGGGGGGATTTACCCCATTTGCTATCCCCTAGATAACATTTACCGGGTTTggatgggccccacattttgtggct                                                                                                             >  2:1247163/1‑217 (MQ=255)
                                                                                                                               ggTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGGTCTGCGCGTGTAGAACGGGGGCTCAGCGAGCCGCATCATCCCTCGGGGGGATTAACCGATTTTTGTATCCGTTAAATACATTTTACCGGTTTTGGCAAGGCCAGCAAATTTGGGGGCTTGCTTAGCCGGGCCtt                                                                                            >  1:536407/1‑238 (MQ=255)
                                                                                                                                        aaaaTTTTGACGGTGGCGGGTGGGAAGTGGTTGCTGGGGGGGAAGAAGGGCGGGGGAAAAAATTAAGCCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGGGTGTAAAACGGTTGCTCAGCGAGCCGCATCATCCCTCGGGTGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGt                                                                                 <  2:906926/239‑1 (MQ=255)
                                                                                                                                        aaaaTCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGt                                                                                 <  1:859797/239‑1 (MQ=255)
                                                                                                                                        aaaaTCTGGACGGGGGCGGCGTGGAAGTGGTTGCTTGCGGTGAAGAGGGGCAGGTGGAAAAATTAATGCAGGGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGt                                                                                 <  1:340488/239‑1 (MQ=255)
                                                                                                                                         aattttgagcGATGGAAGGGGGGAGGTGTTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATGAAGGCAGTGGTAAAAGAGGGGCGGCCCGGGTTCTGCGCGTGTGGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTc                                                                                <  2:368106/230‑1 (MQ=255)
                                                                                                                                         aattcttggctttggactaggggAAGTGGTTGGGTGTGGTGAAGAAGGGACGGTGGAAAAATTAATGTGGGGGTAAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGTTGCTCCGCGAGCCGTATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTc                                                                                <  2:823460/218‑1 (MQ=255)
                                                                                                                                         aaattgtGACGGTGGCAGGGTGGAAGTGGTTGCTTGGGGGGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCGGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTc                                                                                <  1:591246/233‑1 (MQ=255)
                                                                                                                                         aaaTTTTGGCGGTGGACGGGTGGAGGTGGTTGTTTGGGGTGAAGAAGGGCAGGTGGAAAAATTAAGGAATGGGTTAAAGAGTGGCGGCCCGCTTTCTGGGGGTTTAGACCGGGGGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTc                                                                                <  2:335284/239‑1 (MQ=255)
                                                                                                                                         aaaTTTTGACGGTGGGAGCGTTGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAAGGCAGGTGTTAAAGAGGGGCGGCCCGTGTTCTGGGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTc                                                                                <  1:624984/239‑1 (MQ=255)
                                                                                                                                         aaaTCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTAAAGAAGGGCAGGTGGAAAAATTAAGGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTc                                                                                <  1:907540/239‑1 (MQ=255)
                                                                                                                                          tgtgtagaagagaccggcttGGAAGTGGTTGCTTGCGGTGAAAAAGGGCAGGTGGAAAAATTAATTCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTc                                                                                <  1:988285/219‑1 (MQ=255)
                                                                                                                                                     tggaggggTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCATTCTGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAg                                                                     <  1:1195444/231‑1 (MQ=255)
                                                                                                                                                      tgaggggTGGAGGGGTTTGGTTGGGGTGAAGAAGGGCAGGTGGAAAATTTAATGCAGTGGTTAAAGAGGGGCGGCCCGCGTTCTGCGCGTGTAAAACGGGTGCTAAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGc                                                                    <  2:61223/232‑1 (MQ=255)
                                                                                                                                                                     gTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCCGTCTGCGCGTGTAGGAACGGGGCTCAACGAGCGGCAACATCCCCCCGGGGGGTTAACCGAGTTTCGTAGTCGGTAAATACATTTTAACGGGTTTGGGCGGACCGCAATTTTTGGGGGTTTGTTTGCCGGGGCTTTTGGTAAAAAGCGGGATAAATAACgttcctttttataaa                                                     >  2:988285/1‑224 (MQ=255)
                                                                                                                                                                                tGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCGCAGCGGGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTa                                                                                                                                              <  1:1230595/140‑1 (MQ=255)
                                                                                                                                                                                  aagaagGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTAt                                                                                                                                                                <  2:314652/121‑1 (MQ=255)
                                                                                                                                                                                  aagaagGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGAGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTAt                                                                                                                                                                >  1:314652/1‑121 (MQ=255)
                                                                                                                                                                                          gaggtgaaaaaattatgtgtggtgTTAAAGGGGGGCGGCCCGGTTTCTGGGGGTGGAGAACGGGTGCCCAGGAGTCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCa                                <  2:1317053/215‑1 (MQ=255)
                                                                                                                                                                                          cgggTGGAAAATTAAGTGAAGGGGTTAAGGGTGTGCGGCCCGGTTTCTGCGGGTTTAGAACCGGGCGTCGTCGTGTCGGATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCa                                <  2:892202/236‑1 (MQ=255)
                                                                                                                                                                                                     aTTAAGGCATGGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCttt                         <  1:679387/234‑1 (MQ=255)
                                                                                                                                                                                                          tGCAGTGGTTAAAGAGGGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCCCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCa                                <  1:1279688/222‑1 (MQ=255)
                                                                                                                                                                                                               tGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGGGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCACTTTTGGTGGCTTTCTTTGCCGGCCCTTTCTGTAAACGACGGTGTAAAAAACGGCGGTTACCCTTTTTTTTTCCCAAATAATATTCCATCGCCTTTACCAGcccaactc           >  2:679387/1‑231 (MQ=255)
                                                                                                                                                                                                                 gTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTc                                                                                <  1:683677/167‑1 (MQ=255)
                                                                                                                                                                                                                    aaaGAGGGGGGGCCCGCGTTCTGCGCGGGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCg     <  1:749315/239‑1 (MQ=255)
                                                                                                                                                                                                                         ggggcggcCCGCGTTTGGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGgag  <  1:76248/235‑1 (MQ=255)
                                                                                                                                                                                                                         gTGGCGGCCCGTGTTTCGCGCGTGTAGAACGAGTGCCCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTTGCTTGCTTAGCCGGCCCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGgag  <  1:150513/237‑1 (MQ=255)
                                                                                                                                                                                                                         gTGGCGGCCCGCGTTCTGGGCGTGTAGAACGGGTGCTAAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGgag  <  1:737342/237‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                  
CGTCCATTAAGGATAAGGAGAGCAGCAATGTCGAAAGTCTGCATAATTGCCTGGGTTTACGGGCGAGTTCAGGGCGTAGGATTTCGCTACACCACACAGTACGAAGCGAAAAGACTGGGTTTAACCGGGTACGCCAAAAATCTTGACGATGGCAGCGTGGAAGTGGTTGCTTGCGGTGAAGAAGGGCAGGTGGAAAAATTAATGCAGTGGTTAAAGAGTGGCGGCCCGCGTTCTGCGCGTGTAGAACGGGTGCTCAGCGAGCCGCATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAG  >  NC_000913/1030037‑1030490

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