Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,058,454 G→A 5.2% V551I (GTC→ATC)  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,4540GA 5.2% 239.1 / 2.0 77V551I (GTC→ATC) recGjunction‑specific ATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (35/38);  new base A (2/2);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

GCCGCCGAAAGCACCTATGAGGAGCTTGGCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCATTGCGTGCTGCTGTACCACCCGC  >  NC_002947/6058307‑6058601
                                                                                                                                                   |                                                                                                                                                   
gccgccGAAAGCACCTATGAGGAGCTTGGCAGTGCGCTGGGTGAACTGAAGGTGGCGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGt                                                                                                                                                    >  5:139535/1‑149 (MQ=255)
gccgccCAAAGCACCTATGAGGAGCTTGGCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGTCGCATGAAACCGGCAGAAAAAGCAGAGATCATGGCCGAGTACAAGGCAGGCAACCTGCAGCTTCTGGTCGCAACCACgg                                                                                                                                                     >  8:287165/1‑148 (MQ=255)
 ccgccgAAAGCACCTATGAGGAGCTTGGCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGt                                                                                                                                                    >  1:257912/1‑148 (MQ=255)
   gccgAAAGCACCTATGAGGAGCTTGGCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCAt                                                                                                                                                 >  3:368410/1‑149 (MQ=255)
   cccgcacgCACCTCTGACGAGCTTGGCAGTGCGCTGGGTGAACTGCAGGTGGGGCGGATCCATGGGCGCATGAAACCGGCACAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCAt                                                                                                                                                 <  4:226620/142‑1 (MQ=255)
      gAAAGCACCTATGAGGAGCTTGGCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGAAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATc                                                                                                                                                >  7:166563/1‑147 (MQ=255)
         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                                                                                                               <  1:110202/124‑1 (MQ=255)
                                                              ccATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCtg                                                                                                               >  2:110202/1‑124 (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                                                                                                                        >  1:272634/1‑100 (MQ=255)
                                                                             aCCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAAc                                                                                                                        <  2:272634/100‑1 (MQ=255)
                                                                              ccGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGtt                                                                      <  1:19131/149‑1 (MQ=255)
                                                                              ccGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGtt                                                                      <  5:76671/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                                                                                                           >  2:301819/1‑73 (MQ=255)
                                                                                                                     ggCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCtgatga                                                                                                           <  1:301819/73‑1 (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                                                                                                             >  4:331670/3‑51 (MQ=18)
                                                                                                                                         gcccacctCGATGACCGTGGTGGGCGTGGACGTGCCCAACGCCAGCCtgat                                                                                                             <  3:331670/49‑1 (MQ=18)
                                                                                                                                         cGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCATTGCGTGctgctg            >  7:181462/1‑148 (MQ=255)
                                                                                                                                            caccacGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCATTGCGTGCTGCTGTac         >  2:362261/1‑148 (MQ=255)
                                                                                                                                            caccacGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCATTGCGTGCTGCTGTac         >  4:38385/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)
                                                                                                                                                   |                                                                                                                                                   
GCCGCCGAAAGCACCTATGAGGAGCTTGGCAGTGCGCTGGGTGAACTGAAGGTGGGGTTGATCCATGGGCGCATGAAACCGGCAGAAAAAGCCGAGATCATGGCCGAGTTCAAGGCTGGCAACCTGCAGCTTCTGGTCGCCACCACGGTCATCGAGGTGGGCGTGGACGTGCCCAACGCCAGCCTGATGATCATCGAGAACCCCGAGCGCCTGGGCCTTGCGCAGTTGCACCAGCTGCGTGGCAGGGTTGGCCGGGGCAGCGCCGTGAGCCATTGCGTGCTGCTGTACCACCCGC  >  NC_002947/6058307‑6058601

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: