Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 254,223 A→T 15.1% intergenic (+21/+36) prfH → / ← pepD pseudogene, RF‑1 domain family; probable peptide chain release factor/aminoacyl‑histidine dipeptidase (peptidase D)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913254,2230AT15.1% 257.2 / 32.8 99intergenic (+21/+36)prfH/pepDpseudogene, RF‑1 domain family; probable peptide chain release factor/aminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  ref base A (40/44);  new base T (8/7);  total (48/51)
Fisher's exact test for biased strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

GCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCG  >  NC_000913/254127‑254322
                                                                                                |                                                                                                   
gcgAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGc                                                                                                  <  1:799532/100‑1 (MQ=255)
  gAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGa                                                                                                <  1:1042991/100‑1 (MQ=255)
   aaTGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCg                                                                                                 >  2:91300/1‑98 (MQ=255)
   aaTGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGAt                                                                                               <  1:2505832/100‑1 (MQ=255)
   aaTGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGAt                                                                                               <  1:1968776/100‑1 (MQ=255)
   aaTGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGAt                                                                                               <  2:2232441/100‑1 (MQ=255)
    aTGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGAtt                                                                                              >  1:1164673/1‑100 (MQ=255)
     tGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTa                                                                                             <  2:830946/100‑1 (MQ=255)
     tGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTa                                                                                             <  1:865600/100‑1 (MQ=255)
     tGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTa                                                                                             >  1:951823/1‑100 (MQ=255)
      gTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTaa                                                                                            >  1:483440/1‑100 (MQ=255)
       ttCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTaaa                                                                                           >  1:41143/1‑100 (MQ=255)
        tCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAAc                                                                                          >  1:1397341/1‑100 (MQ=255)
         ccATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACg                                                                                         <  2:430836/100‑1 (MQ=255)
         ccATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACg                                                                                         >  1:1445814/1‑100 (MQ=255)
          cATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGc                                                                                        <  2:2463211/100‑1 (MQ=255)
             cACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCAt                                                                                     <  1:1581104/100‑1 (MQ=255)
               ccAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATcc                                                                                   <  1:1843154/100‑1 (MQ=255)
                cAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCg                                                                                  <  1:534034/100‑1 (MQ=255)
                 aGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCgg                                                                                 <  2:1194282/100‑1 (MQ=38)
                 aGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCgg                                                                                 >  2:1140051/1‑100 (MQ=255)
                 aGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCgg                                                                                 <  2:1026146/100‑1 (MQ=255)
                  gATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgc                                                                                <  1:865091/100‑1 (MQ=38)
                  gATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgc                                                                                >  2:307005/1‑100 (MQ=255)
                  gATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgc                                                                                >  1:2176795/1‑100 (MQ=255)
                   aTTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgca                                                                               <  1:596731/100‑1 (MQ=255)
                    ttGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCTCCTTCCGgcag                                                                              <  2:1341514/100‑1 (MQ=38)
                    ttGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgcag                                                                              <  2:704664/100‑1 (MQ=38)
                    ttGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgcag                                                                              <  1:149575/100‑1 (MQ=38)
                    ttGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgcag                                                                              <  2:1423462/100‑1 (MQ=38)
                     tGAACGTGGTAACCCGCTACGGATATTTACAGAGTTGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTAAATCTCCTTCCGgca                                                                               <  2:2138655/98‑1 (MQ=255)
                      gAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagca                                                                            >  1:1003918/1‑100 (MQ=255)
                      gAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAAATGCGGGAAGGCGATTAAAcacagtacag                                                                                 >  2:2447099/1‑86 (MQ=255)
                          gTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatc                                                                        <  2:1219327/100‑1 (MQ=255)
                              cAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaaat                                                                    >  1:117929/1‑100 (MQ=255)
                              cAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaaat                                                                    >  2:90307/1‑100 (MQ=255)
                              cAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaaat                                                                    <  1:1549999/100‑1 (MQ=255)
                              cAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaaat                                                                    <  1:640458/100‑1 (MQ=255)
                              cAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaaat                                                                    >  2:2406766/1‑100 (MQ=255)
                                 cccGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataat                                                                 >  1:2252952/1‑100 (MQ=255)
                                   cgcgACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa                                                               >  2:305295/1‑100 (MQ=255)
                                    gcgACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTAc                                                              <  1:1054509/100‑1 (MQ=255)
                                      gACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACtt                                                            >  1:855025/1‑100 (MQ=255)
                                      atctgACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACtt                                                            <  1:846140/96‑1 (MQ=255)
                                        cGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCg                                                          <  1:107061/100‑1 (MQ=255)
                                         ggACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGc                                                         >  1:572031/1‑100 (MQ=255)
                                          gACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACttagat                                                        >  2:846133/1‑96 (MQ=255)
                                          gACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGcc                                                        <  2:1691333/100‑1 (MQ=255)
                                          gACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGcc                                                        >  2:846243/1‑100 (MQ=255)
                                               ttACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGACAGCAAATCAAATAATTACTTCGCCGGAAt                                                   <  2:1940076/100‑1 (MQ=255)
                                                tACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtt                                                  <  1:305299/100‑1 (MQ=255)
                                                  cAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGACATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttc                                                <  1:41746/100‑1 (MQ=255)
                                                     ggATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttcttt                                             <  1:986421/100‑1 (MQ=255)
                                                      aattgCTTTTCTCGAAGCATAATGAAGGAGGAAACGTCCCTAAGGCGATTAAACGCCATCCGTCAGCAAATCAAATAATTACTTCGCCGGAAtttctttg                                            <  2:2359529/96‑2 (MQ=255)
                                                       aTGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttctttc                                            <  1:1577152/99‑1 (MQ=255)
                                                        tGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGTAAATCAAATAATTACTTCGCCGGAATTTCTTTcag                                          <  1:1978519/100‑1 (MQ=255)
                                                        tGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcag                                          <  2:651709/100‑1 (MQ=255)
                                                         ggCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagc                                         <  1:284699/100‑1 (MQ=255)
                                                         ggCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagc                                         <  1:222904/100‑1 (MQ=255)
                                                           cTTTTATCGAAGGATAATGAAGGATGAAACTGCTGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                                       >  2:116801/1‑100 (MQ=255)
                                                              ttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtc                                    <  1:634495/100‑1 (MQ=255)
                                                                aTCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcag                                  >  1:34828/1‑100 (MQ=255)
                                                                aTCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcag                                  <  1:709675/100‑1 (MQ=255)
                                                                aTCGAAGGATAATGAAGGATGAAACTACCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcag                                  >  2:238950/1‑100 (MQ=255)
                                                                 tCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagt                                 >  2:2336335/1‑100 (MQ=255)
                                                                 tCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagt                                 >  2:2335076/1‑100 (MQ=255)
                                                                   gAAGGATAATGAAGGATGAACCTGCCAGAAGGCGATGAAACGCCATCCGGAAGCAAATCAAATAATTACTACGGCGGAAGTTCTTTCAGCAGTacagt                                 >  2:1409458/1‑98 (MQ=255)
                                                                     aGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagc                             >  2:2519893/1‑100 (MQ=255)
                                                                     aGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagc                             >  2:1546225/1‑100 (MQ=255)
                                                                      ggATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagca                            >  1:1637944/1‑100 (MQ=255)
                                                                      ggATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCACCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagca                            >  2:1700622/1‑100 (MQ=255)
                                                                       gATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagcag                           <  2:151580/100‑1 (MQ=255)
                                                                          aaTGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAgtgt                        <  1:955772/100‑1 (MQ=255)
                                                                            tGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTcc                      >  2:2306081/1‑100 (MQ=255)
                                                                            tGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTcc                      <  1:279432/100‑1 (MQ=255)
                                                                            tGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTGCTTTCAGCAGTTCAGTCAGCAGTGTcc                      >  2:38545/1‑100 (MQ=255)
                                                                             gAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCa                     >  2:2033556/1‑100 (MQ=255)
                                                                              aaGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAg                    <  1:1690486/100‑1 (MQ=255)
                                                                              aaGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAg                    <  2:67067/100‑1 (MQ=255)
                                                                               aGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGt                   <  2:376901/100‑1 (MQ=255)
                                                                               aGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGt                   <  1:2059323/100‑1 (MQ=255)
                                                                                 gATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTaa                 >  1:1077724/1‑100 (MQ=255)
                                                                                  aTGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAAt                >  1:2557001/1‑100 (MQ=255)
                                                                                   tGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATg               >  1:87942/1‑100 (MQ=255)
                                                                                    gAAACTGCCGGAAGGCGATTAAACGCCACCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGa              <  1:1700634/100‑1 (MQ=255)
                                                                                      aaCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGAcn            <  1:1707587/100‑2 (MQ=255)
                                                                                      aaCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGAcc            <  2:2176780/100‑1 (MQ=255)
                                                                                       gcTGCCGGATGGCATTTAATCGCCTTCCGGCAGCAAATAAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTGAGTCAGCAGTGTCCAGTAATGACCt           >  2:1582333/2‑100 (MQ=37)
                                                                                       aCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCt           >  2:615516/1‑100 (MQ=255)
                                                                                        cTGCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTa          >  1:1178324/1‑100 (MQ=38)
                                                                                        cTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTa          >  1:1699065/1‑100 (MQ=255)
                                                                                         tGCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTAc         >  2:1593673/1‑100 (MQ=38)
                                                                                          gCCGGATGGCGTTTAATCGCCTTCCGGCAGTAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTaag        >  2:397917/1‑98 (MQ=38)
                                                                                          gCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg        >  2:1098548/1‑100 (MQ=38)
                                                                                          gCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg        >  1:2059114/1‑100 (MQ=38)
                                                                                          gCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg        >  2:889491/1‑100 (MQ=38)
                                                                                          gCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg        <  1:1404747/100‑1 (MQ=255)
                                                                                           ccGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGc       >  1:316484/1‑100 (MQ=255)
                                                                                            cGGATGGCGTTTAATCGCCTTCCGGCAGTAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCt      >  1:951544/1‑100 (MQ=38)
                                                                                                aGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCg  <  2:363381/100‑1 (MQ=255)
                                                                                                |                                                                                                   
GCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCG  >  NC_000913/254127‑254322

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