Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 512,212 T→G 8.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 8.3% 453.5 / 29.8 132V191G (GTG→GGG) glsAglutaminase 1
Reads supporting (aligned to +/- strand):  ref base T (63/58);  new base G (6/5);  total (69/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.40e-01

Only 100 of 140 total aligned reads displayed.

CTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATG  >  NC_000913/512097‑512357
                                                                                                                   |                                                                                                                                                 
        cGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACgtg                                                                                                                                                  <  3:10590/109‑1 (MQ=255)
        cGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACgtg                                                                                                                                                  >  4:10590/1‑109 (MQ=255)
          aaGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCa                                                                                                                      >  2:597479/1‑135 (MQ=255)
          aaGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCa                                                                                                                      <  1:597479/135‑1 (MQ=255)
              cAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCa                                                                                                     >  3:172297/1‑148 (MQ=255)
                 ccAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                               <  2:534026/135‑1 (MQ=255)
                 ccAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                               >  1:534026/1‑135 (MQ=255)
                     tCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTg                                                                                              >  1:12474/1‑148 (MQ=255)
                         aacaaacaACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACTCTCCTCAATACt                                                                                                                 >  2:839960/1‑125 (MQ=255)
                         aacaaacaACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACTCTCCTCAATACt                                                                                                                 <  1:839960/125‑1 (MQ=255)
                          acaaacaaCCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGcgcg                                                                                         >  2:810767/1‑148 (MQ=255)
                              acaacCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                    <  3:50261/117‑1 (MQ=255)
                              acaacCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                    >  4:50261/1‑117 (MQ=255)
                              acaacCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGc                                                                                     >  4:789994/1‑148 (MQ=255)
                               caaccaacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCt                                                                                    <  3:789994/148‑1 (MQ=255)
                                 accaacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGa                                                                                                            >  2:759832/1‑122 (MQ=255)
                                 accaacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTgg                                                                                  >  4:115134/1‑148 (MQ=255)
                                   caacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTggcg                                                                                >  3:513644/1‑148 (MQ=255)
                                   caacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTggcg                                                                                >  3:66394/1‑148 (MQ=255)
                                    aacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACgtg                                                                                                                                                  >  2:1016263/1‑81 (MQ=255)
                                    aacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACgtg                                                                                                                                                  <  1:1016263/81‑1 (MQ=255)
                                    aacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGctc                                                                                                                          <  1:448510/105‑1 (MQ=255)
                                    aacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGctc                                                                                                                          >  2:448510/1‑105 (MQ=255)
                                    aacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTggcgg                                                                               <  1:1064850/148‑1 (MQ=255)
                                           aTAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGgcg                                                                                           >  3:430752/1‑129 (MQ=255)
                                           aTAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGgcg                                                                                           <  4:430752/129‑1 (MQ=255)
                                           aTAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtgg                                                                        <  2:12474/148‑1 (MQ=255)
                                            tAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                               <  4:1022435/108‑1 (MQ=255)
                                            tAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                               >  3:1022435/1‑108 (MQ=255)
                                            tAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                       <  4:172297/148‑1 (MQ=255)
                                            tAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                       >  3:156644/1‑148 (MQ=255)
                                            tAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                       >  2:556397/1‑148 (MQ=255)
                                             aaCCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGgtggtg                                                                      <  4:513644/148‑1 (MQ=255)
                                             aaCCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGgtggtg                                                                      <  4:66394/148‑1 (MQ=255)
                                              aCCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                    >  2:1076467/1‑101 (MQ=255)
                                              aCCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                    <  1:1076467/101‑1 (MQ=255)
                                                 gggCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTTTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                    <  1:105648/98‑1 (MQ=255)
                                                 gggCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGACAGTGCTCCACGCTCCTCAAt                                                                                                                    >  2:105648/1‑98 (MQ=255)
                                                    ccATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCt                                                                                    <  2:331462/127‑1 (MQ=255)
                                                     cATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCg                                                              <  3:115134/148‑1 (MQ=255)
                                                      aTAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                    >  1:802409/1‑93 (MQ=255)
                                                      aTAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCg                                                              >  2:707834/1‑147 (MQ=255)
                                                      aTAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCg                                                              <  1:707834/147‑1 (MQ=255)
                                                        aGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATc                                                                <  3:174159/143‑1 (MQ=255)
                                                        aGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATc                                                                >  4:174159/1‑143 (MQ=255)
                                                        aGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTg                                                           >  3:597514/1‑148 (MQ=255)
                                                          ccTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCg                                                                                                                                          >  2:213531/1‑67 (MQ=255)
                                                                     tACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAAcgcg                                              <  1:556397/148‑1 (MQ=255)
                                                                         ccgccgGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                               >  2:445278/1‑79 (MQ=255)
                                                                         ccgccgGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                               <  1:445278/79‑1 (MQ=255)
                                                                         ccgccgGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTg                                                                                                             <  2:1093350/81‑1 (MQ=255)
                                                                         ccgccgGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACg                                                        >  1:926337/1‑134 (MQ=255)
                                                                         ccgccgGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACg                                                        <  2:926337/134‑1 (MQ=255)
                                                                           gccgGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGttcttc                                        <  1:1082727/148‑1 (MQ=255)
                                                                              ggATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAgg                                     <  4:597514/148‑1 (MQ=255)
                                                                               gATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGt                                                             <  4:42920/123‑1 (MQ=255)
                                                                               gATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGt                                                             >  3:42920/1‑123 (MQ=255)
                                                                                   tctcTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                       >  4:1123676/1‑109 (MQ=255)
                                                                                   tctcTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                       <  3:1123676/109‑1 (MQ=255)
                                                                                   tctcTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAAcg                                                <  1:97126/132‑1 (MQ=255)
                                                                                   tctcTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAAcg                                                >  2:97126/1‑132 (MQ=255)
                                                                                      ctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAgg                                                                           <  4:51899/102‑1 (MQ=255)
                                                                                      ctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAgg                                                                           >  3:51899/1‑102 (MQ=255)
                                                                                      ctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                          <  2:1125397/119‑1 (MQ=255)
                                                                                      ctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                          >  3:829796/1‑119 (MQ=255)
                                                                                      ctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                          <  4:829796/119‑1 (MQ=255)
                                                                                      ctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                          >  1:1125397/1‑119 (MQ=255)
                                                                                      ctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAAc                             >  3:996459/1‑148 (MQ=255)
                                                                                        aTTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTa                                                                                                                >  4:532417/1‑63 (MQ=255)
                                                                                        aTTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTa                                                                                                                <  3:532417/63‑1 (MQ=255)
                                                                                        aTTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGgtggtg                                                                      <  2:1072301/105‑1 (MQ=255)
                                                                                        aTTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                          <  4:351538/117‑1 (MQ=255)
                                                                                        aTTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                          >  3:351538/1‑117 (MQ=255)
                                                                                         ttGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTg                          >  3:840131/1‑148 (MQ=255)
                                                                                                gCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACa                   >  2:487600/1‑148 (MQ=255)
                                                                                                 cAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGt                                                             >  1:141536/1‑105 (MQ=255)
                                                                                                 cAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACg                            <  3:100927/138‑1 (MQ=255)
                                                                                                 cAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACg                            >  4:100927/1‑138 (MQ=255)
                                                                                                  aaTGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACAtt                 >  3:1083502/1‑148 (MQ=255)
                                                                                                   aTGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTc                >  3:1026800/1‑148 (MQ=255)
                                                                                                    tGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGATGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCt               <  4:840131/148‑1 (MQ=255)
                                                                                                     ggAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCTACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTg              <  4:156644/148‑1 (MQ=255)
                                                                                                         gCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCg          <  1:810767/148‑1 (MQ=255)
                                                                                                          ccTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGa         <  4:1083502/148‑1 (MQ=255)
                                                                                                            tgtgACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGgtgt                                                                     >  3:241140/1‑86 (MQ=255)
                                                                                                            tgtgACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGgtgt                                                                     <  4:241140/86‑1 (MQ=255)
                                                                                                             gtgACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAat      >  1:750643/1‑148 (MQ=255)
                                                                                                             ctgACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTgg                                                                                  >  3:744480/2‑72 (MQ=255)
                                                                                                             ctgACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTgg                                                                                  <  4:744480/71‑1 (MQ=255)
                                                                                                             ctgACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                          >  1:466466/2‑96 (MQ=255)
                                                                                                             ctgACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCa                                       >  4:354730/2‑115 (MQ=255)
                                                                                                             ctgACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCa                                       <  3:354730/114‑1 (MQ=255)
                                                                                                              tgACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAatg     >  2:1045671/1‑148 (MQ=255)
                                                                                                               gACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCt               <  2:240086/137‑1 (MQ=255)
                                                                                                               gACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCt               <  2:587904/137‑1 (MQ=255)
                                                                                                               gACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCt               >  1:587904/1‑137 (MQ=255)
                                                                                                               gACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCt               >  1:240086/1‑137 (MQ=255)
                                                                                                                 cGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTAc                    >  1:255488/1‑130 (MQ=255)
                                                                                                                 cGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAatgatg  <  1:487600/148‑1 (MQ=255)
                                                                                                                   tgtATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGaa                                                                  <  2:872357/82‑1 (MQ=255)
                                                                                                                   |                                                                                                                                                 
CTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATG  >  NC_000913/512097‑512357

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