Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,829,334 G→T 9.1% P134Q (CCG→CAG)  astD ← succinylglutamic semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,829,3340GT 9.1% 179.4 / 8.6 55P134Q (CCG→CAG) astDsuccinylglutamic semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (30/20);  new base T (3/2);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.91e-02

GCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAGCGCCTGACCCGTTTCACGCCCGCCCTGCACCAGGTTCAGCACGCCTGGCGGCAAGCCAGCCTGCTGCCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTG  >  NC_000913/1829072‑1829601
                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                           
gcTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAGCGCCTGACCCGTTTCACGCCCGCCCTGCACCAGGTTCAGCACGCCTGGCGGCAAGCCAGCCTGCTGCCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCa                                                                                                                                                                                                                                                                 >  2:242404/1‑275 (MQ=255)
   aCCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAGCGCCTGACCCGTTTCACGCCCGCCCTGCACCAGGTTCAGCACGCCTGGCGGCAAGCCAGCCTGCTGCCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTTCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACCGCCCAAACACCCCCACCCCCCCGTCCGGGCGAGATCGCAgt                                                                                                                                                                                                                                                               >  4:328648/1‑273 (MQ=255)
               cAAACCGTCGAGATCCTCCAGCGCACTCAGCGCCTGACCCGTTTCACGCCCGCCCTGCACCAGGTTCAGCACGCCTGGCGGCAAGCCAGCCTGCTGCCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCCcgccgcgccg                                                                                                                                                                                                                                                   >  2:333772/1‑274 (MQ=255)
                      tCGAGATCCTCCAGCGCACTCAGCGCCTGACCCGTTTCACGCCCGCCCTGCACCAGGTTCAGCACGCCTGGCGGCAAGCCAGCCTGCTGCCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCa                                                                                                                                                                                                                                            <  2:328652/274‑1 (MQ=255)
                                                   aCCCGTTTCACGCCCGCCCTGCACCAGGTTCAGCACGCCTGGCGGCAAGCCAGCCTGCTGCCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAACTGACCAGGGAAATTATACGGCCCAAACACCACCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCCCCGGTACGAACgg                                                                                                                                                                                                              >  1:59189/1‑274 (MQ=255)
                                                                           cAGGTTCGTCACGCCTGGCGGCAAGCCAGCCTGCTGCCAAAAACGCATTACCGCTTCGCCACTCCACGGGGTCAGTTCTCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGa                                                                                                                                                                                      <  2:386557/275‑1 (MQ=255)
                                                                                  aGCACGCCTGGCGGCAAGCCAGCCTGCTGCCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGGGATACGCCTTAATTGATCTAGCGATTTTAt                                                                                                                                                                                >  1:8202/1‑274 (MQ=255)
                                                                                                 aaGCCAGCCTGCTGCCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCa                                                                                                                                                                 >  2:263292/1‑274 (MQ=255)
                                                                                                                  tAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCg                                                                                                                                                                                                                                                                             >  7:260833/1‑149 (MQ=255)
                                                                                                                    aaCGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCggg                                                                                                                                                                                                                                                                           >  7:211479/1‑149 (MQ=255)
                                                                                                                           tACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCCGCACGGCGTGCTGGCGgtgtt                                                                                                                                                                                                                                                                    <  8:334252/149‑2 (MQ=255)
                                                                                                                                  tCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCt                                                                                                                                                                                                                                                             <  7:394897/149‑1 (MQ=255)
                                                                                                                                    gCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcg                                                                                                                                                                                                                                                           >  7:200707/1‑149 (MQ=255)
                                                                                                                                     ccACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGAt                                                                                                                                                                                     >  2:411205/1‑218 (MQ=255)
                                                                                                                                      ccatCCACGGTGTTAGTTCGCGGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTtct                                                                                                                           <  2:386861/272‑1 (MQ=255)
                                                                                                                                         tCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTa                                                                                                                                                                                 <  4:411200/218‑1 (MQ=255)
                                                                                                                                           cACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgccgcgcc                                                                                                                                                                                                                                                    >  7:223479/1‑149 (MQ=255)
                                                                                                                                                          cTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGtt                                                                                                         >  2:44174/1‑273 (MQ=255)
                                                                                                                                                            gggTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTAc                                                                                                                                                                                                                                   <  8:195901/149‑1 (MQ=255)
                                                                                                                                                              gTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAATGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAgcg                                                                                                                                           <  2:112150/235‑1 (MQ=255)
                                                                                                                                                               tttAAAGGTAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACgct                                                                                                                                                                                                                                >  6:261167/1‑149 (MQ=255)
                                                                                                                                                                    aGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCg                                                                                                                                                                                                                           >  6:168510/1‑149 (MQ=255)
                                                                                                                                                                     gagaaTGGTGTGACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTCAt                                                                                            <  1:92343/272‑1 (MQ=255)
                                                                                                                                                                      aTAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCCGCACGGCGTGCTGGCgg                                                                                                                                                                                                                                                                        <  7:79120/102‑2 (MQ=255)
                                                                                                                                                                      aTAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCCGCACGGCGTGCTGGCgg                                                                                                                                                                                                                                                                        >  8:79120/1‑101 (MQ=255)
                                                                                                                                                                       tAATTGTTTTACCTGCCAGCCATGCCCGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGGGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTg                                                                                          <  2:382516/275‑1 (MQ=255)
                                                                                                                                                                       tAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCCCTGGCAATAATCGCGGTTAATTCGGCTTTATTg                                                                                          >  2:233210/1‑275 (MQ=255)
                                                                                                                                                                          tGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgc                                                                                                                                                                                                                                                          >  8:17237/1‑112 (MQ=255)
                                                                                                                                                                          tGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgc                                                                                                                                                                                                                                                          <  7:17237/112‑1 (MQ=255)
                                                                                                                                                                            gtgtTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACg                                                                                                                                                                                                                   >  6:374619/1‑149 (MQ=255)
                                                                                                                                                                                 aCCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGt                                                                                                                                                                                                              <  5:374619/149‑1 (MQ=255)
                                                                                                                                                                                         gCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGcc                                                                                                                                                                                                      <  7:43772/149‑1 (MQ=255)
                                                                                                                                                                                          cAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTAt                                                                                                                                                                                <  2:74170/170‑1 (MQ=255)
                                                                                                                                                                                             tGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCCATTTTATTGATCATCGCCGTCCCTTCGGGTGCCCCTTCCCCACGCGGCTTACCCGTTTCCCCGGCAAAAATCCCTGTTTATTTNNNCTTTTTTCTTTCCCTCAGTTCCGCAAAg                                                                     >  3:145858/1‑274 (MQ=255)
                                                                                                                                                                                                 ggAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGa                                                                                                                                                                                              >  8:251419/1‑149 (MQ=255)
                                                                                                                                                                                                   aaCGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGata                                                                                                                                                                                            >  8:170340/1‑149 (MQ=255)
                                                                                                                                                                                                    actatatGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGAc                                                             <  2:428872/272‑1 (MQ=255)
                                                                                                                                                                                                         atGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGa                                                                                                                                                                                      <  7:170340/149‑1 (MQ=255)
                                                                                                                                                                                                          tGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGAt                                                                                                                                                                                     <  8:223479/149‑1 (MQ=255)
                                                                                                                                                                                                            tCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGAttt                                                                                                                                                                                   >  6:163109/1‑149 (MQ=255)
                                                                                                                                                                                                                  cGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGACAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGtt                                               <  2:362153/275‑1 (MQ=255)
                                                                                                                                                                                                                   ggCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGat                                                                                                                                                                            >  6:80707/1‑149 (MQ=255)
                                                                                                                                                                                                                        aTGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTAtt                                                                                           >  1:50154/1‑225 (MQ=255)
                                                                                                                                                                                                                          gACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGcc                                                                                                                                                                     >  8:49800/1‑149 (MQ=255)
                                                                                                                                                                                                                               gggAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCAc                                                                                                                                                                >  6:173966/1‑149 (MQ=255)
                                                                                                                                                                                                                                     ttATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGt                                                                                                                                                          >  5:402350/1‑149 (MQ=255)
                                                                                                                                                                                                                                        taCGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGc                                                                                                                                                       >  5:435124/1‑149 (MQ=255)
                                                                                                                                                                                                                                            gCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCt                                                                                                                                                   >  5:244471/1‑149 (MQ=255)
                                                                                                                                                                                                                                             cccAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGt                                                                                                                                                                                                              >  4:97221/1‑89 (MQ=255)
                                                                                                                                                                                                                                             cccAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCCTGACGTTCAGCAAATGAAAGCCGCGACCAGcgcg                    >  2:101765/1‑275 (MQ=255)
                                                                                                                                                                                                                                                 aaCACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATa                                                                                                                                                                                                          <  2:97222/89‑1 (MQ=255)
                                                                                                                                                                                                                                                   cacCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAgc                                                                                                                                            <  5:163109/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    atcGCCCGCACGGCGTGCTGGCGGTGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAgcg                                                                                                                                           >  5:282022/3‑149 (MQ=255)
                                                                                                                                                                                                                                                    atcGCCCGCACGGCGTGCTGGCGGTGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAgcg                                                                                                                                           >  8:195692/3‑149 (MQ=255)
                                                                                                                                                                                                                                                       gcccgcccGCCGTGCGGGCGTTTTCCCAGGCTCCCCGCGCCGTCCTGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGTTATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAAcgccgc          <  1:145858/275‑1 (MQ=255)
                                                                                                                                                                                                                                                             aCGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTAcc                                                                                                                                  >  6:148725/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                ccGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCCGTTCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCCCGGAAACGCCGCACGGGATg  >  2:282440/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                           
GCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAGCGCCTGACCCGTTTCACGCCCGCCCTGCACCAGGTTCAGCACGCCTGGCGGCAAGCCAGCCTGCTGCCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTG  >  NC_000913/1829072‑1829601

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