Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,829,322 A→C 8.8% L138R (CTG→CGG)  astD ← succinylglutamic semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,829,3220AC 8.8% 177.5 / 10.8 57L138R (CTG→CGG) astDsuccinylglutamic semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (32/20);  new base C (3/2);  total (35/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 = 2.59e-01

GCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAGCGCCTGACCCGTTTCACGCCCGCCCTGCACCAGGTTCAGCACGCCTGGCGGCAAGCCAGCCTGCTGCCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGC  >  NC_000913/1829072‑1829593
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                               
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)
                                                                                                      agccTGCTGCCATAAACGCATTACCGCTGCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAAGGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCa                                                                                                                                                                                                                                                                                 <  8:99182/149‑1 (MQ=255)
                                                                                                                 aTAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGc                                                                                                                                                                                                                                                                      >  8:394897/1‑149 (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)
                                                                                                                                                     gTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGc                                                                                                                                                                                                                                                                           <  6:32769/108‑1 (MQ=255)
                                                                                                                                                     gTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGc                                                                                                                                                                                                                                                                           >  5:32769/1‑108 (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                                                                                                                                                                                                                                                                >  8:79120/1‑101 (MQ=255)
                                                                                                                                                                      aTAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCCGCACGGCGTGCTGGCgg                                                                                                                                                                                                                                                                <  7:79120/102‑2 (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)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                               
GCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAGCGCCTGACCCGTTTCACGCCCGCCCTGCACCAGGTTCAGCACGCCTGGCGGCAAGCCAGCCTGCTGCCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGC  >  NC_000913/1829072‑1829593

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