Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,695,950 C→T 6.2% R73H (CGC→CAC)  ligD ← DNA ligase D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,695,9500CT 6.2% 243.1 / 3.6 81R73H (CGC→CAC) ligDDNA ligase D
Reads supporting (aligned to +/- strand):  ref base C (40/35);  new base T (3/2);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-01

CGTAGGCTTCATGCGCGTCGCCTTCAGGCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGCCCGTTTGC  >  NC_002947/3695804‑3696088
                                                                                                                                                  |                                                                                                                                            
cGTAGGCTTCATGCGCGTCGCCTTCAGGCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGgcggcg                                                                                                                                             >  2:55187/1‑148 (MQ=255)
 gTAGGCTTCATGCTCGTCGAATTCAGGCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCa                                                                                                                                           >  5:396237/1‑149 (MQ=255)
    ggCTTCATGCGCGTCGCCTTCAGGCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCGCCGCCTGGCGGTGCatg                                                                                                                                         <  5:11955/148‑3 (MQ=255)
    ggCTTCATGCGCGTCGCCTTCAGGCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCAcc                                                                                                                                         <  7:217868/148‑1 (MQ=255)
     gCTTCATGCGCGTCGCCTTCAGGCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATACACCGCCAGGCGGCGCACCtt                                                                                                                                       >  5:13686/1‑149 (MQ=255)
                    ccTTCAGGCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACgg                                                                                                                        <  3:293571/149‑1 (MQ=255)
                      ttCAGGCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCAAGCGCGTGATCCTCAACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGc                                                                                                                      >  1:196717/1‑149 (MQ=255)
                          ggCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCTCCAGGCGGCGAACCTTGGGGTCCAACGACGGGCCttt                                                                                                                  >  3:13840/1‑149 (MQ=255)
                             tCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACgg                                                                                                                        <  4:200880/140‑1 (MQ=255)
                             tCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACgg                                                                                                                        >  3:200880/1‑140 (MQ=255)
                               ccAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTggg                                                                                                               <  8:176877/147‑1 (MQ=255)
                               ccAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGa                                                                                                             >  7:114538/1‑149 (MQ=255)
                                 aGATACCTCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTAGGGGTCCAACGACGGGCCTTTGGGGAta                                                                                                           >  7:265258/1‑148 (MQ=255)
                                   aTACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGc                                                                                                         >  5:191125/1‑149 (MQ=255)
                                      gcgaggTCGCATACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATTGCCCTCTAAATCGGCTTAGTCCATCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGccc                                                                                                       <  8:265258/144‑1 (MQ=255)
                                       cgcgGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAg                                                                                                     >  8:241518/1‑149 (MQ=255)
                                       cgcgGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAg                                                                                                     >  1:243305/1‑149 (MQ=255)
                                        gcgGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCGCCTTGGGGTCCAACGACGGGCCGTTGGGGATCGCCCAGc                                                                                                    >  3:32005/1‑149 (MQ=255)
                                        gcgGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGc                                                                                                    <  2:243305/149‑1 (MQ=255)
                                        gcgGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAAACGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGCCGGCGCACCTTGGGGGCCAACGACGGGCCTTTGGGGATCGACCAGc                                                                                                    >  6:212514/1‑149 (MQ=255)
                                         cgGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATACACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCt                                                                                                   <  6:13686/149‑1 (MQ=255)
                                           gTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTg                                                                                                                                      >  5:298357/1‑112 (MQ=255)
                                           ctCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTg                                                                                                                                      <  6:298357/111‑1 (MQ=255)
                                              ccAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCtttt                                                                                                >  1:186373/1‑147 (MQ=255)
                                              ccAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTc                                                                                               <  6:396237/148‑1 (MQ=255)
                                              ccAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTc                                                                                               <  3:55379/148‑1 (MQ=255)
                                              ccAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCa                                                                                              <  2:186373/149‑1 (MQ=255)
                                                aGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAgg                                                                                            <  1:56280/149‑1 (MQ=255)
                                                aGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAgg                                                                                            <  8:114538/149‑1 (MQ=255)
                                                  aCGCTCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGGAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGt                                                                                          <  4:207129/149‑1 (MQ=255)
                                                          gTCCCCCGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGAACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATc                                                                                   >  4:363096/1‑147 (MQ=255)
                                                           tCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAg                                                                                <  1:55187/149‑1 (MQ=255)
                                                            ccccTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCGCCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGc                                                                               <  4:32005/149‑1 (MQ=255)
                                                             cccTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCt                                                                              >  6:69341/1‑149 (MQ=255)
                                                              ccTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTc                                                                             >  2:19822/1‑149 (MQ=255)
                                                               cTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTc                                                                             >  5:392020/1‑148 (MQ=255)
                                                                  cACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAg‑g                                                                        >  2:115551/1‑149 (MQ=255)
                                                                        aGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGTGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAg                                                                  >  8:50250/1‑149 (MQ=255)
                                                                            gCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGc                                                                               >  4:330493/1‑133 (MQ=255)
                                                                            gCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGc                                                                               <  3:330493/133‑1 (MQ=255)
                                                                            gCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCg                                                               >  2:102095/1‑148 (MQ=255)
                                                                             ccTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTa                                                             <  6:392020/149‑1 (MQ=255)
                                                                              cTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTggg                                                                                                               <  4:123712/100‑1 (MQ=255)
                                                                              cTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTggg                                                                                                               >  3:123712/1‑100 (MQ=255)
                                                                               ttCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGt                                                           >  5:68338/1‑149 (MQ=255)
                                                                                 ccGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGt                                                           <  7:50250/147‑1 (MQ=255)
                                                                                         tGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACt                                                 <  7:241518/149‑1 (MQ=255)
                                                                                          ggCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTggg                                                                                                               <  5:18348/88‑1 (MQ=255)
                                                                                          ggCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTggg                                                                                                               >  6:18348/1‑88 (MQ=255)
                                                                                          ggCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCa                                                                                              <  2:252493/105‑1 (MQ=255)
                                                                                          ggCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCa                                                                                              >  1:252493/1‑105 (MQ=255)
                                                                                          ggCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCg                                                                                  <  7:225549/116‑1 (MQ=255)
                                                                                          ggCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCg                                                                                  >  8:225549/1‑116 (MQ=255)
                                                                                             ccTCGAAATCGGCGTAGTCCAGCGGGTGCTCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGcg                                             >  4:16059/1‑149 (MQ=255)
                                                                                             ccTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCg                                                                      >  7:140808/1‑124 (MQ=255)
                                                                                             ccTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCg                                                                      <  8:140808/124‑1 (MQ=255)
                                                                                             ccTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGcg                                             >  5:321980/1‑149 (MQ=255)
                                                                                                cGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGcgtcg                                          >  3:38768/1‑149 (MQ=255)
                                                                                                 gataTCGGCGTAGTCGAGCGGGTGATCCTCCACATGCACCGCCAGGCGGAGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCGAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGcgtcgt                                         <  2:196717/146‑1 (MQ=255)
                                                                                                        gCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTg                                  <  1:25300/149‑1 (MQ=255)
                                                                                                        gCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTg                                  <  1:18979/149‑1 (MQ=255)
                                                                                                         cGTAGGCCAGCGGGTGATCCTCCACTTGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGCCGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑TTCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCTTCGTGTTTCTgg                                 <  5:212514/149‑1 (MQ=255)
                                                                                                           tAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGAt                               <  6:191125/149‑1 (MQ=255)
                                                                                                            aGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATg                              >  4:90326/1‑149 (MQ=255)
                                                                                                                      gTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCa                    <  4:38768/149‑1 (MQ=255)
                                                                                                                         aTCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAgcg                 >  8:154203/1‑149 (MQ=255)
                                                                                                                         aTCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAgc                  <  6:68338/148‑1 (MQ=255)
                                                                                                                         aTCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAgc                  <  5:69341/148‑1 (MQ=255)
                                                                                                                            ctccACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCAt              >  3:394953/1‑149 (MQ=255)
                                                                                                                            ctccACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCACCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCAt              >  3:273595/1‑149 (MQ=255)
                                                                                                                            ctccACATGAACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCa               <  3:363096/148‑1 (MQ=255)
                                                                                                                            cccaaggtGCGCCGCCTGGCGGTGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCAt              >  4:398617/8‑149 (MQ=255)
                                                                                                                               cacaTGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAGACCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCAt              <  1:115551/147‑1 (MQ=255)
                                                                                                                                acaTGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGc          >  5:291122/1‑149 (MQ=255)
                                                                                                                                 ggtGCGCCGCCTGGCGGTGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGcc         >  6:249240/3‑149 (MQ=255)
                                                                                                                                  gtGCGCCGCCTGGCGGTGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCAATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGAc                                                  <  1:259439/107‑1 (MQ=255)
                                                                                                                                  gtGCGCCGCCTGGCGGTGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGAc                                                  >  2:259439/2‑107 (MQ=255)
                                                                                                                                      aCCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGCCCGttt    >  5:48426/1‑149 (MQ=255)
                                                                                                                                      aCCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGCCCGttt    <  6:48426/149‑1 (MQ=255)
                                                                                                                                        cGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGAGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGCCCGTTTGc  >  5:30437/1‑149 (MQ=255)
                                                                                                                                        cGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTAAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGCCCGTTTGc  >  1:179923/1‑149 (MQ=255)
                                                                                                                                         gCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGCCCGTTTGc  >  4:333425/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
CGTAGGCTTCATGCGCGTCGCCTTCAGGCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCC‑ATCGAGCTCCAG‑GCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGCCCGTTTGC  >  NC_002947/3695804‑3696088

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