Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,071,673 C→T 6.2% P31P (CCG→CCA)  rutD ← putative reactive intermediate detoxifying aminoacrylate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,071,6730CT 6.2% 203.0 / 2.9 65P31P (CCG→CCA) rutDputative reactive intermediate detoxifying aminoacrylate hydrolase
Reads supporting (aligned to +/- strand):  ref base C (33/27);  new base T (2/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

CTCAATCCCTGCGGCTACCAGCGCCTGATGCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTC  >  NC_000913/1071526‑1071817
                                                                                                                                                   |                                                                                                                                                
cTCACTCCCTGCGGCTACCAGCGCCTGATGCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGc                                                                                                                                                  <  5:118630/148‑1 (MQ=255)
cTCAATCCCTGCGGCTACCAGCGCCTGATGCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCg                                                                                                                                                 <  1:330488/149‑1 (MQ=255)
     tCCCTGCGGCTACCAGCGCCTGATGCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAAc                                                                                                                                            <  7:391384/149‑1 (MQ=255)
       ccTGCGGCTACCAGCGCCTGATGCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCa                                                                                                                                          <  2:211197/149‑1 (MQ=255)
       ccTGCGGCTACCAGCGCCTGATGCAGTTCCGCTGCCATCGGGGCGATACTGTAATCTGCTGCCAGCGTGTCGGGATTATTGCCGGCGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAAc                                                                                                                                            <  1:105552/147‑1 (MQ=255)
       ccTGCGGCTACCAGCGCCTGATGCAGTTACGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGGCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGAGACTCCTGCTCCAGCACCGACAGGTGCGGTAAc                                                                                                                                            >  5:74226/1‑147 (MQ=255)
                aCCAGCGCCTGATGCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGTGGTAACCAGTAACTGcc                                                                                                                                 >  5:365391/1‑149 (MQ=255)
                  cAGCGCCTGATGCAGTTCCGCTGCCATCTGGGCGATACTGTAAGCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGATGGTCGTAACAGACTACCTGATAGTCATGCTCCAGTACCGCCACTTGCGGTAACCAGTAACTGCCGc                                                                                                                               >  1:106781/1‑149 (MQ=255)
                        cTGATGCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTAcccc                                                                                                                         >  8:413193/1‑149 (MQ=255)
                           aTGCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTAccccc                                                                                                                        <  2:20699/147‑1 (MQ=255)
                             gCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTAcc                                                                                                                           >  8:199709/1‑142 (MQ=255)
                             gCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTAcc                                                                                                                           <  7:199709/142‑1 (MQ=255)
                                   ccGCTGCCATCTGTGCGATACTGTAATCTTCTTCCAGCGTGTCGGGATTATTGCGGGTGCCGCGCTGGTCTTAACAGACTACCTGATACTCCTTCTCCAGCACCGCCAGTTGCGGTAAACAGTAACTGCAGCTACCCCCAAGACCCGa                                                                                                               <  6:193722/148‑1 (MQ=255)
                                       tGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATc                                                                                                           >  2:327167/1‑148 (MQ=255)
                                         ccATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGTGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAac                                                                                                        <  6:365391/149‑1 (MQ=255)
                                         ccATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAac                                                                                                        >  4:294339/1‑149 (MQ=255)
                                                        tGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGTGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCa                                                                                         >  3:92344/1‑149 (MQ=255)
                                                        tGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCa                                                                                         >  2:300753/1‑149 (MQ=255)
                                                            aTCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCa                                                                                      <  7:413193/148‑1 (MQ=255)
                                                               ttctGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCa                                                                                   >  5:386957/1‑148 (MQ=255)
                                                                 ctGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGc                                                                                    <  1:50570/145‑1 (MQ=255)
                                                                      aGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAggggg                                                                           >  6:257547/1‑149 (MQ=255)
                                                                       gCGTTTCGGGATTATTGACGGTGCCGCTCTGGTCGTAACATACTACCTGATAATCCTGCTCCAGCACCGCCAGTTGCGGTAAATAGTAATTGACTCTAACCCCAAGAAACGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGa                                                                          <  6:386957/149‑1 (MQ=255)
                                                                       gCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTACCCAGTAACTGCCGCTACCCCCCAGCCCCGAAATCAACCCCACGACGGGCGCAGCAGCATAAgggg                                                                            >  2:381060/1‑147 (MQ=255)
                                                                         gtgtCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAgg                                                                        <  1:327167/149‑1 (MQ=255)
                                                                            tCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGt                                                                       >  6:86303/1‑147 (MQ=255)
                                                                            tCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGt                                                                       >  7:188250/1‑147 (MQ=255)
                                                                            tCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTg                                                                      <  1:300753/148‑1 (MQ=255)
                                                                              gggATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGTGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTgaga                                                                   <  4:92344/149‑1 (MQ=255)
                                                                                  ttattGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTg                                                                >  4:123034/1‑148 (MQ=255)
                                                                                     ttGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGaaa                                                             <  3:294339/148‑1 (MQ=255)
                                                                                     ttGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGaaa                                                             <  5:86303/148‑1 (MQ=255)
                                                                                       gCCTGTGCCGCGCTGGTCGTAACAGACTACCTGATAGTCCTGCTCCAGTTCAGCGAGTTGAGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATGATCACCAATACGGGCTCATCAGCATAAGGGGGAGGTGAGAGTGaaatt                                                           <  2:106781/148‑3 (MQ=255)
                                                                                             gCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATcg                                                    >  1:264486/1‑149 (MQ=255)
                                                                                               cgcgCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATcgcg                                                  >  1:342834/1‑149 (MQ=255)
                                                                                                    tGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGc                                                <  4:308341/146‑1 (MQ=255)
                                                                                                    tGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGc                                                >  3:308341/1‑146 (MQ=255)
                                                                                                            aCAGCCTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAACTCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGc                                      <  1:381060/148‑1 (MQ=255)
                                                                                                              aGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGa                                    >  4:205435/1‑148 (MQ=255)
                                                                                                                  tACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCaa                                                                                                         <  7:153135/75‑1 (MQ=255)
                                                                                                                  tACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCaa                                                                                                         >  8:153135/1‑75 (MQ=255)
                                                                                                                     cTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGgc                                                                                            <  1:184767/85‑1 (MQ=255)
                                                                                                                     cTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGgc                                                                                            >  2:184767/1‑85 (MQ=255)
                                                                                                                        aTACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCCCTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAAtt                         >  1:222559/1‑149 (MQ=255)
                                                                                                                        aTACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCaa                           <  2:264486/147‑1 (MQ=255)
                                                                                                                           cTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATtgtg                      <  8:188250/149‑1 (MQ=255)
                                                                                                                            tCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGACATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTgg                     >  7:381598/1‑149 (MQ=255)
                                                                                                                            tCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAggg                                                                             <  5:314875/93‑1 (MQ=255)
                                                                                                                            tCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAggg                                                                             >  6:314875/1‑93 (MQ=255)
                                                                                                                            tCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACt                                          <  2:173943/128‑1 (MQ=255)
                                                                                                                            tCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACt                                          >  1:173943/1‑128 (MQ=255)
                                                                                                                            tCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTgg                     >  6:294152/1‑149 (MQ=255)
                                                                                                                              cTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGc                    >  6:218876/1‑148 (MQ=255)
                                                                                                                              cTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACATGGCGATATCCGCAATTGTTGc                    >  6:336385/1‑148 (MQ=255)
                                                                                                                               tGCTCCCGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGt                                                                 <  4:10690/102‑1 (MQ=255)
                                                                                                                               tGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGa                  >  4:173917/1‑149 (MQ=255)
                                                                                                                                gCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCTTCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGAt                 >  6:203644/1‑149 (MQ=255)
                                                                                                                                 cTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGAtt                <  5:218876/149‑1 (MQ=255)
                                                                                                                                    cAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGAc                                                                                                                   >  6:234993/1‑47 (MQ=255)
                                                                                                                                    cAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGAc                                                                                                                   <  5:234993/47‑1 (MQ=255)
                                                                                                                                         ccGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGTTTTCCACCAg        >  1:153242/1‑149 (MQ=255)
                                                                                                                                              aGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGAttt               <  4:103251/137‑1 (MQ=255)
                                                                                                                                              aGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGAttt               >  3:103251/1‑137 (MQ=255)
                                                                                                                                               gTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGt   >  6:150869/1‑148 (MQ=255)
                                                                                                                                               gTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTc  <  3:205435/149‑1 (MQ=255)
                                                                                                                                                ttGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGt   >  2:128109/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CTCAATCCCTGCGGCTACCAGCGCCTGATGCAGTTCCGCTGCCATCTGGGCGATACTGTAATCTTCTGCCAGCGTGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTC  >  NC_000913/1071526‑1071817

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