Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,082,631 C→T G434S (GGC→AGC)  speA ← biosynthetic arginine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,082,6310CT99.6% 1054.8 / ‑5.3 271G434S (GGC→AGC) speAbiosynthetic arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (145/125);  total (146/125)
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.38e-01

Only 100 of 275 total aligned reads displayed.

CCAGCTGCTTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGG  >  NC_000913/3082541‑3082720
                                                                                          |                                                                                         
ccAGCTGCTTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGa                                                                                         <  2:1370069/93‑1 (MQ=255)
ccAGCTGCTTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGa                                                                                         <  1:3628429/93‑1 (MQ=255)
ccAGCTGCTTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGa                                                                                         >  2:5473458/1‑93 (MQ=255)
 cAGCTGCTTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGat                                                                                        >  1:1418069/1‑93 (MQ=255)
  aGCTGCTTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGata                                                                                       >  2:2532888/1‑93 (MQ=255)
  aGCTGCTTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGata                                                                                       >  2:2099538/1‑93 (MQ=255)
    ctgctTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATg                                                                                     >  1:3445484/1‑93 (MQ=255)
     tgctTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGa                                                                                    <  1:752655/93‑1 (MQ=255)
         tttGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGt                                                                                <  2:4275725/93‑1 (MQ=255)
          ttGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTc                                                                               <  1:5697355/93‑1 (MQ=255)
           tGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCg                                                                              <  1:3073432/93‑1 (MQ=255)
             cACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTg                                                                            >  1:5978222/1‑93 (MQ=255)
              aCTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTga                                                                           >  1:4801979/1‑92 (MQ=255)
              aCTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGc                                                                           <  1:5068924/93‑1 (MQ=255)
                ttCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAg                                                                         >  2:5135816/1‑93 (MQ=255)
                  cATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATTTCGTGCAGAt                                                                       >  2:2349734/1‑93 (MQ=255)
                   aTGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATc                                                                      <  1:2726048/93‑1 (MQ=255)
                    tGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATcc                                                                     <  1:4882184/93‑1 (MQ=255)
                     ggCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCa                                                                    >  2:4141124/1‑93 (MQ=255)
                       cacaTGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATc                                                                  <  1:4586731/93‑1 (MQ=255)
                       cacaTGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATc                                                                  <  2:4644280/93‑1 (MQ=255)
                       cacaTGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATc                                                                  >  2:1005124/1‑93 (MQ=255)
                          aTGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGa                                                               >  1:4943102/1‑93 (MQ=255)
                          aTGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGa                                                               <  1:6237283/93‑1 (MQ=255)
                             cTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTg                                                            <  2:1530734/93‑1 (MQ=255)
                              tCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGt                                                           >  2:5855979/1‑93 (MQ=255)
                               cAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTc                                                          >  2:1726148/1‑93 (MQ=255)
                               cAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTc                                                          >  2:5837150/1‑93 (MQ=255)
                                aaaTAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCg                                                         >  2:1113016/1‑93 (MQ=255)
                                aaaTAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCg                                                         >  1:3890169/1‑93 (MQ=255)
                                 aaTAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCgt                                                        <  2:3429989/93‑1 (MQ=255)
                                 aaTAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCgt                                                        >  2:1943381/1‑93 (MQ=255)
                                 aaTAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCgt                                                        >  2:1618361/1‑93 (MQ=255)
                                  aTAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCgtg                                                       >  2:6097189/1‑93 (MQ=255)
                                   tAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCgtgt                                                      >  1:2798318/1‑93 (MQ=255)
                                   tAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCgtgt                                                      <  2:3009435/93‑1 (MQ=255)
                                      aGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAAc                                                   <  1:2070743/93‑1 (MQ=255)
                                      aGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAAc                                                   >  1:6158692/1‑93 (MQ=255)
                                      aGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAAc                                                   >  2:4689936/1‑93 (MQ=255)
                                       gctgctCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAAcc                                                  <  2:4388012/93‑1 (MQ=255)
                                       gctgctCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAAcc                                                  >  1:5623741/1‑93 (MQ=255)
                                       gctgctCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAAcc                                                  <  1:5228136/93‑1 (MQ=255)
                                       gctgctCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGCAAcc                                                  >  2:5064863/1‑93 (MQ=255)
                                         tgctCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCAt                                                >  2:1826616/1‑93 (MQ=255)
                                         tgctCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCAt                                                >  1:1971345/1‑93 (MQ=255)
                                          gctCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCAtt                                               <  2:3956428/93‑1 (MQ=255)
                                          gctCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCAtt                                               >  1:4216680/1‑93 (MQ=255)
                                           ctCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTc                                              <  1:3486025/93‑1 (MQ=255)
                                           ctCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTc                                              <  1:1801432/93‑1 (MQ=255)
                                           ctCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTc                                              <  1:1672866/93‑1 (MQ=255)
                                             cAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCAc                                            <  1:3946022/93‑1 (MQ=255)
                                              aGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACg                                           <  2:286284/93‑1 (MQ=255)
                                                cccATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCa                                         <  1:5567446/93‑1 (MQ=255)
                                                  cATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCaga                                       <  2:5989430/93‑1 (MQ=255)
                                                  cATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCaga                                       >  2:2089909/1‑93 (MQ=255)
                                                  cATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCaga                                       <  1:257318/93‑1 (MQ=255)
                                                   aTGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCagag                                      >  1:2206494/1‑93 (MQ=255)
                                                    tGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAgaga                                     >  2:6478105/1‑93 (MQ=255)
                                                    tGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAgaga                                     >  2:482646/1‑93 (MQ=255)
                                                     gCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGaa                                    >  1:645157/1‑93 (MQ=255)
                                                        cGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACgg                                 >  1:346175/1‑93 (MQ=255)
                                                         gTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGc                                >  1:811278/1‑93 (MQ=255)
                                                          ttCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCg                               >  2:5799898/1‑93 (MQ=255)
                                                           tCTTGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGa                              <  2:1576913/93‑1 (MQ=255)
                                                              tGCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGtt                           >  1:5967191/1‑93 (MQ=255)
                                                               gCAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTc                          <  1:6305394/93‑1 (MQ=255)
                                                                cAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTcc                         <  1:6410791/93‑1 (MQ=255)
                                                                cAGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTcc                         <  1:4529185/93‑1 (MQ=255)
                                                                 aGGCTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTccc                        <  2:3307233/93‑1 (MQ=255)
                                                                    cTAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGt                     <  1:3751167/93‑1 (MQ=255)
                                                                     tAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGtt                    <  1:5415107/93‑1 (MQ=255)
                                                                     tAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGtt                    <  1:1071841/93‑1 (MQ=255)
                                                                     tAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGtt                    <  1:4296175/93‑1 (MQ=255)
                                                                     tAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGtt                    <  1:6034861/93‑1 (MQ=255)
                                                                     tAAAGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGtt                    <  2:3176577/93‑1 (MQ=255)
                                                                       aaGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGGTCg                  >  1:5869166/1‑93 (MQ=255)
                                                                        aGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCgg                 >  2:1301198/1‑92 (MQ=255)
                                                                        aGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGt                 <  2:6467401/93‑1 (MQ=255)
                                                                        aGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGt                 >  1:1713777/1‑93 (MQ=255)
                                                                        aGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGt                 >  2:596514/1‑93 (MQ=255)
                                                                        aGATGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGt                 >  2:6342621/1‑93 (MQ=255)
                                                                          aTGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGc               <  1:5099545/93‑1 (MQ=255)
                                                                          aTGCCGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGc               >  1:2808405/1‑93 (MQ=255)
                                                                             ccGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCAtc            <  2:3874046/93‑1 (MQ=255)
                                                                              cGGAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCAtct           >  1:2268389/1‑93 (MQ=255)
                                                                               ggAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCAtctc          >  2:3720841/1‑93 (MQ=255)
                                                                                gAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTcc         <  1:4040608/93‑1 (MQ=255)
                                                                                gAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTcc         <  2:531499/93‑1 (MQ=255)
                                                                                gAAGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTcc         <  2:3626762/93‑1 (MQ=255)
                                                                                 aaGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCt        <  2:705951/93‑1 (MQ=255)
                                                                                 aaGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCt        >  1:1121213/1‑93 (MQ=255)
                                                                                 aaGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCt        <  2:4513693/93‑1 (MQ=255)
                                                                                 aaGAGTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCt        <  2:3229511/93‑1 (MQ=255)
                                                                                  agagTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTg       >  2:5564662/1‑93 (MQ=255)
                                                                                  agagTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTg       >  2:4905521/1‑93 (MQ=255)
                                                                                    agTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGcc     <  2:3542580/93‑1 (MQ=255)
                                                                                     gTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCa    >  1:4683301/1‑93 (MQ=255)
                                                                                     gTAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCa    >  1:301651/1‑93 (MQ=255)
                                                                                      tAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAg   >  2:2777117/1‑93 (MQ=255)
                                                                                      tAGCTGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAg   <  2:802878/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCAGCTGCTTTTGCACTTCATGGCACATGCTCAAATAAAGCTGCTCAGCCCATGCACGTTCTTGCAGGCTAAAGATGCCGGAAGAGTAGCCGATATGAATGTCGTGCAGATCCATCTGACTGTCGTGTAACCATTCACGCAGAGAACGGCGAGTTCCCGGTTCGTGCATCTCCTGCCAGG  >  NC_000913/3082541‑3082720

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