Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,673,161 G→A 6.5% G118D (GGC→GAC)  rhlE‑II → ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,673,1610GA 6.5% 179.8 / 2.3 61G118D (GGC→GAC) rhlE‑IIATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base G (22/35);  new base A (2/2);  total (24/37)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CCCGACAAATCACAGCGCCACGGCCCGCGCCAGCCTCGCGTGCTGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGTGGAAATCCTGGTACTTGATGAAGCCGATCGCA  >  NC_002947/5673013‑5673301
                                                                                                                                                    |                                                                                                                                            
cccGACAAATCACAGCGCCACGGCCCGCGCCAGCCTCGCGTGCTGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGc                                                                                                                                             <  1:1228915/150‑1 (MQ=255)
  cGACAAATCACAGCGCCACGGCCCGCGCCAGCCTCGCGTGCTGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCAt                                                                                                                                           <  2:1063649/150‑1 (MQ=255)
      aaaTCACAGCGCCACGGCCCGCGCCAGCCTCGCGTGCTGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAAc                                                                                                                                       >  2:986942/1‑150 (MQ=255)
      aaaTCACAGCGCCACGGCCCGCGCCAGCCTCGCGTGCTGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAAc                                                                                                                                       >  1:308153/1‑150 (MQ=255)
         tCACAGCGCCACGGCCCGCGCCAGCCTCGCGTGCTGATCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCg                                                                                                                                    >  2:225876/1‑150 (MQ=255)
         acacaGCGCCACGGCCCGCGCCAGCCTCGCGTGCTGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCg                                                                                                                                    <  1:1284344/149‑1 (MQ=255)
           acaGCGCCACGGCCCGCGCCAGCCTCGCGTGCTGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCa                                                                                                                                  <  2:172701/150‑1 (MQ=255)
             aGCGCCACGGCCCGCGCCAGCCTCGCGTGCTGGTCCTGACGCCGGCCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGa                                                                                                                                <  2:726775/150‑1 (MQ=255)
                gcCACGGCCCGCGCCAGCCTCGCGTGCTGGTCCTGACGCCGATCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCGTCGGCGGCGTTGGCATGAACCCGCAGATcc                                                                                                                             <  2:500654/150‑1 (MQ=255)
                  cACGGCCCGCGCCAGCCTCGCGTGCTGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAg                                                                                                                           >  1:1001471/1‑150 (MQ=255)
                  cACGGCCCGCGCCAGCCTCGCGTGCTGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAg                                                                                                                           >  2:725922/1‑150 (MQ=255)
                          gcgcCAGCCTCGCGTGCTGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGc                                                                                                                   <  1:932835/150‑1 (MQ=255)
                          gcgcCAGCCTCGCGTGCTGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGc                                                                                                                   >  1:1095754/1‑150 (MQ=255)
                             ccAGCCTCGCGTGCTGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCaa                                                                                                                <  1:725922/150‑1 (MQ=255)
                               aGCCTCGCGTGCTGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAg                                                                                                              <  1:249344/150‑1 (MQ=255)
                                   tCGCGTGCTGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGt                                                                                                          >  2:46445/1‑150 (MQ=255)
                                   tCGCGTGCTGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGt                                                                                                          >  1:646704/1‑150 (MQ=255)
                                   tCGCGTGCTGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGCTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGt                                                                                                          >  2:930953/1‑150 (MQ=255)
                                       gTGCTGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGAt                                                                                                      <  1:698453/150‑1 (MQ=255)
                                           tGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGc                                                                                                  <  1:46445/150‑1 (MQ=255)
                                            ggTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCt                                                                                                 <  1:225876/150‑1 (MQ=255)
                                             gTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTg                                                                                                <  2:308153/150‑1 (MQ=255)
                                                       cGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGcc                                                                                      <  2:511092/150‑1 (MQ=255)
                                                       cGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGcc                                                                                      <  1:345655/150‑1 (MQ=255)
                                                           ccGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAgg                                                                                  <  1:39169/150‑1 (MQ=255)
                                                                            aGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTc                                                                        >  1:609953/1‑143 (MQ=255)
                                                                            aGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTc                                                                        <  2:609953/143‑1 (MQ=255)
                                                                               tGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCg                                                              <  2:1095754/150‑1 (MQ=255)
                                                                                  aTGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGcc                                                           >  2:211977/1‑150 (MQ=255)
                                                                                     aCAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGACATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGAt                                                                                                      <  2:740797/104‑1 (MQ=255)
                                                                                     aCAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGACATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGAt                                                                                                      >  1:740797/1‑104 (MQ=255)
                                                                                     aCAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGACATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGAt                                                                                                      <  2:740450/104‑1 (MQ=255)
                                                                                     aCAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGACATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGAt                                                                                                      >  1:740450/1‑104 (MQ=255)
                                                                                         tttCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCcaaggcaa                                                    >  1:438445/1‑150 (MQ=255)
                                                                                            cAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTgg                                                                                               >  2:585824/1‑104 (MQ=255)
                                                                                            cAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTgg                                                                                               <  1:585824/104‑1 (MQ=255)
                                                                                                gTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGAc                                             <  1:565162/150‑1 (MQ=255)
                                                                                                    aTGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTgg                                         >  2:728055/1‑150 (MQ=255)
                                                                                                        ccGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGc                                                                                                  >  1:101867/1‑89 (MQ=255)
                                                                                                        ccGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGc                                                                                                  <  2:101867/89‑1 (MQ=255)
                                                                                                         cGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCAc                                    >  2:156891/1‑150 (MQ=255)
                                                                                                           tGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGt                                  <  2:524466/150‑1 (MQ=255)
                                                                                                             aCCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGTgg                                <  2:646704/150‑1 (MQ=255)
                                                                                                               cTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGTGGaa                              >  2:441590/1‑150 (MQ=255)
                                                                                                                tGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGTGGaaa                             <  2:1001471/150‑1 (MQ=255)
                                                                                                                  ccACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGTCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTgg                                         <  1:1085531/136‑1 (MQ=255)
                                                                                                                  ccACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGTCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTgg                                         >  2:1085531/1‑136 (MQ=255)
                                                                                                                      tGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGTGGAAATCCTgg                       <  2:756005/150‑1 (MQ=255)
                                                                                                                      tGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGTGGAAATCCTgg                       <  1:211977/150‑1 (MQ=255)
                                                                                                                      tGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGTGGAAATCCTgg                       <  2:481175/150‑1 (MQ=255)
                                                                                                                       ggTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCTAAAGGCGTGGATTTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGTGGAAATCCTGGt                      <  1:930953/150‑1 (MQ=255)
                                                                                                                        gTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGTGGAAATcc                          <  2:719052/145‑1 (MQ=255)
                                                                                                                        gTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTAGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGTGGAAATcc                          >  1:719052/1‑145 (MQ=255)
                                                                                                                          cAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGTGGAAATCCTGGTACt                   <  1:504336/150‑1 (MQ=255)
                                                                                                                            gcgcCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGTGGAAATCCTGGTACTtg                 >  1:1160359/1‑150 (MQ=255)
                                                                                                                              gcCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGTGGAAATCCTGGTACTtgat               <  2:1021195/150‑1 (MQ=255)
                                                                                                                                           gcggcgTTGGCATGAACCCGCAGATCTAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGt                                  <  1:919936/118‑1 (MQ=255)
                                                                                                                                           gcggcgTTGGCATGAACCCGCAGATCTAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGt                                  >  2:919936/1‑118 (MQ=255)
                                                                                                                                           gcggcgTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGTGGAAATCCTGGTACTTGATGAAGCCGATCGCa  <  1:986942/150‑1 (MQ=255)
                                                                                                                                                 ttGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGTGGAAATCCTGGTACTTGATGAAGCCGa       <  1:329862/139‑1 (MQ=255)
                                                                                                                                                 ttGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGTGGAAATCCTGGTACTTGATGAAGCCGa       >  2:329862/1‑139 (MQ=255)
                                                                                                                                                    gCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGAc                                                                     >  1:760521/1‑74 (MQ=255)
                                                                                                                                                    gCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGAc                                                                     <  2:760521/74‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
CCCGACAAATCACAGCGCCACGGCCCGCGCCAGCCTCGCGTGCTGGTCCTGACGCCGACCCGCGAACTGGCAGCCCAGGTGCATGACAGTTTCAAGGTCTATGCCCGTGACCTGCCACTGGTCAGCGCCTGCATCTTCGGCGGCGTTGGCATGAACCCGCAGATCCAGGCCATTGCCAAAGGCGTTGATGTGCTGGTCGCCTGCCCAGGCCGCCTGCTCGACCTCGCCGGCCAAGGCAAGGTCGACCTGGCCCACGTGGAAATCCTGGTACTTGATGAAGCCGATCGCA  >  NC_002947/5673013‑5673301

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