Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,695,938 A→G 6.4% V77A (GTG→GCG)  ligD ← DNA ligase D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,695,9380AG 6.4% 221.2 / 4.4 78V77A (GTG→GCG) ligDDNA ligase D
Reads supporting (aligned to +/- strand):  ref base A (37/36);  new base G (3/2);  total (40/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 = 7.43e-01

CTTTGGCGTAGGCTTCATGCGCGTCGCCTTCAGGCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGCCCGTTT  >  NC_002947/3695798‑3696086
                                                                                                                                            |                                                                                                                                                      
cTTTGGCGTAGGCTTCATGCGCGTCGCCTTCAGGCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAgg                                                                                                                                                <  4:46818/149‑1 (MQ=255)
      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)
                                                                gTCCCCCGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGAACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATc                                                                                 >  4:363096/1‑147 (MQ=255)
                                                                 tCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAg                                                                              <  1:55187/149‑1 (MQ=255)
                                                                  ccccTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCGCCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGc                                                                             <  4:32005/149‑1 (MQ=255)
                                                                   cccTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCt                                                                            >  6:69341/1‑149 (MQ=255)
                                                                    ccTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTc                                                                           >  2:19822/1‑149 (MQ=255)
                                                                     cTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTc                                                                           >  5:392020/1‑148 (MQ=255)
                                                                        cACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAgg                                                                      >  2:115551/1‑149 (MQ=255)
                                                                              aGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGTGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAg                                                                >  8:50250/1‑149 (MQ=255)
                                                                                  gCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGc                                                                             >  4:330493/1‑133 (MQ=255)
                                                                                  gCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGc                                                                             <  3:330493/133‑1 (MQ=255)
                                                                                  gCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCg                                                             >  2:102095/1‑148 (MQ=255)
                                                                                   ccTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTa                                                           <  6:392020/149‑1 (MQ=255)
                                                                                    cTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTggg                                                                                                             >  3:123712/1‑100 (MQ=255)
                                                                                    cTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTggg                                                                                                             <  4:123712/100‑1 (MQ=255)
                                                                                     ttCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGt                                                         >  5:68338/1‑149 (MQ=255)
                                                                                       ccGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGt                                                         <  7:50250/147‑1 (MQ=255)
                                                                                               tGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACt                                               <  7:241518/149‑1 (MQ=255)
                                                                                                ggCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTggg                                                                                                             >  6:18348/1‑88 (MQ=255)
                                                                                                ggCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTggg                                                                                                             <  5:18348/88‑1 (MQ=255)
                                                                                                ggCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCa                                                                                            >  1:252493/1‑105 (MQ=255)
                                                                                                ggCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCa                                                                                            <  2:252493/105‑1 (MQ=255)
                                                                                                ggCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCg                                                                                <  7:225549/116‑1 (MQ=255)
                                                                                                ggCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCg                                                                                >  8:225549/1‑116 (MQ=255)
                                                                                                   ccTCGAAATCGGCGTAGTCCAGCGGGTGCTCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGcg                                           >  4:16059/1‑149 (MQ=255)
                                                                                                   ccTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCg                                                                    >  7:140808/1‑124 (MQ=255)
                                                                                                   ccTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCg                                                                    <  8:140808/124‑1 (MQ=255)
                                                                                                   ccTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGcg                                           >  5:321980/1‑149 (MQ=255)
                                                                                                      cGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGcgtcg                                        >  3:38768/1‑149 (MQ=255)
                                                                                                       gataTCGGCGTAGTCGAGCGGGTGATCCTCCACATGCACCGCCAGGCGGAGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCGAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGcgtcgt                                       <  2:196717/146‑1 (MQ=255)
                                                                                                              gCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTg                                <  1:18979/149‑1 (MQ=255)
                                                                                                              gCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTg                                <  1:25300/149‑1 (MQ=255)
                                                                                                               cGTAGGCCAGCGGGTGATCCTCCACTTGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGCCGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCTTCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCTTCGTGTTTCTgg                               <  5:212514/149‑1 (MQ=255)
                                                                                                                 tAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGAt                             <  6:191125/149‑1 (MQ=255)
                                                                                                                  aGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATg                            >  4:90326/1‑149 (MQ=255)
                                                                                                                            gTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCa                  <  4:38768/149‑1 (MQ=255)
                                                                                                                               aTCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAgcg               >  8:154203/1‑149 (MQ=255)
                                                                                                                               aTCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAgc                <  6:68338/148‑1 (MQ=255)
                                                                                                                               aTCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAgc                <  5:69341/148‑1 (MQ=255)
                                                                                                                                  ctccACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCAt            >  3:394953/1‑149 (MQ=255)
                                                                                                                                  ctccACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCACCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCAt            >  3:273595/1‑149 (MQ=255)
                                                                                                                                  ctccACATGAACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCa             <  3:363096/148‑1 (MQ=255)
                                                                                                                                  cccaaggtGCGCCGCCTGGCGGTGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCAt            >  4:398617/8‑149 (MQ=255)
                                                                                                                                     cacaTGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGACCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCAt            <  1:115551/147‑1 (MQ=255)
                                                                                                                                      acaTGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGc        >  5:291122/1‑149 (MQ=255)
                                                                                                                                       ggtGCGCCGCCTGGCGGTGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGcc       >  6:249240/3‑149 (MQ=255)
                                                                                                                                        gtGCGCCGCCTGGCGGTGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCAATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGAc                                                <  1:259439/107‑1 (MQ=255)
                                                                                                                                        gtGCGCCGCCTGGCGGTGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGAc                                                >  2:259439/2‑107 (MQ=255)
                                                                                                                                            aCCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGCCCGttt  <  6:48426/149‑1 (MQ=255)
                                                                                                                                            aCCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGCCCGttt  >  5:48426/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                      
CTTTGGCGTAGGCTTCATGCGCGTCGCCTTCAGGCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGCCCGTTT  >  NC_002947/3695798‑3696086

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: