Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,263,132 G→C 6.7% S425R (AGC→AGG)  PP_1105 ← ATP‑dependent DNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,263,1320GC 6.7% 228.8 / 4.8 71S425R (AGC→AGG) PP_1105ATP‑dependent DNA ligase
Reads supporting (aligned to +/- strand):  ref base G (38/28);  new base C (2/3);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCG  >  NC_002947/1262987‑1263280
                                                                                                                                                 |                                                                                                                                                    
tGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCt                                                                                                                                                    >  8:346060/1‑148 (MQ=255)
         aaTAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGt                                                                                                                                            <  7:248853/147‑1 (MQ=255)
           tAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACCCTTTCAGCGTCGAt                                                                                                                                        <  8:49828/149‑1 (MQ=255)
             ggCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATc                                                                                                                                       <  6:309328/148‑1 (MQ=255)
              gCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATctt                                                                                                                                     >  8:38414/1‑149 (MQ=255)
               ccTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGCCAAAGGTGTAGTCGCTGTACAGACTGGCACTGCGGCCGTGGCCGCGTTCGGCATAGATCAGCACCGCATCGACCCTTTCAGCGTCGAtgcgg                                                                                                                                    <  2:325792/149‑5 (MQ=255)
                cTTGGCAAAGGGCACCAACGCCCGCCCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCt                                                                                                                                                    <  5:326793/132‑1 (MQ=255)
                cTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCt                                                                                                                                                    >  6:326793/1‑132 (MQ=255)
                         gggCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACCCTGAAAGGGTCGATCttccacttccac                                                                                                                           >  5:55758/1‑148 (MQ=255)
                          ggCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCttccacttccac                                                                                                                           >  6:368532/1‑147 (MQ=255)
                          ggCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTccacca                                                                                                                         <  1:262581/149‑1 (MQ=255)
                            cACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACCCTTTCAGCGTCGAt                                                                                                                                        <  4:361162/132‑1 (MQ=255)
                            cACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACCCTTTCAGCGTCGAt                                                                                                                                        >  3:361162/1‑132 (MQ=255)
                             aCCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGAt                                                                                                                      <  8:93295/149‑1 (MQ=255)
                                 aCGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGGCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGccc                                                                                                                  <  4:101711/149‑1 (MQ=255)
                                   gCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGAGCTTCCACTTCCACCAGATGCCCAt                                                                                                                >  2:331403/1‑149 (MQ=255)
                                        aCCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCt                                                                                                           >  5:287693/1‑149 (MQ=255)
                                            ggggTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGt                                                                                                       <  7:346060/149‑1 (MQ=255)
                                             gggTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTg                                                                                                      >  4:17693/1‑149 (MQ=255)
                                                   cTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGcccc                                                                                                >  2:199787/1‑149 (MQ=255)
                                                     gCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAg                                                                                                                                               <  3:282446/100‑1 (MQ=255)
                                                     gCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAg                                                                                                                                               >  4:282446/1‑100 (MQ=255)
                                                           gCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTa                                                                                        >  4:178803/1‑149 (MQ=255)
                                                               cccATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGAAGCCCATGTCCTTGGTGCGCCCCACACCGTACAgg                                                                                    >  6:144143/1‑149 (MQ=255)
                                                                ccATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAggg                                                                                   <  6:287693/149‑1 (MQ=255)
                                                                 cATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTc                                                                                                                                           <  7:20737/92‑1 (MQ=255)
                                                                 cATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTc                                                                                                                                           >  8:20737/1‑92 (MQ=255)
                                                                  aTACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTg                                                                                                         <  3:93403/125‑1 (MQ=255)
                                                                  aTACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCAACCAGATGCCCATGTCCTTg                                                                                                         >  4:93403/1‑125 (MQ=255)
                                                                     cGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTcg                                                                              <  3:17693/149‑1 (MQ=255)
                                                                      ggCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTcgc                                                                             <  4:187728/149‑1 (MQ=255)
                                                                      ggCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTcgc                                                                             >  6:301682/1‑149 (MQ=255)
                                                                        cAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGt                                                                           >  1:161206/1‑149 (MQ=255)
                                                                         aaaGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGtt                                                                          <  5:59592/149‑1 (MQ=255)
                                                                           aGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGc                                                                        <  5:144143/149‑1 (MQ=255)
                                                                               gtAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGt                                                                                                       <  4:301233/114‑1 (MQ=255)
                                                                               gtAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGt                                                                                                       >  3:301233/1‑114 (MQ=255)
                                                                                tAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAg                                                                   <  5:368532/149‑1 (MQ=255)
                                                                                  gTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGc                                                                  >  7:172080/1‑148 (MQ=255)
                                                                                       tGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGc                                                            >  6:187815/1‑149 (MQ=255)
                                                                                         tACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGTGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGccc                                                          >  3:244893/1‑149 (MQ=255)
                                                                                           cAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCtt                                                        >  2:215683/1‑149 (MQ=255)
                                                                                            aGACTGGCCCTGCGGCCGTGGCCGCGTTGGGTATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTc                                                       >  7:195847/1‑149 (MQ=255)
                                                                                            aGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCAt                                                                                                                <  5:329889/92‑1 (MQ=255)
                                                                                            aGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCAt                                                                                                                >  6:329889/1‑92 (MQ=255)
                                                                                            aGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTc                                                       >  2:192396/1‑149 (MQ=255)
                                                                                                  gCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCCGCCCTTCCACGcc                                                 >  4:352003/1‑149 (MQ=255)
                                                                                                        cGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCtt                                                                      <  3:109788/122‑1 (MQ=255)
                                                                                                        cGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCtt                                                                      >  4:109788/1‑122 (MQ=255)
                                                                                                          gCCGTGGCGGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGTGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGc                                         <  6:55758/149‑1 (MQ=255)
                                                                                                             acggCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGc                                                                                  >  2:336227/3‑105 (MQ=255)
                                                                                                             acggCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGc                                                                                  <  1:336227/103‑1 (MQ=255)
                                                                                                                 ccGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGt                                                                                                                                            <  8:274942/43‑1 (MQ=255)
                                                                                                                 ccGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGt                                                                                                                                            >  7:274942/1‑43 (MQ=255)
                                                                                                                 ccGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTg                                  >  3:183132/1‑149 (MQ=255)
                                                                                                                         ggCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCa                                                              >  3:176497/1‑113 (MQ=255)
                                                                                                                         ggCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCa                                                              <  4:176497/113‑1 (MQ=255)
                                                                                                                           cATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTcc                                                      <  3:309171/119‑1 (MQ=255)
                                                                                                                           cATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTgg                        >  3:126147/1‑149 (MQ=255)
                                                                                                                           cATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTAAGAATCAGCCCTTCCACGCCCAGGATGCGTGACTGTTCCCGCTgg                        >  6:318769/1‑149 (MQ=255)
                                                                                                                           cATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCACATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTcc                                                      >  4:309171/1‑119 (MQ=255)
                                                                                                                            aTAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCACCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGc                       >  5:247674/1‑149 (MQ=255)
                                                                                                                                aTCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCggg                    >  4:164610/1‑148 (MQ=255)
                                                                                                                                  cAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGa                                     >  2:256606/1‑129 (MQ=255)
                                                                                                                                  cAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGa                                     <  1:256606/129‑1 (MQ=255)
                                                                                                                                   aGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAg                <  5:187815/149‑1 (MQ=255)
                                                                                                                                     cACCGCATCGACGCTGAAAGGGACGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGt              >  7:238570/1‑149 (MQ=255)
                                                                                                                                       ccGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCg            <  5:301682/149‑1 (MQ=255)
                                                                                                                                          cATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGcc         >  8:259732/1‑149 (MQ=255)
                                                                                                                                           aTCGACGATGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCACCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGACCAGGCTGCTTGACTGATCGCGCTCGCGGGCAACGTCGGcc         >  8:343406/1‑148 (MQ=255)
                                                                                                                                            tCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCa       <  8:238570/149‑1 (MQ=255)
                                                                                                                                                 gCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCg  >  4:147783/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCG  >  NC_002947/1262987‑1263280

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: ‑