Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,058,461 A→T 5.2% E553V (GAG→GTG)  recG → junction‑specific ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,058,4610AT 5.2% 241.0 / 2.5 77E553V (GAG→GTG) recGjunction‑specific ATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base A (32/40);  new base T (2/2);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

AGCACCTATGAGGAGCTTGGCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCATTGCGTGCTGCTGTACCACCCGCCGCTGTCG  >  NC_002947/6058316‑6058609
                                                                                                                                                 |                                                                                                                                                    
aGCACCTATGAGGAGCTTGGCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGt                                                                                                                                                   >  3:127125/1‑149 (MQ=255)
  cACCTATGAGGAGCTTGGCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCCTCAGCTTCTGGTCGCCACCACGGTCATCGAGGTgg                                                                                                                                                 <  6:90284/149‑1 (MQ=255)
     cTATGAGGAGCTTGGCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCg                                                                                                                                              >  4:269303/1‑149 (MQ=255)
      tATGAGGAGCTTGGCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCCCCACGGGCAACGCGGTGGGCGt                                                                                                                                             >  8:307704/1‑149 (MQ=255)
          aggagCTTGGCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGa                                                                                                                                          <  1:28796/148‑1 (MQ=255)
           ggcgcctGGCAGTGCGCGGGGTGTAATGAAGGGGGGGTTTATCCATGGGCGCATGAAACCGGCATAAAAAGCCGATATCATGGCCGAGTTCAAGGCTGGCAATCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGAc                                                                                                                                         <  8:57348/142‑1 (MQ=255)
           ggagCTTGGCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACg                                                                                                                                        >  3:366288/1‑149 (MQ=255)
             agCTTGGCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTg                                                                                                                                            >  6:103448/1‑143 (MQ=255)
             agCTTGGCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTg                                                                                                                                            <  5:103448/143‑1 (MQ=255)
              gCTTGGCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGc                                                                                                                                     >  1:7982/1‑149 (MQ=255)
               cTTGGCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGcc                                                                                                                                    <  8:52370/149‑1 (MQ=255)
                  ggCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGccc                                                                                                                                   <  2:257912/147‑1 (MQ=255)
                     aGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGc                                                                                                                              >  6:76671/1‑149 (MQ=255)
                         cgcTGGGTGAACTGAAGGTGGGGTTGATCCTTGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGc                                                                                                                          >  1:31092/1‑149 (MQ=255)
                          gcTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGTCGTGGACGTGCCCAACGCCAGc                                                                                                                          >  3:49088/1‑148 (MQ=255)
                          gcTGGGTGAACTGAAGGTGGGGTAGATCCATGGGCGCAAGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGcc                                                                                                                         <  8:309576/149‑1 (MQ=255)
                                    cTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGatcat                                                                                                               >  3:66647/1‑149 (MQ=255)
                                     tGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGatcatc                                                                                                              >  1:115387/1‑149 (MQ=255)
                                       aaGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGatcatc                                                                                                              >  5:42236/1‑147 (MQ=255)
                                       aaGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGatcatc                                                                                                              <  2:115387/147‑1 (MQ=255)
                                       aaGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCga                                                                                                            <  2:7982/149‑1 (MQ=255)
                                       aaGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCg                                                                                                             <  3:269303/148‑1 (MQ=255)
                                         ggTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGCCGTCGACCTGCCCAACCCCAGCCTGATGATCATCgag                                                                                                           >  3:294197/1‑148 (MQ=255)
                                         cgtcgggTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACAACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCgag                                                                                                           >  2:19131/5‑148 (MQ=255)
                                            ggggTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGcc                                                                                                                         >  6:124872/1‑131 (MQ=255)
                                               gTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCg                                                                                                    <  7:287165/149‑1 (MQ=255)
                                                ttGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGa                                                                                                   <  2:155659/149‑1 (MQ=255)
                                                  gATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAgc                                                                                                 >  8:216773/1‑149 (MQ=255)
                                                    tCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAgcgc                                                                                               <  4:368410/149‑1 (MQ=255)
                                                     ccATGGGCGCATGAATCCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAgcgc                                                                                               <  6:139535/148‑1 (MQ=255)
                                                     ccATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCtg                                                                                                                       >  2:110202/1‑124 (MQ=255)
                                                     ccATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCtg                                                                                                                       <  1:110202/124‑1 (MQ=255)
                                                      cATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCt                                                                                             >  2:58590/1‑149 (MQ=255)
                                                       aTGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTg                                                                                            <  8:49193/149‑1 (MQ=255)
                                                        tGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAAc                                                                                                        >  6:321048/1‑136 (MQ=255)
                                                        tGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAAc                                                                                                        <  5:321048/136‑1 (MQ=255)
                                                        tGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTgg                                                                                           <  7:216773/149‑1 (MQ=255)
                                                          ggCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTggg                                                                                          <  4:200296/148‑1 (MQ=255)
                                                             gcATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCtt                                                                                      >  6:186799/1‑149 (MQ=255)
                                                                    aCCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAAc                                                                                                                                <  2:272634/100‑1 (MQ=255)
                                                                    aCCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAAc                                                                                                                                >  1:272634/1‑100 (MQ=255)
                                                                     ccGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGtt                                                                              <  5:76671/149‑1 (MQ=255)
                                                                     ccGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGtt                                                                              <  1:19131/149‑1 (MQ=255)
                                                                       ggCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGc                                                                            <  2:31092/149‑1 (MQ=255)
                                                                       ggCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACACCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGc                                                                            <  6:125849/149‑1 (MQ=255)
                                                                         cAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCAc                                                                          <  4:66647/149‑1 (MQ=255)
                                                                              aaaGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCt                                                                     >  6:195684/1‑149 (MQ=255)
                                                                                  acgtgaTCATGGGCGAGTTCAAGGCTGGCTACCTGCAGCTTCTGGTCGCCACCACTGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGt                                                                 <  4:49088/145‑1 (MQ=255)
                                                                                   cGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGTGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTg                                                                >  1:316456/1‑149 (MQ=255)
                                                                                   cGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTg                                                                <  6:374195/149‑1 (MQ=255)
                                                                                                   ttCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCgggg                                                <  5:186799/149‑1 (MQ=255)
                                                                                                   ttCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCgggg                                                <  8:166563/149‑1 (MQ=255)
                                                                                                          caggCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGcc                                         <  7:307704/147‑1 (MQ=255)
                                                                                                          cTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGcc                                         <  5:195684/149‑1 (MQ=255)
                                                                                                           tGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCg                                        <  4:366288/149‑1 (MQ=255)
                                                                                                            ggCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCtgatga                                                                                                                   <  1:301819/73‑1 (MQ=255)
                                                                                                            ggCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCtgatga                                                                                                                   >  2:301819/1‑73 (MQ=255)
                                                                                                              cAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGa                                     <  1:58590/149‑1 (MQ=255)
                                                                                                                  cTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCa                                 >  1:181582/1‑149 (MQ=255)
                                                                                                                  cTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCa                                 >  8:215773/1‑149 (MQ=255)
                                                                                                                   tGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTCGGCCGGGGCAGCGCCGTGAGCCAt                                >  5:120696/1‑149 (MQ=255)
                                                                                                                    gCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGtt                                                        <  1:137848/124‑1 (MQ=255)
                                                                                                                    gCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGtt                                                        >  2:137848/1‑124 (MQ=255)
                                                                                                                        cTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCATTGCGt                           <  6:42236/149‑1 (MQ=255)
                                                                                                                             cacgcccacctCGATGACCGTGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGc                                                                      >  4:212907/6‑101 (MQ=38)
                                                                                                                             cacgcccacctCGATGACCGTGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGc                                                                      <  3:212907/96‑1 (MQ=38)
                                                                                                                                gcccacctCGATGACCGTGGTGGGCGTGGACGTGCCCAACGCCAGCCtgat                                                                                                                     <  3:331670/49‑1 (MQ=18)
                                                                                                                                gcccacctCGATGACCGTGGTGGGCGTGGACGTGCCCAACGCCAGCCtgat                                                                                                                     >  4:331670/3‑51 (MQ=18)
                                                                                                                                cGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCATTGCGTGctgctg                    >  7:181462/1‑148 (MQ=255)
                                                                                                                                   caccacGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCATTGCGTGCTGCTGTac                 >  4:38385/1‑148 (MQ=255)
                                                                                                                                   caccacGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCATTGCGTGCTGCTGTac                 >  2:362261/1‑148 (MQ=255)
                                                                                                                                     ccacGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCATTGCGTGCTGCTGTaccac              <  8:181462/149‑1 (MQ=255)
                                                                                                                                         ggTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCATTGCGTGCTGCTGTACCACccgc          >  6:297978/1‑149 (MQ=255)
                                                                                                                                         ggTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCATTGCGTGCTGCTGTACCACccgc          >  5:332464/1‑149 (MQ=255)
                                                                                                                                           taaTCGAGGTGTGCGTGGAAGTGCCCAATGCCAGCCTGATGATCATCGAGAACCGCGAGCGCGTGGGCCTTGCGCAGTTGCACCTGCTGGGTGGCAGGGTTGGCCGGGGCAGGGCCGTGAGCCATTGCGTGCTGCTGTACCACccgccg        <  4:294197/147‑1 (MQ=255)
                                                                                                                                           tCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAAGCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCATTGCGTGCTGCTGTACCACccgccg        <  3:37843/149‑1 (MQ=255)
                                                                                                                                           tCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCATTGCGTGCTGCTGTACCACccgccg        <  5:297978/149‑1 (MQ=255)
                                                                                                                                                 aGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCATTGCGTGCTGCTGTACCACCCGCCGCTGTCg  <  6:332464/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AGCACCTATGAGGAGCTTGGCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCATTGCGTGCTGCTGTACCACCCGCCGCTGTCG  >  NC_002947/6058316‑6058609

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: