Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,623,199 A→G 6.1% intergenic (+35/‑211) PP_5517 → / → PP_2295 hypothetical protein/antirestriction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,623,1990AG 6.1% 190.8 / 2.2 65intergenic (+35/‑211)PP_5517/PP_2295hypothetical protein/antirestriction protein
Reads supporting (aligned to +/- strand):  ref base A (31/30);  new base G (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

TTCGGGGATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGTTAC  >  NC_002947/2623056‑2623345
                                                                                                                                               |                                                                                                                                                    
ttCGGGGATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCgg                                                                                                                                                 >  5:202111/1‑149 (MQ=255)
   ggggATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGcc                                                                                                                                              >  2:282028/1‑149 (MQ=255)
   ggggATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGc                                                                                                                                               <  6:202111/148‑1 (MQ=255)
      gATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTTGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCtt                                                                                                                                           >  5:297790/1‑149 (MQ=255)
       aTGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCCTACGCAATCACACCTCAAGGGGCCCAATCGGGCCCttt                                                                                                                                          <  5:278251/149‑1 (MQ=255)
            tgtgGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTcgc                                                                                                                                     <  8:332945/149‑1 (MQ=255)
            tgtgGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAAGGGCCCGATTGGGCCCcttgaggt                                                                                                                                     <  6:343044/149‑6 (MQ=255)
                   gtgtTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTa                                                                                                                              >  8:47433/1‑149 (MQ=255)
                              gAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACA‑‑TCGATTCTg                                                                                                                 <  2:324796/149‑1 (MQ=255)
                                    gTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGAAGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TAGATTCTGGAGAAt                                                                                                           >  7:315720/1‑149 (MQ=255)
                                           gTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAgg                                                                                                    >  7:315916/1‑149 (MQ=255)
                                            tCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAggg                                                                                                   >  2:76716/1‑149 (MQ=255)
                                             cACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAAGGGCCCGATTGGGCCCcttgag                                                                                                                                       >  7:20842/1‑111 (MQ=255)
                                             cACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAAGGGCCCGATTGGGCCCcttgag                                                                                                                                       <  8:20842/114‑4 (MQ=255)
                                              aCCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGc                                                                                                 >  3:347956/1‑149 (MQ=255)
                                                cttccttTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCt                                                                                                >  7:128278/1‑148 (MQ=255)
                                                      ttCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCt                                                                                                <  6:356362/142‑1 (MQ=255)
                                                      ttCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCt                                                                                                >  5:356362/1‑142 (MQ=255)
                                                         tcCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGCCTGTACG‑‑TCGAtt                                                                                                                    >  2:105700/1‑119 (MQ=255)
                                                         tcCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGCCTGTACG‑‑TCGAtt                                                                                                                    <  1:105700/119‑1 (MQ=255)
                                                                 gTCGAGCCTTCTACACCTACACCAAGTTTGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCcgc                                                                              <  6:297790/149‑1 (MQ=255)
                                                                 gTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCcgc                                                                              <  7:47433/149‑1 (MQ=255)
                                                                 gTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCcgc                                                                              <  8:128278/149‑1 (MQ=255)
                                                                  tCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGa                                                                                        <  5:12180/138‑1 (MQ=255)
                                                                  tCGAGCCTTCTACACCTACACAAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGa                                                                                        >  6:12180/1‑138 (MQ=255)
                                                                        cTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGa                                                                                    >  5:190242/1‑136 (MQ=255)
                                                                        cTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGa                                                                                    <  6:190242/136‑1 (MQ=255)
                                                                             acacctacaccAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAgcgcgc                                                                  >  4:205060/1‑149 (MQ=255)
                                                                                     accAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTa                                                          >  7:263215/1‑149 (MQ=255)
                                                                                       cAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACg                                                        <  3:205060/149‑1 (MQ=255)
                                                                                        aaGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGccc                                                                                          >  8:151023/1‑114 (MQ=255)
                                                                                        aaGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGccc                                                                                          <  7:151023/114‑1 (MQ=255)
                                                                                         aGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGaa                                                      <  8:263215/149‑1 (MQ=255)
                                                                                               tGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGCTTTGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAaca                                                  <  1:282028/149‑1 (MQ=255)
                                                                                                  ggAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCa                                             >  8:103924/1‑149 (MQ=255)
                                                                                                   gAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCAt                                            >  1:21967/1‑149 (MQ=255)
                                                                                                        gcTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTa                                                          >  2:843/1‑130 (MQ=255)
                                                                                                        gcTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTa                                                          <  1:843/130‑1 (MQ=255)
                                                                                                        gcTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTcc                                        >  3:89593/1‑148 (MQ=255)
                                                                                                         cTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTTGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCca                                       >  7:17054/1‑148 (MQ=255)
                                                                                                         cTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCcac                                      <  1:76716/149‑1 (MQ=255)
                                                                                                             tgCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAAt                                  <  8:315720/148‑1 (MQ=255)
                                                                                                             cGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAAt                                  <  7:247669/149‑1 (MQ=255)
                                                                                                                aaGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTcc                                        <  4:199784/140‑1 (MQ=255)
                                                                                                                aaGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTcc                                        >  3:199784/1‑140 (MQ=255)
                                                                                                                aaGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGa                               <  7:103924/149‑1 (MQ=255)
                                                                                                                  gTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGaaa                             <  8:315916/149‑1 (MQ=255)
                                                                                                                   tCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGaaa                             >  2:53932/1‑148 (MQ=255)
                                                                                                                         gCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGt                                               <  8:299976/124‑1 (MQ=255)
                                                                                                                         gCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGt                                               >  7:299976/1‑124 (MQ=255)
                                                                                                                                cacCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCAt               <  8:17054/149‑1 (MQ=255)
                                                                                                                                 acCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCa                      <  2:81572/141‑1 (MQ=255)
                                                                                                                                 acCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCa                      >  1:81572/1‑141 (MQ=255)
                                                                                                                                 acCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATg              >  5:354266/1‑149 (MQ=255)
                                                                                                                                   gaaaaagggCCCGATTGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATggg            >  1:254333/7‑149 (MQ=255)
                                                                                                                                   cTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATggg            <  3:273856/149‑1 (MQ=255)
                                                                                                                                   cTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATggg            >  4:273856/1‑149 (MQ=255)
                                                                                                                                   cTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGGAGAAGCGCGCCAAAATTACGAAAAGACGTCATTTCGACAAATGGCTATAAATCACTTGCTTggg            >  2:97150/1‑149 (MQ=255)
                                                                                                                                    tCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATgggg           >  7:133746/1‑149 (MQ=255)
                                                                                                                                           gCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGtt    <  4:347956/149‑1 (MQ=255)
                                                                                                                                            cccAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGaaataaat                        >  5:23934/1‑128 (MQ=255)
                                                                                                                                            cccAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGaaataaat                        <  6:23934/128‑1 (MQ=255)
                                                                                                                                             ccAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATggg            <  8:191238/139‑1 (MQ=255)
                                                                                                                                             ccAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATggg            >  7:191238/1‑139 (MQ=255)
                                                                                                                                             ccAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGTTAc  <  1:53932/149‑1 (MQ=255)
                                                                                                                                               aaTCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTAAGAAAACACGTCATTTCAACaaag                                  >  6:171624/1‑114 (MQ=255)
                                                                                                                                               aaTCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCAGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGaaat                            <  6:181338/121‑1 (MQ=255)
                                                                                                                                               aaTCGCGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCAGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGaaat                            >  5:181338/1‑121 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TTCGGGGATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACG‑‑TCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGTTAC  >  NC_002947/2623056‑2623345

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