Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 512,219 C→A 9.7% 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 9.7% 495.2 / 33.1 143T193T (ACC→ACAglsAglutaminase 1
Reads supporting (aligned to +/- strand):  ref base C (67/62);  new base A (7/7);  total (74/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-01

Only 100 of 149 total aligned reads displayed.

CTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGG  >  NC_000913/512097‑512366
                                                                                                                          |                                                                                                                                                   
          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)
                  cAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                        <  4:18409/134‑1 (MQ=255)
                     tCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTg                                                                                                       >  1:12474/1‑148 (MQ=255)
                      cGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                        <  4:464159/130‑1 (MQ=255)
                         aacaaacaACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACTCTCCTCAATACt                                                                                                                          <  1:839960/125‑1 (MQ=255)
                          acaaacaaCCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGcgcg                                                                                                  >  2:810767/1‑148 (MQ=255)
                             aacaacCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                             <  2:60995/118‑1 (MQ=255)
                              acaacCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                             >  4:164253/1‑117 (MQ=255)
                              acaacCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                             <  3:50261/117‑1 (MQ=255)
                               caaccaacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCt                                                                                             <  3:789994/148‑1 (MQ=255)
                                 accaacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTTTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTgg                                                                                           >  4:977541/1‑148 (MQ=255)
                                 accaacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGa                                                                                                                     <  1:759832/122‑1 (MQ=255)
                                   caacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTggcg                                                                                         >  3:513644/1‑148 (MQ=255)
                                   caacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTggcg                                                                                         >  3:66394/1‑148 (MQ=255)
                                   caacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTggcg                                                                                         >  3:1113955/1‑148 (MQ=255)
                                    aacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGctc                                                                                                                                   >  2:448510/1‑105 (MQ=255)
                                    aacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTggcgg                                                                                        <  1:1064850/148‑1 (MQ=255)
                                        tCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTg                                                                                            <  1:538465/140‑1 (MQ=255)
                                           aTAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGa                                                                                                                     >  3:795510/1‑112 (MQ=255)
                                           aTAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGa                                                                                                                     <  4:795510/112‑1 (MQ=255)
                                           aTAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGgcg                                                                                                    >  3:430752/1‑129 (MQ=255)
                                           aTAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGgcg                                                                                                    <  4:430752/129‑1 (MQ=255)
                                            tAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                        <  4:1022435/108‑1 (MQ=255)
                                            tAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                        >  3:1022435/1‑108 (MQ=255)
                                            tAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                                >  3:156644/1‑148 (MQ=255)
                                            tAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                                >  2:556397/1‑148 (MQ=255)
                                             aaCCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGgtggtg                                                                               <  4:66394/148‑1 (MQ=255)
                                             aaCCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGgtggtg                                                                               <  4:513644/148‑1 (MQ=255)
                                             aaCCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGgtggtg                                                                               <  4:1113955/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                                                                                             >  1:331462/1‑127 (MQ=255)
                                                    ccATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCt                                                                                             <  2:331462/127‑1 (MQ=255)
                                                     cATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCg                                                                       <  3:115134/148‑1 (MQ=255)
                                                      aTAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                             >  1:802409/1‑93 (MQ=255)
                                                      aTAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                             <  2:802409/93‑1 (MQ=255)
                                                      aTAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCg                                                                       <  1:707834/147‑1 (MQ=255)
                                                      aTAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCg                                                                       >  2:707834/1‑147 (MQ=255)
                                                        aGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATc                                                                         <  3:174159/143‑1 (MQ=255)
                                                        aGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATc                                                                         >  4:174159/1‑143 (MQ=255)
                                                          ccTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCg                                                                                                                                                   >  2:213531/1‑67 (MQ=255)
                                                          ccTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCg                                                                                                                                                   <  1:213531/67‑1 (MQ=255)
                                                          ccTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                                >  4:1047094/1‑134 (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                                                                      >  3:42920/1‑123 (MQ=255)
                                                                                   tctcTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                                <  4:359356/109‑1 (MQ=255)
                                                                                   tctcTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                                >  4:1123676/1‑109 (MQ=255)
                                                                                   tctcTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAAcg                                                         >  2:97126/1‑132 (MQ=255)
                                                                                   tctcTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAAcg                                                         <  1:97126/132‑1 (MQ=255)
                                                                                      ctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAgg                                                                                    <  4:51899/102‑1 (MQ=255)
                                                                                      ctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                                   >  1:1125397/1‑119 (MQ=255)
                                                                                      ctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                                   >  3:829796/1‑119 (MQ=255)
                                                                                      ctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                                   <  2:1125397/119‑1 (MQ=255)
                                                                                      ctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAAc                                      >  3:996459/1‑148 (MQ=255)
                                                                                        aTTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGctcc                                                                                                                                  <  3:519651/54‑1 (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)
                                                                                         ttGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTg                                   >  3:840131/1‑148 (MQ=255)
                                                                                                gCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACa                            >  2:487600/1‑148 (MQ=255)
                                                                                                 cAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGt                                                                      <  2:141536/105‑1 (MQ=255)
                                                                                                 cAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACg                                     <  3:100927/138‑1 (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)
                                                                                                          ccTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGa                  >  1:475045/1‑148 (MQ=255)
                                                                                                            tgtgACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGgtgt                                                                              <  4:241140/86‑1 (MQ=255)
                                                                                                            tgtgACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGgtgt                                                                              >  3:241140/1‑86 (MQ=255)
                                                                                                             ctgACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTgg                                                                                           <  4:744480/71‑1 (MQ=255)
                                                                                                             ctgACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                                   <  2:466466/95‑1 (MQ=255)
                                                                                                             ctgACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGa                                                                   >  1:466466/2‑96 (MQ=255)
                                                                                                             ctgACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCa                                                <  3:354730/114‑1 (MQ=255)
                                                                                                              tgACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAatg              >  2:1045671/1‑148 (MQ=255)
                                                                                                               gACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCt                        >  1:240086/1‑137 (MQ=255)
                                                                                                               gACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCt                        >  1:587904/1‑137 (MQ=255)
                                                                                                               gACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCt                        <  2:587904/137‑1 (MQ=255)
                                                                                                                 cGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTAc                             <  2:255488/130‑1 (MQ=255)
                                                                                                                  gggtATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAatgatga          <  1:1045671/146‑1 (MQ=255)
                                                                                                                   tgtATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGaa                                                                           >  1:872357/1‑82 (MQ=255)
                                                                                                                   tgtATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGaa                                                                           <  2:872357/82‑1 (MQ=255)
                                                                                                                   tgtATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACa                            >  3:63468/1‑129 (MQ=255)
                                                                                                                     tataCACGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                        <  3:357442/35‑1 (MQ=255)
                                                                                                                     tataCACGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                        >  4:357442/1‑35 (MQ=255)
                                                                                                                       taCCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGaa     >  1:658920/1‑148 (MQ=255)
                                                                                                                        aCCCGTCAGTGCTCCACGCTGCTCAGTACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGAc                                                                  >  3:200373/1‑86 (MQ=255)
                                                                                                                         cccGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGgtggtg                                                                               >  2:885075/1‑72 (MQ=255)
                                                                                                                         cccGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGgtggtg                                                                               >  4:811587/1‑72 (MQ=255)
                                                                                                                         cccGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAgg   >  3:998608/1‑148 (MQ=255)
                                                                                                                          |                                                                                                                                                   
CTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGGAAGGG  >  NC_000913/512097‑512366

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: