Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,952,942 A→G T154A (ACC→GCC)  ilvE → branched‑chain amino acid aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,952,9420AG98.4% 201.0 / NA 61T154A (ACC→GCC) ilvEbranched‑chain amino acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (31/29);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAGGCGAAAACCTGTTT  >  NC_000913/3952798‑3953085
                                                                                                                                                |                                                                                                                                               
gTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGCTCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCAt                                                                                                                                             >  2:108878/1‑149 (MQ=255)
      aTGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGAc                                                                                                                                       >  1:717746/1‑149 (MQ=255)
      aTGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGAACCAAACGCCATCCCGAc                                                                                                                                       >  1:249162/1‑149 (MQ=255)
          gAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGAcggcg                                                                                                                                   <  2:657075/149‑1 (MQ=255)
          gAGGAACCCCGCCTGCGGGATCCTCAACCGACGTGATTATCGCTGCCTTCCCGTGGGCACCGTATGTGGGCGCAGCAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGAcggcg                                                                                                                                   <  2:263591/149‑1 (MQ=255)
           aGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGAcggcgg                                                                                                                                  <  1:627921/149‑1 (MQ=255)
             tAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATGTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCa                                                                                                                                >  1:628044/1‑149 (MQ=255)
             tAAACCCGCCAGCGGGATACTCAACCGACGTGATTAGCGCGGCTTTCCCGTGGGGGGCGTAGCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCa                                                                                                                                >  2:695481/1‑149 (MQ=255)
                aCCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGGCGGCGGCaaag                                                                                                                             >  2:750108/1‑148 (MQ=255)
                 cccGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTATCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCaaacg                                                                                                                            >  2:646558/1‑147 (MQ=255)
                  ccGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGc                                                                                                                           >  2:476514/1‑149 (MQ=255)
                   cGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGcc                                                                                                                          >  1:166769/1‑149 (MQ=255)
                      cAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCggt                                                                                                                       <  1:740591/149‑1 (MQ=255)
                       aGCGGGATACTCGACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCggtg                                                                                                                      <  2:836176/149‑1 (MQ=255)
                       aGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCGGCACCAAACGCCATCCCGACGGCGGCAAAAGCCggtg                                                                                                                      <  2:190655/149‑1 (MQ=255)
                                 tCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACct                                                                                                            <  1:355976/149‑1 (MQ=255)
                                    aCCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctc                                                                                                         >  1:872823/1‑149 (MQ=255)
                                         cGTGAGTAGCGCGGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCTCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAGCCGCGCAGCATCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTccc                                                                                                    <  2:120082/149‑1 (MQ=255)
                                           tGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctg                                                                                                  >  2:895231/1‑149 (MQ=255)
                                            gATTATCCCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgc                                                                                                 <  1:786948/149‑1 (MQ=255)
                                             aTTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgct                                                                                                >  1:860456/1‑149 (MQ=255)
                                               tATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtgg                                                                                              <  2:655507/149‑1 (MQ=255)
                                                aTCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCTGGTAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggt                                                                                             <  1:695481/149‑1 (MQ=255)
                                                  cGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTCACTACCTCTCTTCCCTGCtggtgg                                                                                           >  1:303614/1‑149 (MQ=255)
                                                      gctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTTACTACCTCTCTTCCCTGCTGGTGGGTAg                                                                                       <  1:523492/149‑1 (MQ=255)
                                                          tCCCGTGGGGAGCGTATCGGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGaa                                                                                   <  1:936439/149‑1 (MQ=255)
                                                              gTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGACgcgc                                                                               >  1:784955/1‑149 (MQ=255)
                                                              gTGGGGAGCGTATCTGGGCGCAGAAGCACTGGAGCAGGGGATCGATGGGATGGTTTCCTCCTGTCACCTCGCAGCACTAAACACGAACCCGACGGCGGCAAGAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTTGGTAGCGAAgcgc                                                                               >  2:924/1‑149 (MQ=255)
                                                                              ggCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCa                                                               <  1:652285/149‑1 (MQ=255)
                                                                                 gcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCCGGa                                                            >  2:429155/1‑149 (MQ=255)
                                                                                          cTGGAGCAGGGGATCGCTGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTCACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAATGTacccc                                                   <  2:628044/149‑5 (MQ=255)
                                                                                                   gggATCGATGCGATGGTTTCCTCCGGGAGCCGGGCAGCACCAAACGCCAGCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGAtgt                                          <  2:772155/149‑1 (MQ=255)
                                                                                                   gggATCGAGGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGAtgt                                          <  1:13586/149‑1 (MQ=255)
                                                                                                    ggATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGAtgtg                                         >  2:98310/1‑149 (MQ=255)
                                                                                                     gATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGCTGTGa                                        >  1:552770/1‑149 (MQ=255)
                                                                                                      aTCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAAGAAGGTATCGCGCTGGATGTGaa                                       <  1:519515/149‑1 (MQ=255)
                                                                                                         gctgcgacgGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAATCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACgg                                    <  2:427869/147‑1 (MQ=255)
                                                                                                         gatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACgg                                    >  1:588179/1‑149 (MQ=255)
                                                                                                         gatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACgg                                    >  1:479450/1‑149 (MQ=255)
                                                                                                         gatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCCCGCTGGATGTGGACgg                                    >  1:6220/1‑149 (MQ=255)
                                                                                                          acgcggtgGTTGCCTCCTGGAACCGCGCCGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGt                                   <  2:846866/147‑1 (MQ=255)
                                                                                                            gcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTta                                 >  2:573796/1‑149 (MQ=255)
                                                                                                              gatgGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTtatg                               >  2:4366/1‑148 (MQ=255)
                                                                                                               atgGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTatat                              >  1:808269/1‑149 (MQ=255)
                                                                                                                   tttCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAACAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATctct                          <  1:128024/149‑1 (MQ=255)
                                                                                                                   tttCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACgg                                    <  2:575112/139‑1 (MQ=255)
                                                                                                                   tttCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACgg                                    >  1:575112/1‑139 (MQ=255)
                                                                                                                   tttCCGCCCGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATctct                          <  1:7431/149‑1 (MQ=255)
                                                                                                                        tcctGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGGAACTACCTCTCTTGCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGGACGGGTATATCTGTGGAgg                     >  1:198738/1‑149 (MQ=255)
                                                                                                                          cgggAACCGCGCAGCAACAAACGCCAACCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGgcg                   <  1:87218/147‑1 (MQ=255)
                                                                                                                             gAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAg                >  1:202364/1‑149 (MQ=255)
                                                                                                                                 cgcgCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAGGCGa            >  2:445261/1‑149 (MQ=255)
                                                                                                                                    gcagcaCCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAGGCGaaaa         >  1:104416/1‑149 (MQ=255)
                                                                                                                                      agcaCCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTACCGAAGCGCGCCGCCACGGTTATCAGGACGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCCGGCGAAAAcc       >  2:38379/1‑149 (MQ=255)
                                                                                                                                         cccAAACGCCATCCCGACGGCGGCAAAAGCCGGGGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAGGCGAAAACCTGt    <  2:930226/148‑1 (MQ=255)
                                                                                                                                         aCCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAGGCGAAAACCTGt    <  2:588179/149‑1 (MQ=255)
                                                                                                                                          ccAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAGGCGAAAACCTGtt   >  2:168036/1‑149 (MQ=255)
                                                                                                                                           cAAACGCCATGCCGGCGGCGGCAAAAGCCGGCGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTTTCAGGAAGGTATCGCGCTGGAGGTGAACGGTTATATCTCTGAAGGCGCAGGCGAAAACCTGttt  <  1:175103/149‑1 (MQ=255)
                                                                                                                                           cAAACGCCATCCCGACGGTGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTATGAAGGCGCAGGCGAAAACCTGttc  >  1:693264/1‑148 (MQ=255)
                                                                                                                                           cAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAGGCGAAAACCTGttt  <  1:750108/149‑1 (MQ=255)
                                                                                                                                           cAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAGGCGAAAACCTGttt  <  2:808269/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAGGCGAAAACCTGTTT  >  NC_000913/3952798‑3953085

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: ‑