Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 512,212 T→G 6.3% V191G (GTG→GGG)  glsA → glutaminase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913512,2120TG 6.3% 624.5 / 28.5 181V191G (GTG→GGG) glsAglutaminase 1
Reads supporting (aligned to +/- strand):  ref base T (88/80);  new base G (6/6);  total (95/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

Only 100 of 194 total aligned reads displayed.

CAGTTATCGCCTTAGAGTTGCATGGCGGCAAACCGCTTTCGCCACTGGTAAATGCTGGCGCTATTGCCACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTGGCGGTCGTCCCTGGAGTGATGGGAATTGCCGCGTTCTCACCACCG  >  NC_000913/511939‑512486
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                  
      tCGCCTTAGAGTTGCATGGCGGCAAACCGCTTTCGCCACTGGTAAATGCTGGCGCTATTGCCACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATAcc                                                                                                                                                                                                                                                                             <  2:403317/275‑1 (MQ=255)
             agagTTGCATGGCGGCAAACCGCTTTCGCCACTGGTAAATGCTGGCGCTATTGCCACCACCAGCCTGATCAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCAATAGCCTGGCTGCTGTAATCCGCCGGATATCGCTATTGAGATGCAATGGAAGCCTGTGACGTGTATACCCgctgg                                                                                                                                                                                                                                                                       >  1:57168/1‑270 (MQ=255)
                    cATGGCGGCAAACCGCTTTCGCCACTGGTAAATGCTGGCGCTATTGCCACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCAc                                                                                                                                                                                                                                                               >  2:214162/1‑275 (MQ=255)
                           gcAAACCGCTTTCGCCACTGGTAAATGCTGGCGCTATTGCCACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATAccc                                                                                                                                                                                                                                                                            >  2:235161/1‑255 (MQ=255)
                               aCCGCTTTCGCCACTGGTAAATGCTGGCGCTATTGCCACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCa                                                                                                                                                                                                                                                                        <  4:235160/255‑1 (MQ=255)
                                                              aTTGCCACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATACCCGGGCCATAGCCTGGCGGCTGTACTCCGCCGACTCACTCTATTGGGCTCCACGGGAGGCCTGTTACTTCTATACCCGTCAGTGCTCCAAGCGCCTCAATAATATTGAACTGGCAACCCtttgggcaccccg                                                                                                                                                                                                                     >  4:505878/1‑264 (MQ=255)
                                                               ttGCCACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTCTTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCCCGCTCCTCAATCCTATTGAACTGGCAACGCTTGGCGCGACGCTg                                                                                                                                                                                                                    >  2:356513/1‑275 (MQ=255)
                                                                         caGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCCCTATTGTGATGCAATGGAAACCTGTGACGTGTATACCCGTCAGTGCTCCACGCCCCCCAATACTATTGAACTGGNAACCCNNGGCCCCCCGCCggcgcccggg                                                                                                                                                                                                           >  3:10529/1‑268 (MQ=255)
                                                                            ccTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGCTGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCCTGGCGAGACGCTGGCGGCAggtggt                                                                                                                                                                                                        >  2:59905/1‑274 (MQ=255)
                                                                                                                       ttACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGt                                                                                                                                                            <  2:69550/275‑1 (MQ=255)
                                                                                                                                 aaCAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTTGGAATCCGTTGACGCTTAAACGCGTTCTTAAGGCCGACAACGTCCCGTACAtt                                                                                                                                                  >  2:426855/1‑275 (MQ=255)
                                                                                                                                        aCTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCg                                                                                                                                           >  2:141570/1‑275 (MQ=255)
                                                                                                                                        aCTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACCCTGGCTGCCGGTGGTTTGAATCCGTTAACGCATAAACGCGTTCTTCAGGCCGCCACCTGGCCGTCCATTCTGGCCg                                                                                                                                           >  4:375977/1‑275 (MQ=255)
                                                                                                                                           ggcggggtggCACGTAGCGCCTTCTGACGAAGCACACCAAGCGGAACAAACAACCAACTTCCGGACCCGTGCCATTGCCTGGCTGCTGTGTGCCGCCGGATATCGCTATTGTGAAGCAATGGAAGCCTGTGACGGGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGACGGCAGGTGGTGTGAATCCGGTGACGCATAAACGCGTTCGGCAGGCCGACAACGTGCCGTACATTCTGGCCGaga                                                                                                                                        <  4:59904/266‑3 (MQ=255)
                                                                                                                                                                      cGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACgtg                                                                                                                                                                                                                                                                                   <  7:10590/109‑1 (MQ=255)
                                                                                                                                                                        aaGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCa                                                                                                                                                                                                                                                       >  6:597479/1‑135 (MQ=255)
                                                                                                                                                                        aaGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCa                                                                                                                                                                                                                                                       <  5:597479/135‑1 (MQ=255)
                                                                                                                                                                         aGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGGATGGTCGCTcc                                                                                                           >  2:283036/1‑274 (MQ=255)
                                                                                                                                                                            cAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCa                                                                                                                                                                                                                                      >  7:172297/1‑148 (MQ=255)
                                                                                                                                                                             aaCCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCg                                                                                                                                                        <  2:178288/225‑1 (MQ=255)
                                                                                                                                                                               ccAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                                                                                                                                                >  5:534026/1‑135 (MQ=255)
                                                                                                                                                                               ccAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                                                                                                                                                <  6:534026/135‑1 (MQ=255)
                                                                                                                                                                                cAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                                                                                                                                                <  8:18409/134‑1 (MQ=255)
                                                                                                                                                                                cAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                                                                                                                                                >  7:18409/1‑134 (MQ=255)
                                                                                                                                                                                 aGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCGGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTgg                                                                                                  <  2:53543/275‑1 (MQ=255)
                                                                                                                                                                                   tCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTg                                                                                                                                                                                                                               >  5:12474/1‑148 (MQ=255)
                                                                                                                                                                                       aacaaacaCCCCACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCCGGTATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTa                                                                                             <  2:505879/274‑1 (MQ=255)
                                                                                                                                                                                       aacaaacaACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACTCTCCTCAATACt                                                                                                                                                                                                                                                  >  6:839960/1‑125 (MQ=255)
                                                                                                                                                                                            acaacCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                                                                                                                                                     <  7:50261/117‑1 (MQ=255)
                                                                                                                                                                                            acaacCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                                                                                                                                                     >  8:50261/1‑117 (MQ=255)
                                                                                                                                                                                               accaacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGa                                                                                                                                                                                                                                             <  5:759832/122‑1 (MQ=255)
                                                                                                                                                                                                  aacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACgtg                                                                                                                                                                                                                                                                                   <  5:1016263/81‑1 (MQ=255)
                                                                                                                                                                                                  aacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACgtg                                                                                                                                                                                                                                                                                   >  6:1016263/1‑81 (MQ=255)
                                                                                                                                                                                                  aacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGctc                                                                                                                                                                                                                                                           <  5:448510/105‑1 (MQ=255)
                                                                                                                                                                                                  aacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTggcgg                                                                                                                                                                                                                <  5:1064850/148‑1 (MQ=255)
                                                                                                                                                                                                         aTAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGa                                                                                                                                                                                                                                             <  8:795510/112‑1 (MQ=255)
                                                                                                                                                                                                         aTAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGa                                                                                                                                                                                                                                             >  7:795510/1‑112 (MQ=255)
                                                                                                                                                                                                         aTAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGgcg                                                                                                                                                                                                                            >  7:430752/1‑129 (MQ=255)
                                                                                                                                                                                                         aTAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGgcg                                                                                                                                                                                                                            <  8:430752/129‑1 (MQ=255)
                                                                                                                                                                                                         aTAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtgg                                                                                                                                                                                                         <  6:12474/148‑1 (MQ=255)
                                                                                                                                                                                                          tAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                                                                                                                                                <  8:1022435/108‑1 (MQ=255)
                                                                                                                                                                                                          tAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                                                                                                                                                        >  6:556397/1‑148 (MQ=255)
                                                                                                                                                                                                          tAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                                                                                                                                                        >  7:156644/1‑148 (MQ=255)
                                                                                                                                                                                                          tAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                                                                                                                                                        <  8:172297/148‑1 (MQ=255)
                                                                                                                                                                                                            aCCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                                                                                                                                                     <  5:1076467/101‑1 (MQ=255)
                                                                                                                                                                                                               gggCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGACAGTGCTCCACGCTCCTCAAt                                                                                                                                                                                                                                                     >  6:105648/1‑98 (MQ=255)
                                                                                                                                                                                                                  ccATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCt                                                                                                                                                                                                                     >  5:331462/1‑127 (MQ=255)
                                                                                                                                                                                                                    aTAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                                                                                                                                                     <  6:802409/93‑1 (MQ=255)
                                                                                                                                                                                                                    aTAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCg                                                                                                                                                                                               <  5:707834/147‑1 (MQ=255)
                                                                                                                                                                                                                    aTAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCg                                                                                                                                                                                               >  6:707834/1‑147 (MQ=255)
                                                                                                                                                                                                                        ccTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCg                                                                                                                                                                                                                                                                           >  6:213531/1‑67 (MQ=255)
                                                                                                                                                                                                                        ccTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCg                                                                                                                                                                                                                                                                           <  5:213531/67‑1 (MQ=255)
                                                                                                                                                                                                                        ccTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                                                                                                                                                        >  8:1047094/1‑134 (MQ=255)
                                                                                                                                                                                                                        ccTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGAc                                                                                                                                                                                          >  6:1082727/1‑148 (MQ=255)
                                                                                                                                                                                                                          tGGCGGCTGTACTCCGCCGGATACCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTCTACCCGTCAGTGCGCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTgg                                                         <  1:96089/275‑1 (MQ=255)
                                                                                                                                                                                                                                   tACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAAcgcg                                                                                                                                                                               <  5:556397/148‑1 (MQ=255)
                                                                                                                                                                                                                                       ccgccgGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                                                                                                                                                >  6:445278/1‑79 (MQ=255)
                                                                                                                                                                                                                                       ccgccgGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTg                                                                                                                                                                                                                                              >  5:1093350/1‑81 (MQ=255)
                                                                                                                                                                                                                                         gccggtatgctccATTGGGTTGCAATGGAAGCCTGTGACGTGTCTACCCGTGAGTGCTCCACTCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACCCGGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTGGCGgt                                          <  4:198956/266‑1 (MQ=255)
                                                                                                                                                                                                                                           cgGATATTTCTATTGTGATGAAATGGAAGCCTGTGACGTGTATACCCGCCAGTGCTCCACGCTCCTAAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTGGCGgtc                                         <  4:446829/274‑1 (MQ=255)
                                                                                                                                                                                                                                                 tctcTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                                                                                                                                                        <  7:1123676/109‑1 (MQ=255)
                                                                                                                                                                                                                                                 tctcTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                                                                                                                                                        >  7:359356/1‑109 (MQ=255)
                                                                                                                                                                                                                                                 tctcTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                                                                                                                                                        >  8:1123676/1‑109 (MQ=255)
                                                                                                                                                                                                                                                 tctcTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                                                                                                                                                        <  8:359356/109‑1 (MQ=255)
                                                                                                                                                                                                                                                 tctcTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAAcg                                                                                                                                                                                 >  6:97126/1‑132 (MQ=255)
                                                                                                                                                                                                                                                 tctcTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAAcg                                                                                                                                                                                 <  5:97126/132‑1 (MQ=255)
                                                                                                                                                                                                                                                   tctATTGGGATGCAATGGAAGCCTGTGACGTTGATACCCGTAAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTGGCGGTCGTCCCTGGa                                <  1:10529/275‑1 (MQ=255)
                                                                                                                                                                                                                                                    ctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAgg                                                                                                                                                                                                            <  8:51899/102‑1 (MQ=255)
                                                                                                                                                                                                                                                    ctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAgg                                                                                                                                                                                                            >  7:51899/1‑102 (MQ=255)
                                                                                                                                                                                                                                                    ctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                                                                                                                                                           <  6:1125397/119‑1 (MQ=255)
                                                                                                                                                                                                                                                    ctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAAc                                                                                                                                                              >  7:996459/1‑148 (MQ=255)
                                                                                                                                                                                                                                                      aTTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTa                                                                                                                                                                                                                                                 >  8:532417/1‑63 (MQ=255)
                                                                                                                                                                                                                                                      aTTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTa                                                                                                                                                                                                                                                 <  7:532417/63‑1 (MQ=255)
                                                                                                                                                                                                                                                      aTTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGgtggtg                                                                                                                                                                                                       <  6:1072301/105‑1 (MQ=255)
                                                                                                                                                                                                                                                      aTTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                                                                                                                                                           <  8:351538/117‑1 (MQ=255)
                                                                                                                                                                                                                                                      aTTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                                                                                                                                                           >  7:351538/1‑117 (MQ=255)
                                                                                                                                                                                                                                                       ttGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTg                                                                                                                                                           >  7:840131/1‑148 (MQ=255)
                                                                                                                                                                                                                                                            aTGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTGGCGGTCGTCCCTGGAGTGATggg                        >  1:114352/1‑274 (MQ=255)
                                                                                                                                                                                                                                                               cAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACg                                                                                                                                                             <  7:100927/138‑1 (MQ=255)
                                                                                                                                                                                                                                                                aaTGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACAtt                                                                                                                                                  >  7:1083502/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                 aTGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTc                                                                                                                                                 >  7:1026800/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                 aTGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCCAAAGTGGTGTAGGTGGAGGTCTTCTGGCGGTTGTTCCTGGAGTGAGGGGACTTg                   >  2:434652/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                 aTGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTGGCGgtcgtc                                      >  4:151189/1‑255 (MQ=255)
                                                                                                                                                                                                                                                                  tGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGATGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCt                                                                                                                                                <  8:840131/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                   ggAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCTACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTg                                                                                                                                               <  8:156644/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCg                                                                                                                                           <  5:810767/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ccTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGa                                                                                                                                          <  8:1083502/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ccTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGa                                                                                                                                          >  5:475045/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                          tgtgACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTgg                                              >  4:422135/1‑238 (MQ=255)
                                                                                                                                                                                                                                                                           gtgACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAat                                                                                                                                       >  6:913238/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                           ctgACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTgg                                                                                                                                                                                                                   <  8:744480/71‑1 (MQ=255)
                                                                                                                                                                                                                                                                           ctgACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTgg                                                                                                                                                                                                                   >  7:744480/2‑72 (MQ=255)
                                                                                                                                                                                                                                                                           ctgACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                                                                                                                                                           <  6:466466/95‑1 (MQ=255)
                                                                                                                                                                                                                                                                             gACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCt                                                                                                                                                >  5:240086/1‑137 (MQ=255)
                                                                                                                                                                                                                                                                               cGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTAc                                                                                                                                                     <  6:255488/130‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 tgtATACCCTTCAGTGCTCCACGCGCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGAGGATGGGAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTGGAGGTCGTCCCTGGAGTGATGGGAATTGCCGCGTTCTCACCACCg  <  1:37020/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 tgtATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGaa                                                                                                                                                                                                   >  5:872357/1‑82 (MQ=255)
                                                                                                                                                                                                                                                                                 tgtATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACa                                                                                                                                                    >  7:63468/1‑129 (MQ=255)
                                                                                                                                                                                                                                                                                 tgtATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACa                                                                                                                                                    <  8:63468/129‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 tgtATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGGATCGTGTTGGTTTCCCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTGGCGGTCGTCCCTGGCGTGATGGGAATTGCCGGGTTCTCCCCACCg  >  4:122719/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                  
CAGTTATCGCCTTAGAGTTGCATGGCGGCAAACCGCTTTCGCCACTGGTAAATGCTGGCGCTATTGCCACCACCAGCCTGATTAACGCTGAAAATGTTGAACAACGCTGGCAGCGAATTTTACATATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGGCTGTATGGTCGCTCCGGTGACTGGGCGTATCGTGTTGGTTTACCGGGCAAAAGCGGTGTAGGTGGCGGTATTCTGGCGGTCGTCCCTGGAGTGATGGGAATTGCCGCGTTCTCACCACCG  >  NC_000913/511939‑512486

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