Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 512,219 C→A 7.5% T193T (ACC→ACA glsA → glutaminase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913512,2190CA 7.5% 676.8 / 31.7 194T193T (ACC→ACAglsAglutaminase 1
Reads supporting (aligned to +/- strand):  ref base C (93/86);  new base A (8/7);  total (101/93)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-01

Only 100 of 202 total aligned reads displayed.

TCGCCTTAGAGTTGCATGGCGGCAAACCGCTTTCGCCACTGGTAAATGCTGGCGCTATTGCCACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTGGCGGTCGTCCCTGGAGTGATGGGAATTGCCGCGTTCTCACCACCGCTGGA  >  NC_000913/511945‑512491
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                
tCGCCTTAGAGTTGCATGGCGGCAAACCGCTTTCGCCACTGGTAAATGCTGGCGCTATTGCCACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATAcc                                                                                                                                                                                                                                                                                  <  2:403317/275‑1 (MQ=255)
   ccTTAGAGTTGCATGGCGGAGAACGGCTTTCGCCACTGGTAAATGCGGGCGCTATTGCCACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATAGCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGt                                                                                                                                                                                                                                                                               <  1:43056/275‑1 (MQ=255)
             gCATGGCGGCAAACCGCTTTCGCCACTGGTAAATGCTGGCGCTATTGCCACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCa                                                                                                                                                                                                                                                                     <  2:122888/275‑1 (MQ=255)
               aaggcggcAAACCCCTTTCGCCACTGGTAAATGCTTGCTCTATTGCAACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGATATCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACg                                                                                                                                                                                                                                                                   <  2:528627/273‑1 (MQ=255)
                        aaacagtttGGCCACTGGTAAAGGCTGGCGGTATTGCCACAACCAGCCTGATTAACGCGGAAAATGGTGAACAACGCTGGCAGCGAATTTTACATATACAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGGCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGTCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCaa                                                                                                                                                                                                                                                           <  4:502937/268‑1 (MQ=255)
                         aCCGCTTTCGCCACTGGTAAATGCTGGCGCTATTGCCACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCa                                                                                                                                                                                                                                                                             <  4:235160/255‑1 (MQ=255)
                             cTTTCGCCACTGGTAAATGCTGGCGCTATTGCCACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACt                                                                                                                                                                                                                                                       >  2:198957/1‑273 (MQ=255)
                                 cGCCACTGGTAAATGCTTGCGCTATTGCCACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGa                                                                                                                                                                                                                                                  <  1:68662/274‑1 (MQ=255)
                                                   gcgcTATTGCCACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATCTCCCTATTGTGATTCAATGGAAGCCTTTGACGCTTATACCCGTCAGTGCTCAACGCTCCCAACTCCCATTGAACTGGACACCCTTGgcgga                                                                                                                                                                                                                               >  4:68534/1‑273 (MQ=255)
                                                         ttGCCACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTCTTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCCCGCTCCTCAATCCTATTGAACTGGCAACGCTTGGCGCGACGCTg                                                                                                                                                                                                                         >  2:356513/1‑275 (MQ=255)
                                                                   caGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCCCTATTGTGATGCAATGGAAACCTGTGACGTGTATACCCGTCAGTGCTCCACGCCCCCCAATACTATTGAACTGGNAACCCNNGGCCCCCCGCCggcgcccggg                                                                                                                                                                                                                >  3:10529/1‑268 (MQ=255)
                                                                                                  cGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATaaa                                                                                                                                                                                        >  2:228052/1‑267 (MQ=255)
                                                                                                                 ttACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGt                                                                                                                                                                 <  2:69550/275‑1 (MQ=255)
                                                                                                                           aaCAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTTGGAATCCGTTGACGCTTAAACGCGTTCTTAAGGCCGACAACGTCCCGTACAtt                                                                                                                                                       >  2:426855/1‑275 (MQ=255)
                                                                                                                               gCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTgg                                                                                                                                                   >  2:257461/1‑275 (MQ=255)
                                                                                                                                  aCTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACCCTGGCTGCCGGTGGTTTGAATCCGTTAACGCATAAACGCGTTCTTCAGGCCGCCACCTGGCCGTCCATTCTGGCCg                                                                                                                                                >  4:375977/1‑275 (MQ=255)
                                                                                                                                                        ctctctGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATa                                                                                                                                                                                                                                                         >  6:1064850/1‑148 (MQ=255)
                                                                                                                                                                  aaGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCa                                                                                                                                                                                                                                                            >  6:597479/1‑135 (MQ=255)
                                                                                                                                                                  aaGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCa                                                                                                                                                                                                                                                            <  5:597479/135‑1 (MQ=255)
                                                                                                                                                                         ccAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                                                                                                                                                     >  5:534026/1‑135 (MQ=255)
                                                                                                                                                                         ccAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                                                                                                                                                     <  6:534026/135‑1 (MQ=255)
                                                                                                                                                                              cGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                                                                                                                                                     <  8:464159/130‑1 (MQ=255)
                                                                                                                                                                              cGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                                                                                                                                                     >  7:464159/1‑130 (MQ=255)
                                                                                                                                                                                 aacaaacaACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACTCTCCTCAATACt                                                                                                                                                                                                                                                       >  6:839960/1‑125 (MQ=255)
                                                                                                                                                                                 aacaaacaACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACTCTCCTCAATACt                                                                                                                                                                                                                                                       <  5:839960/125‑1 (MQ=255)
                                                                                                                                                                                      acaacCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                                                                                                                                                          <  7:164253/117‑1 (MQ=255)
                                                                                                                                                                                      acaacCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                                                                                                                                                          >  8:164253/1‑117 (MQ=255)
                                                                                                                                                                                      acaacCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGc                                                                                                                                                                                                                           >  8:789994/1‑148 (MQ=255)
                                                                                                                                                                                       caaccaacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCt                                                                                                                                                                                                                          <  7:789994/148‑1 (MQ=255)
                                                                                                                                                                                         accaacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTgg                                                                                                                                                                                                                        >  8:115134/1‑148 (MQ=255)
                                                                                                                                                                                           caacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTggcg                                                                                                                                                                                                                      >  7:66394/1‑148 (MQ=255)
                                                                                                                                                                                           caacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTggcg                                                                                                                                                                                                                      >  7:1113955/1‑148 (MQ=255)
                                                                                                                                                                                           caacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTggcg                                                                                                                                                                                                                      >  7:513644/1‑148 (MQ=255)
                                                                                                                                                                                            aacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGctc                                                                                                                                                                                                                                                                <  5:448510/105‑1 (MQ=255)
                                                                                                                                                                                            aacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGctc                                                                                                                                                                                                                                                                >  6:448510/1‑105 (MQ=255)
                                                                                                                                                                                             acTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTggcggc                                                                                                                                                                                                                    >  6:975605/1‑148 (MQ=255)
                                                                                                                                                                                                   aTAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtgg                                                                                                                                                                                                              <  6:12474/148‑1 (MQ=255)
                                                                                                                                                                                                    tAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                                                                                                                                                     >  7:1022435/1‑108 (MQ=255)
                                                                                                                                                                                                    tAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                                                                                                                                                             >  6:556397/1‑148 (MQ=255)
                                                                                                                                                                                                     aaCCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGgtggtg                                                                                                                                                                                                            <  8:1113955/148‑1 (MQ=255)
                                                                                                                                                                                                     aaCCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGgtggtg                                                                                                                                                                                                            <  8:513644/148‑1 (MQ=255)
                                                                                                                                                                                                     aaCCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGgtggtg                                                                                                                                                                                                            <  8:66394/148‑1 (MQ=255)
                                                                                                                                                                                                      aCCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                                                                                                                                                          >  6:1076467/1‑101 (MQ=255)
                                                                                                                                                                                                         gggCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTTTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                                                                                                                                                          <  5:105648/98‑1 (MQ=255)
                                                                                                                                                                                                            ccATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCt                                                                                                                                                                                                                          >  5:331462/1‑127 (MQ=255)
                                                                                                                                                                                                            ccATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCt                                                                                                                                                                                                                          <  6:331462/127‑1 (MQ=255)
                                                                                                                                                                                                             cATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCg                                                                                                                                                                                                    <  7:115134/148‑1 (MQ=255)
                                                                                                                                                                                                              aTAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                                                                                                                                                          >  5:802409/1‑93 (MQ=255)
                                                                                                                                                                                                              aTAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                                                                                                                                                          <  6:802409/93‑1 (MQ=255)
                                                                                                                                                                                                              aTAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCg                                                                                                                                                                                                    >  6:707834/1‑147 (MQ=255)
                                                                                                                                                                                                               tAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGTCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGaa                                                                                                                                  <  4:428112/212‑1 (MQ=255)
                                                                                                                                                                                                                aGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATc                                                                                                                                                                                                      >  8:174159/1‑143 (MQ=255)
                                                                                                                                                                                                                aGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATc                                                                                                                                                                                                      <  7:174159/143‑1 (MQ=255)
                                                                                                                                                                                                                aGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTg                                                                                                                                                                                                 >  7:597514/1‑148 (MQ=255)
                                                                                                                                                                                                                  ccTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                                                                                                                                                             <  7:1047094/134‑1 (MQ=255)
                                                                                                                                                                                                                  ccTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                                                                                                                                                             >  8:1047094/1‑134 (MQ=255)
                                                                                                                                                                                                                             tACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAAcgcg                                                                                                                                                                                    <  5:556397/148‑1 (MQ=255)
                                                                                                                                                                                                                                 ccgccgGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTg                                                                                                                                                                                                                                                   <  6:1093350/81‑1 (MQ=255)
                                                                                                                                                                                                                                 ccgccgGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACg                                                                                                                                                                                              <  6:926337/134‑1 (MQ=255)
                                                                                                                                                                                                                                   gccggtatgctccATTGGGTTGCAATGGAAGCCTGTGACGTGTCTACCCGTGAGTGCTCCACTCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACCCGGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTGGCGgt                                               <  4:198956/266‑1 (MQ=255)
                                                                                                                                                                                                                                   gccgGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGttcttc                                                                                                                                                                              <  5:1082727/148‑1 (MQ=255)
                                                                                                                                                                                                                                   gccgGATATCTCGATTGTGAGTGAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATATTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTGGCGgt                                               <  1:142067/275‑1 (MQ=255)
                                                                                                                                                                                                                                      ggATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAgg                                                                                                                                                                           <  8:597514/148‑1 (MQ=255)
                                                                                                                                                                                                                                       gATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGt                                                                                                                                                                                                   <  8:42920/123‑1 (MQ=255)
                                                                                                                                                                                                                                       gATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGt                                                                                                                                                                                                   >  7:42920/1‑123 (MQ=255)
                                                                                                                                                                                                                                           tctcTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAAcg                                                                                                                                                                                      <  5:97126/132‑1 (MQ=255)
                                                                                                                                                                                                                                             tctATTGGGATGCAATGGAAGCCTGTGACGTTGATACCCGTAAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTGGCGGTCGTCCCTGGa                                     <  1:10529/275‑1 (MQ=255)
                                                                                                                                                                                                                                              ctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                                                                                                                                                                >  5:1125397/1‑119 (MQ=255)
                                                                                                                                                                                                                                              ctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAAc                                                                                                                                                                   >  7:996459/1‑148 (MQ=255)
                                                                                                                                                                                                                                                aTTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGctcc                                                                                                                                                                                                                                                               <  7:519651/54‑1 (MQ=255)
                                                                                                                                                                                                                                                aTTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGctcc                                                                                                                                                                                                                                                               >  8:519651/1‑54 (MQ=255)
                                                                                                                                                                                                                                                aTTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGgtggtg                                                                                                                                                                                                            >  5:1072301/1‑105 (MQ=255)
                                                                                                                                                                                                                                                aTTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                                                                                                                                                                <  8:351538/117‑1 (MQ=255)
                                                                                                                                                                                                                                                aTTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                                                                                                                                                                >  7:351538/1‑117 (MQ=255)
                                                                                                                                                                                                                                                 ttGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTg                                                                                                                                                                >  7:840131/1‑148 (MQ=255)
                                                                                                                                                                                                                                                         cAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGt                                                                                                                                                                                                   >  5:141536/1‑105 (MQ=255)
                                                                                                                                                                                                                                                         cAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACg                                                                                                                                                                  >  8:100927/1‑138 (MQ=255)
                                                                                                                                                                                                                                                           aTGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCCAAAGTGGTGTAGGTGGAGGTCTTCTGGCGGTTGTTCCTGGAGTGAGGGGACTTg                        >  2:434652/1‑274 (MQ=255)
                                                                                                                                                                                                                                                           aTGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTGGCGgtcgtc                                           >  4:151189/1‑255 (MQ=255)
                                                                                                                                                                                                                                                            tGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGATGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCt                                                                                                                                                     <  8:840131/148‑1 (MQ=255)
                                                                                                                                                                                                                                                               aaGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTGGCGGTCGTCCCTg                                       <  2:151190/255‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ccTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGa                                                                                                                                               >  5:475045/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                     gtgACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAat                                                                                                                                            >  5:750643/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                      tgACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGNTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCNAANGCGgnggng                                                                  >  4:570703/1‑218 (MQ=255)
                                                                                                                                                                                                                                                                      tgACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAatg                                                                                                                                           >  6:1045671/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                       gACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCt                                                                                                                                                     <  6:240086/137‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCt                                                                                                                                                     >  5:240086/1‑137 (MQ=255)
                                                                                                                                                                                                                                                                          gtgtATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTgg                                                              <  2:570704/221‑1 (MQ=255)
                                                                                                                                                                                                                                                                           tgtATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGGATCGTGTTGGTTTCCCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTGGCGGTCGTCCCTGGCGTGATGGGAATTGCCGGGTTCTCCCCACCg       >  4:122719/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                             tataCACGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                                                                                                                                                     >  8:357442/1‑35 (MQ=255)
                                                                                                                                                                                                                                                                             tataCACGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                                                                                                                                                     <  7:357442/35‑1 (MQ=255)
                                                                                                                                                                                                                                                                              ataCCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGa                                                                                                                                   <  6:475045/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                               taCCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGaa                                                                                                                                  >  5:658920/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                                aCCCGTCAGTGCTCCACGCTGCTCAGTACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGAc                                                                                                                                                                                               >  7:200373/1‑86 (MQ=255)
                                                                                                                                                                                                                                                                                aCCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGAc                                                                                                                                                                                               <  8:200373/86‑1 (MQ=255)
                                                                                                                                                                                                                                                                                aCCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTGGCGGTCGTCCCTGGAGTGATGGGAATTGCCGCGTTCTCACCCCCgtggga  >  2:331516/1‑270 (MQ=255)
                                                                                                                                                                                                                                                                                 cccGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGgtggtg                                                                                                                                                                                                            <  7:811587/72‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 cccGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGgtggtg                                                                                                                                                                                                            <  5:885075/72‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 cccGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGgtggtg                                                                                                                                                                                                            >  6:885075/1‑72 (MQ=255)
                                                                                                                                                                                                                                                                                 cccGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGgtggtg                                                                                                                                                                                                            >  8:811587/1‑72 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                
TCGCCTTAGAGTTGCATGGCGGCAAACCGCTTTCGCCACTGGTAAATGCTGGCGCTATTGCCACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTGGCGGTCGTCCCTGGAGTGATGGGAATTGCCGCGTTCTCACCACCGCTGGA  >  NC_000913/511945‑512491

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