Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,953,953 A→T 6.3% L684H (CTT→CAT)  recB ← exonuclease V (RecBCD complex), beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,953,9530AT 6.3% 196.5 / 3.8 64L684H (CTT→CAT) recBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base T (2/2);  total (30/34)
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.76e-01

ACGCATTTGTTGGCTGGAGGCATTACTGTCTGGCTCGAGGATATGTTGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGTT  >  NC_000913/2953820‑2954100
                                                                                                                                     |                                                                                                                                                   
aCGCATTTGTTGGCTGGAGGCATTACTGTCTGGCTCGAGGATATGTTGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCg                                                                                                                                       >  7:65173/1‑148 (MQ=255)
   cATTTGTTGGCTGGAGGCATTACTGTCTGGCTCGAGGATATGTTGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGAcgc                                                                                                                                                >  5:278538/1‑136 (MQ=255)
   cATTTGTTGGCTGGAGGCATTACTGTCTGGCTCGAGGATATGTTGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGAcgc                                                                                                                                                <  6:278538/136‑1 (MQ=255)
      ttgttgGCTGGAGGCATTACTGTCTGGCTCGAGGATATGTTGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGcc                                                                                                                                 >  7:146191/1‑148 (MQ=255)
        gttgGCTGGAGGCATTACTGTCTGGCTCGAGGATATGTTGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCg                                                                                                                                <  8:146191/147‑1 (MQ=255)
           ggCTGGAGGCATTACTGTCTGGCTCGAGGATATGTTGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTg                                                                                                                             <  1:366389/147‑1 (MQ=255)
             cTGGAGGCATTACTGGCTGGCTCGAGGATATGTTGCGATAACCAGCGTACCACCGCAGGGTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTACGACGCCGCTCACCGCCTGCCGTTGCCa                                                                                                                          <  4:182453/148‑1 (MQ=255)
                gAGGCATTACTGTCTGGCTCGAGGATATGTTGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGc                                                                                                                        >  7:39869/1‑147 (MQ=255)
                    cATTACTGTCTGGCTCGAGGATATGTTGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGtt                                                                                                                   <  4:221315/148‑1 (MQ=255)
                     aTTACTGTCTGGCTCGAGGATATGTTGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCCGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGttt                                                                                                                  <  1:140608/148‑1 (MQ=255)
                       tACTGTCTGGCTCGAGGATATGTTGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGtttt                                                                                                                 >  7:407446/1‑147 (MQ=255)
                       tACTGTCTGGCTCGACGATATGTTGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGtttt                                                                                                                 >  7:329047/1‑147 (MQ=255)
                        aCTGTCTGGCTCGAGGATATGTTGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGtttt                                                                                                                 >  2:257045/1‑146 (MQ=255)
                        aCTGTCTGGCTCGAGGATATGTTGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCATTGCCAGCAAGtttt                                                                                                                 >  3:72739/1‑146 (MQ=255)
                           gTCTGGCTCGAGGATATGTTGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGc                                                                                                             >  4:370077/1‑147 (MQ=255)
                            tCTGGCTCGAGGATATGTTGCGATAACCAGCGTACCAGCGGATGGTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGTTTTGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCaa                                                                                                           <  8:127047/148‑1 (MQ=255)
                             cTGGCTCGAGGATATGTTGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAAt                                                                                                          >  1:35843/1‑148 (MQ=255)
                                gCTCGAGGATATGTTGCGATACCCAGCGGACCAGCGCATGTTCACTTTCCAGCTGCGGTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTGCAGCATTGt                                                                                                        <  4:281632/147‑1 (MQ=255)
                                gCTCGAGGATATGTTGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCATTGCCAGCAAGTTTTCAGCAATGtt                                                                                                       <  4:72739/148‑1 (MQ=255)
                                              tGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACAt                                                                                           >  7:337147/1‑146 (MQ=255)
                                                  aTAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAgcgc                                                                                     >  3:333417/1‑148 (MQ=255)
                                                       cAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTa                                                                                                      >  3:221371/1‑126 (MQ=255)
                                                       cAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTa                                                                                                      <  4:221371/126‑1 (MQ=255)
                                                        aGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCgcag                                                                               >  3:42367/1‑148 (MQ=255)
                                                              ccAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCgg                                                                         >  2:367603/1‑148 (MQ=255)
                                                               cAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGtt                                                                                                       >  7:259576/1‑117 (MQ=255)
                                                               cAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGtt                                                                                                       <  8:259576/117‑1 (MQ=255)
                                                               cAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGc                                                                        <  8:407446/148‑1 (MQ=255)
                                                                 gcgcATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCa                                                                       >  1:55096/1‑147 (MQ=255)
                                                                      tGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGc                                                                 <  8:71660/148‑1 (MQ=255)
                                                                       gTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTATGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGc                                                                 >  3:231524/1‑147 (MQ=255)
                                                                       gTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTATGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGc                                                                 >  1:335063/1‑147 (MQ=255)
                                                                       gTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTAcgcg                                                                                                  <  3:91044/114‑1 (MQ=255)
                                                                       gTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTAcgcg                                                                                                  >  4:91044/1‑114 (MQ=255)
                                                                       gTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGc                                                                 >  7:299169/1‑147 (MQ=255)
                                                                        ttCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCTTCGGCATAACGcc                                                                <  3:370077/147‑1 (MQ=255)
                                                                        ttCACTATCCAGCTGCGTTCCGGCTTATTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGGTGCCAGCTAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCTGCTTTACGCCa                                                               <  1:341955/148‑1 (MQ=255)
                                                                          cACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACg                                                             >  6:257153/1‑148 (MQ=255)
                                                                           aCTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACg                                                             <  8:180942/147‑1 (MQ=255)
                                                                            cTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTATGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGtt                                                           <  2:335063/148‑1 (MQ=255)
                                                                            cTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTATGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGtt                                                           <  4:231524/148‑1 (MQ=255)
                                                                              ttCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTg                                                         >  3:234099/1‑148 (MQ=255)
                                                                               tCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGttt                                                          <  8:39869/146‑1 (MQ=255)
                                                                                  aGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTgcgc                                                      >  5:303022/1‑147 (MQ=255)
                                                                                       cGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCGCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGAtt                                                <  8:329047/148‑1 (MQ=255)
                                                                                          tCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGc                                             >  6:348993/1‑148 (MQ=255)
                                                                                           ccGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGcc                                            >  2:384374/1‑148 (MQ=255)
                                                                                           ccGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGc                                             <  2:35843/147‑1 (MQ=255)
                                                                                            cGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATc                                                                                          >  8:3345/1‑101 (MQ=255)
                                                                                            cGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATc                                                                                          <  7:3345/101‑1 (MQ=255)
                                                                                            cGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCg                                           <  7:14760/148‑1 (MQ=255)
                                                                                                        gCAGTTCGCTGATCTGCAAGATATCTGTAAGACGCCGCTCACCTCCTGCCTGTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATcta                                <  1:53739/147‑3 (MQ=255)
                                                                                                        gCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGa                                <  5:257153/147‑1 (MQ=255)
                                                                                                        gCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGa                                <  4:42367/147‑1 (MQ=255)
                                                                                                        gCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGa                                <  4:333417/147‑1 (MQ=255)
                                                                                                         cAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGaa                               >  6:333003/1‑147 (MQ=255)
                                                                                                              cGCTGATATGCAAGATCTCGGTAAGACGCCGCTCAGCGCGTGCCGTTGCCAGCAAGTTTTCAGCTATGTTACGCGCCGACATCATCGCCCGCAGCATCGGCATTACGCCACGTTTGCGACAGATTTGCCGATAACCATCGAAGTCtt                          <  8:337147/147‑1 (MQ=255)
                                                                                                               gCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCg                        <  1:257045/148‑1 (MQ=255)
                                                                                                                  gATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGAc                      >  1:391733/1‑147 (MQ=255)
                                                                                                                  gATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGAc                      >  5:314221/1‑147 (MQ=255)
                                                                                                                   atatGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGAc                      <  2:55096/146‑1 (MQ=255)
                                                                                                                             tatCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATccc           <  3:34615/147‑1 (MQ=255)
                                                                                                                             tatCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCa          <  5:348993/148‑1 (MQ=255)
                                                                                                                               tCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACg        <  1:384374/148‑1 (MQ=255)
                                                                                                                                     aGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGt   <  2:391733/147‑1 (MQ=255)
                                                                                                                                     aGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATCtt  >  1:163013/1‑148 (MQ=255)
                                                                                                                                     |                                                                                                                                                   
ACGCATTTGTTGGCTGGAGGCATTACTGTCTGGCTCGAGGATATGTTGCGATAACCAGCGTACCAGCGCATGTTCACTTTCCAGCTGCGTTCCGGCTTCTTGTAGCAGTTCGCTGATATGCAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGTT  >  NC_000913/2953820‑2954100

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: