Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 23,255 G→A 5.4% D289N (GAT→AAT)  ileS → isoleucyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091323,2550GA 5.4% 243.9 / 3.2 75D289N (GAT→AAT) ileSisoleucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (38/32);  new base A (2/2);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

GGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCA  >  NC_000913/23109‑23403
                                                                                                                                                   |                                                                                                                                                    
ggACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGAt                                                                                                                                                    >  3:217730/1‑149 (MQ=255)
        acGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCAATTACACCAt                                                                                                                                            >  2:706787/1‑149 (MQ=255)
                 aCTCTGCCTGCCAACCGCGCAATCTCTATTGCTCCAGATTTTGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCAc                                                                                                                                   <  4:518704/149‑1 (MQ=255)
                  ctctGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACg                                                                                                                                  >  2:801957/1‑149 (MQ=255)
                  ctctGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACg                                                                                                                                  <  1:339435/149‑1 (MQ=255)
                         tgccAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAg                                                                                                                           >  1:633676/1‑149 (MQ=255)
                         tgccAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCAATTACACCATTCTCGGCACGGTAAAAg                                                                                                                           >  4:411785/1‑149 (MQ=255)
                           ccAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGt                                                                                                                         >  1:301284/1‑149 (MQ=255)
                                cgcgCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGa                                                                                                                    >  3:206135/1‑149 (MQ=255)
                                 gcgcAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAg                                                                                                                   >  2:2942/1‑149 (MQ=255)
                                     attctcTATGGCACCAGTTTCGAACAATGCGCTGGTGCAAATCGCCGGCCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAGCGTAATGCAGCGTATCGGCGTGACCCATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGa                                                                                                              <  1:105616/147‑1 (MQ=255)
                                     aaTCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTg                                                                                                               <  4:827442/149‑1 (MQ=255)
                                      aTCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGa                                                                                                              <  1:117013/149‑1 (MQ=255)
                                      aTCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGa                                                                                                              <  3:836942/149‑1 (MQ=255)
                                      aTCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGa                                                                                                              <  1:561910/149‑1 (MQ=255)
                                      aTCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCAATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGa                                                                                                              <  3:411785/149‑1 (MQ=255)
                                      aTCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTGGAAAGCGTAATGCAGCGTATCGGCGGGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGa                                                                                                              >  2:432372/1‑149 (MQ=255)
                                         tctATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAgct                                                                                                           >  3:838276/1‑149 (MQ=255)
                                           tATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGTCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAgctgc                                                                                                         >  4:298556/1‑149 (MQ=255)
                                             ttGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATGACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGctgctg                                                                                                       >  4:569104/1‑149 (MQ=255)
                                                 aCCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCt                                                                                                   >  3:775853/1‑149 (MQ=255)
                                                 aCCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCt                                                                                                   <  1:258735/149‑1 (MQ=255)
                                                  ccAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCtt                                                                                                  >  1:504555/1‑149 (MQ=255)
                                                  ccAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGCTTGAGAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAGAAGGTGCGGAGCTTGAGCTGCTGCGCtt                                                                                                  <  2:627378/149‑1 (MQ=255)
                                                    aGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTa                                                                                                >  3:710058/1‑149 (MQ=255)
                                                    aGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTa                                                                                                <  1:333622/149‑1 (MQ=255)
                                                    aGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTa                                                                                                <  2:552942/149‑1 (MQ=255)
                                                     gATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTAc                                                                                               <  1:741076/149‑1 (MQ=255)
                                                           gACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATcc                                                                                         <  3:92135/149‑1 (MQ=255)
                                                                 tGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAgctgct                                                                                                        >  3:671249/1‑129 (MQ=255)
                                                                 tGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAgctgct                                                                                                        <  4:671249/129‑1 (MQ=255)
                                                                  gcgcTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTAt                                                                                   <  2:256716/149‑1 (MQ=255)
                                                                   cgcTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATg                                                                                  <  2:547946/149‑1 (MQ=255)
                                                                    gcTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATgg                                                                                 >  1:815462/1‑149 (MQ=255)
                                                                     cTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATggg                                                                                <  3:179504/149‑1 (MQ=255)
                                                                         tGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTc                                                                            <  3:603512/149‑1 (MQ=255)
                                                                          gCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCg                                                                           >  2:680191/1‑149 (MQ=255)
                                                                           cAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGa                                                                          >  1:26040/1‑149 (MQ=255)
                                                                            aGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAGTGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGACTTCGAc                                                                         >  2:55289/1‑149 (MQ=255)
                                                                             gATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACg                                                                        >  2:325823/1‑149 (MQ=255)
                                                                              aTCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGt                                                                       >  2:315522/1‑149 (MQ=255)
                                                                                           ccGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCAATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTc                                                          <  1:706787/149‑1 (MQ=255)
                                                                                            cGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCg                                                         <  4:451291/149‑1 (MQ=255)
                                                                                            cGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCg                                                         <  3:457539/149‑1 (MQ=255)
                                                                                            cGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCg                                                         <  1:876514/149‑1 (MQ=255)
                                                                                              tGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGc                                                       <  2:633676/149‑1 (MQ=255)
                                                                                               gATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCg                                                      <  2:504555/149‑1 (MQ=255)
                                                                                                    tGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCAc                                                 <  4:775853/149‑1 (MQ=255)
                                                                                                      gCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGt                                               >  1:822006/1‑149 (MQ=255)
                                                                                                      gCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGt                                               <  4:717742/149‑1 (MQ=255)
                                                                                                            gATCTGGTTGAAAGCGTAATGCAGCGTATTGGCGAGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCCCGTTACCCt                                         >  2:31263/1‑149 (MQ=255)
                                                                                                            gATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCt                                         <  1:801957/149‑1 (MQ=255)
                                                                                                             aTCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTg                                        >  2:128997/1‑149 (MQ=255)
                                                                                                              tCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTgg                                       >  3:157220/1‑149 (MQ=255)
                                                                                                              tCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTgg                                       >  1:866338/1‑149 (MQ=255)
                                                                                                                tGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGAt                                     >  1:633245/1‑149 (MQ=255)
                                                                                                                tGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGAt                                     >  1:668416/1‑149 (MQ=255)
                                                                                                                tGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGAt                                     >  1:747140/1‑149 (MQ=255)
                                                                                                                 ggTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATg                                    >  2:431731/1‑149 (MQ=255)
                                                                                                                  gTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGc                                   >  1:323515/1‑149 (MQ=255)
                                                                                                                   ttGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTaaaa                                                                                                                            <  1:288646/59‑1 (MQ=255)
                                                                                                                   ttGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTaaaa                                                                                                                            >  2:288646/1‑59 (MQ=255)
                                                                                                                    tGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCg                                 <  3:416440/149‑1 (MQ=255)
                                                                                                                            tAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGt                         <  4:185134/149‑1 (MQ=255)
                                                                                                                            gaaTGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGa                                      <  4:36582/135‑1 (MQ=255)
                                                                                                                            gaaTGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGa                                      >  3:36582/2‑136 (MQ=255)
                                                                                                                                gCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCg                     >  1:198650/1‑149 (MQ=255)
                                                                                                                                 cAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGt                    <  1:478622/149‑1 (MQ=255)
                                                                                                                                 cAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGt                    >  4:663346/1‑149 (MQ=255)
                                                                                                                                  aGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGtt                   >  1:380838/1‑149 (MQ=255)
                                                                                                                                  aGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGtt                   <  1:325823/149‑1 (MQ=255)
                                                                                                                                  aGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGtt                   <  3:858638/149‑1 (MQ=255)
                                                                                                                                   gCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTc                  >  3:668920/1‑149 (MQ=255)
                                                                                                                                                aCCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTgg     >  1:49688/1‑149 (MQ=255)
                                                                                                                                                 ccGGTTACACCGTTGTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCGTCCGTTTAAGGGCTTAGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGc    >  2:168533/1‑149 (MQ=255)
                                                                                                                                                   gATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCa  >  4:806302/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCA  >  NC_000913/23109‑23403

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