Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,739,851 T→C 7.4% F50L (TTC→CTC)  avtA → valine‑pyruvate aminotransferase; transaminase C; alanine‑valine transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,739,8510TC 7.4% 135.1 / 2.8 54F50L (TTC→CTC) avtAvaline‑pyruvate aminotransferase; transaminase C; alanine‑valine transaminase
Reads supporting (aligned to +/- strand):  ref base T (25/25);  new base C (2/2);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

ATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCGAGAAGTTGGGTTGGGATATCGAACCACAGAA  >  NC_000913/3739704‑3739999
                                                                                                                                                   |                                                                                                                                                     
aTGACATCCCCCCTTTTTGGGGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACtt                                                                                                                                                      <  2:239999/149‑1 (MQ=255)
    cATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAg                                                                                                                                                  <  2:760070/149‑1 (MQ=255)
         tcCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCt                                                                                                                                             >  2:264309/1‑149 (MQ=255)
         tcCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATGATGCTCGGCGGCGGGAAGCCGGCGCAGATCCCGGGAATGCAGGACTACTTCCAGACGCt                                                                                                                                             >  1:748613/1‑149 (MQ=255)
          tccTTTTTGGTCACAAATTTCCCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTa                                                                                                                                            <  2:473085/148‑1 (MQ=255)
          cccTTTTTGGCGACAACTTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTa                                                                                                                                            <  2:16759/149‑1 (MQ=255)
                        aaaTTTCCCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCt                                                                                                                              <  1:197541/149‑1 (MQ=255)
                        aaaTTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAACTGCAGGACTACTTCCAGACGCTACTGACCGACATGCt                                                                                                                              >  1:44181/1‑149 (MQ=255)
                        aaaTTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCt                                                                                                                              >  1:525016/1‑149 (MQ=255)
                         aatctccccGCCCCTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTg                                                                                                                             <  1:222935/143‑1 (MQ=255)
                          aTTTACCCGCCACGCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTgg                                                                                                                            <  1:173912/149‑1 (MQ=255)
                                 cGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGGCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGCTCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTgg                                                                                                                     >  2:466686/1‑149 (MQ=255)
                                 cGCCACCCCGGCACTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTgg                                                                                                                     <  1:415277/149‑1 (MQ=255)
                                              ttACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGTCGGCGGTAATCCGGCGCAGATCCCTGAAATGCAGGACTACTTCCAGACGCGACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGAt                                                                                                        <  1:518821/149‑1 (MQ=255)
                                                 cGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCa                                                                                                     <  1:37842/149‑1 (MQ=255)
                                                 cGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCa                                                                                                     <  1:652855/149‑1 (MQ=255)
                                                   cTGTTGATGGAATATCTGAACGACGGTTTACGCACGCCTGGCGCGACTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACt                                                                                                   <  2:92957/149‑1 (MQ=255)
                                                   cTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACt                                                                                                   <  2:309670/149‑1 (MQ=255)
                                                     gTTGATGGAAGATCTGAACGACGGTTTACGCCCGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGCAATGCAGGAATACTTCCAGCCGATACTGCACGACATGCTGCAGCGTGGCAAAGCGACTGATGccgtgt                                                                                                 >  1:362423/1‑144 (MQ=255)
                                                         aTGGAAGAGCTGAACGACTGTTTACGCACGCCGGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑aa                                                                                            <  2:645020/149‑1 (MQ=255)
                                                                   tGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACCTCCAGACGCTACTGACCGACATGTTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGt                                                                                  >  2:450808/1‑149 (MQ=255)
                                                                   tGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACCTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGt                                                                                  >  2:780500/1‑149 (MQ=255)
                                                                              tgatccacgcCTGGCGCGATTCTGCCGGGCGGCGGTAATCCGGCGCAGAGCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCGGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGaa                                                                       <  1:59216/144‑1 (MQ=255)
                                                                                    acgcCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGAAGGACTAGTTCCAGACTCTACTGACCGAGATGGTGGAGTTAGTCAAAGCAACGGAGGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGa                                                                 >  2:768273/1‑149 (MQ=255)
                                                                                        cTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTa                                                             >  1:757938/1‑149 (MQ=255)
                                                                                         tGGGGCGATTATGCCCTGCGGCGGTAATTCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGGCCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTAc                                                            <  2:399589/149‑1 (MQ=255)
                                                                                            cgcgATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTca                                                         <  2:760304/149‑1 (MQ=255)
                                                                                             gcgATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTGCTcac                                                        >  2:47882/1‑149 (MQ=255)
                                                                                              cgATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTcaca                                                       <  2:529660/149‑1 (MQ=255)
                                                                                              cgATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTcaca                                                       <  2:45359/149‑1 (MQ=255)
                                                                                              cgATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTcaca                                                       <  2:27845/149‑1 (MQ=255)
                                                                                                    gggtgggcggcggTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCtt                                                 <  2:602428/144‑1 (MQ=255)
                                                                                                      cTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGc                                               >  2:58657/1‑149 (MQ=255)
                                                                                                       tCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGACAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGcc                                              >  2:287876/1‑149 (MQ=255)
                                                                                                          gcggcggTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGa                                           >  2:383853/1‑149 (MQ=255)
                                                                                                            ggcggTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCCACTCACACTGCTTGCCGGAAt                                         >  1:517054/1‑149 (MQ=255)
                                                                                                                 tAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGGCCACAGGGGAAAACGGAGCTACTCAAACGACTTGCCGGAAtgctgc                                    >  2:686197/1‑149 (MQ=255)
                                                                                                                   aTCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTgcgc                                  >  2:309636/1‑149 (MQ=255)
                                                                                                                      cGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCgag                               >  1:693503/1‑149 (MQ=255)
                                                                                                                      cGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCgag                               >  2:756218/1‑149 (MQ=255)
                                                                                                                      cGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCgag                               >  1:743021/1‑149 (MQ=255)
                                                                                                                      cGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCAACTGACCGACATGCCGGAAAGTGGAAAAGCGACTGCTGCACTGTGT‑AACTACGACGGTCCACAGGGGCAAACCGAGCTACGCACACTGCTTGCCGGAATGCGGcgcggc                               >  1:698951/1‑147 (MQ=255)
                                                                                                                            aGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCGAGAAgttg                         <  2:207208/149‑1 (MQ=255)
                                                                                                                            aGATCACGGGAATGCGGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAGGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCGAGAAgttg                         <  2:561518/149‑1 (MQ=255)
                                                                                                                             gATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCGAGAAgttgg                        >  2:704476/1‑149 (MQ=255)
                                                                                                                             gATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCGAGAAgttgg                        <  1:264309/149‑1 (MQ=255)
                                                                                                                               tCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCGAGAAgttgggt                      >  1:711972/1‑149 (MQ=255)
                                                                                                                                cccGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCGAGAAgttgggtt                     >  1:755504/1‑149 (MQ=255)
                                                                                                                                    gAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGGTGCACTGGGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTAGCCACACTGCTTGCCGGAATGCTGCGCGAGAAGTTGGGGTGGGa                 >  1:280683/1‑149 (MQ=255)
                                                                                                                                        tGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAGGTGGCAAAGCGACTGATGCACTGTGT‑AACGACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCGAGAAGTTGGGTTGGGATATc             >  1:575167/1‑149 (MQ=255)
                                                                                                                                        tGCAGGACTACTGCCAGACGCTACTGACGGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCGACGCACACTGCTTGCCGGAATGCTGCGCGAGAAGTTGGGTTGGGATATc             <  2:716169/149‑1 (MQ=255)
                                                                                                                                           aGGACTACCTCCAGACGCTACTGACCGACATGTTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCGAGAAGTTGGGTTGGGATATCGaa          <  1:450808/149‑1 (MQ=255)
                                                                                                                                           aGGACTACCTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCGAGAAGTTGGGTTGGGATATCGaa          <  1:780500/149‑1 (MQ=255)
                                                                                                                                           aGGACCACTTCCAGACGCTACTGACCGACATGCGGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGCCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCGAGAAGTTGGGTTGGGATATCGaa          <  1:250928/149‑1 (MQ=255)
                                                                                                                                             gACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCGAGAAGTTGGGTTGGGATATCGGAcc        >  2:583834/1‑149 (MQ=255)
                                                                                                                                               ctactTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCGAGAAGTTGGGTTGGGATATCGAACcac      >  2:417391/1‑149 (MQ=255)
                                                                                                                                               ctactTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCGAGAAGTTGGGTTGGGATATCGAACcac      >  1:397486/1‑149 (MQ=255)
                                                                                                                                                  ctTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCGAGAAGTTGGGTTGGGATATCGAACCACAGa   <  1:768273/149‑1 (MQ=255)
                                                                                                                                                  ctTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGAACTGAGTCAAATGAGCTGTTGCACA‑GACAAGCCGGGGCTACTCACACTGCTTGGCGTGATGCTGCGGGAGAAGTTGGGGTGGGATATCTAACCACAGa   >  1:386251/1‑149 (MQ=255)
                                                                                                                                                   ttCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGGCCTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCGAGAAGTTGGGTTGGGCTATCGAACCACAGaa  >  1:572441/1‑149 (MQ=255)
                                                                                                                                                   ttCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGGAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGGGAGAAGTTGGGTTGGGATATCGAACCACAGaa  >  1:16967/1‑149 (MQ=255)
                                                                                                                                                   ttCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCGAGAAGTTGGGTTGGGATATCGAACCACAGaa  >  2:720212/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGGACTACTTCCAGACGCTACTGACCGACATGCTGGAAAGTGGCAAAGCGACTGATGCACTGTGT‑AACTACGACGGTCCACAGGGGAAAACGGAGCTACTCACACTGCTTGCCGGAATGCTGCGCGAGAAGTTGGGTTGGGATATCGAACCACAGAA  >  NC_000913/3739704‑3739999

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑