Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,133,979 C→G 9.2% A351P (GCA→CCA) ‡ PP_3638 ← acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,133,9790CG 9.2% 268.7 / 14.4 87A351P (GCA→CCA) ‡PP_3638acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (37/42);  new base G (6/2);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

CTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAG  >  NC_002947/4133838‑4134117
                                                                                                                                               |                                                                                                                                          
cTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGt                                                                                                                                     <  4:156962/149‑1 (MQ=255)
cTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGt                                                                                                                                     <  4:151385/149‑1 (MQ=255)
   gggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGtt                                                                                                                                    <  6:231143/147‑1 (MQ=255)
   gggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTg                                                                                                                                   <  8:102555/148‑1 (MQ=255)
   gggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGc                                                                                                                                  <  5:147325/149‑1 (MQ=255)
     gCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTg                                                                                                                                   >  4:2570/1‑146 (MQ=255)
     gCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTg                                                                                                                                   <  3:2570/146‑1 (MQ=255)
      ccAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGt                                                                                                                               >  7:350261/1‑149 (MQ=255)
      ccAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGt                                                                                                                               >  2:96139/1‑149 (MQ=255)
         aGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGc                                                                                                                                  >  4:79704/1‑143 (MQ=255)
         aGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGc                                                                                                                                  <  3:79704/143‑1 (MQ=255)
         aGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCg                                                                                                                            <  8:334271/149‑1 (MQ=255)
                tCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGgcgcgc                                                                                                                      >  2:365433/1‑148 (MQ=255)
                     aCCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCCCGCTCA‑‑TGCAACGCCTGCCTCGCGTTGCAGTTCGGCGCGCAAttcc                                                                                                                <  4:7324/149‑3 (MQ=255)
                       cTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGTAGTTCGGCGCGCAATTg                                                                                                                 >  6:340406/1‑146 (MQ=255)
                       cTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGTAGTTCGGCGCGCAATTg                                                                                                                 <  5:340406/146‑1 (MQ=255)
                        tCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCAcc                                                                                                             >  6:69861/1‑149 (MQ=255)
                        tCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCAc                                                                                                              >  5:109162/1‑148 (MQ=255)
                          aCCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGACCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAg                                                                                                           >  4:45746/1‑149 (MQ=255)
                           ccTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAgg                                                                                                          <  7:46084/149‑1 (MQ=255)
                            cTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCAGGCAGGCGTTGc                                                                                                                                  >  6:37234/1‑124 (MQ=255)
                            cTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCAGGCAGGCGTTGc                                                                                                                                  <  5:37234/124‑1 (MQ=255)
                               gTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTgg                                                                                                      <  6:129816/149‑1 (MQ=255)
                                     gCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                >  4:288537/1‑149 (MQ=255)
                                         cAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGa                                                                                             <  2:75594/148‑1 (MQ=255)
                                          aGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATc                                                                                           <  2:277092/149‑1 (MQ=255)
                                            gTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCgg                                                                                         <  6:109162/149‑1 (MQ=255)
                                               gAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAgg                                                                                                          <  3:3606/129‑1 (MQ=255)
                                               gAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAgg                                                                                                          >  4:3606/1‑129 (MQ=255)
                                                aaCTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCAc                                                                                      >  5:232049/1‑148 (MQ=255)
                                                aaCTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCAc                                                                                      >  7:71766/1‑148 (MQ=255)
                                                      gcCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                <  5:145944/132‑1 (MQ=255)
                                                      gcCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                >  6:145944/1‑132 (MQ=255)
                                                         aGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGAc                                                                            <  3:288537/149‑1 (MQ=255)
                                                               tgctgcCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTgg                                                                                                      <  5:356482/117‑1 (MQ=255)
                                                               tgctgcCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTgg                                                                                                      >  6:356482/1‑117 (MQ=255)
                                                                  tgcCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCAcc                                                                                                             >  8:167011/1‑107 (MQ=255)
                                                                  tgcCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCAcc                                                                                                             <  7:167011/107‑1 (MQ=255)
                                                                    ccAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcg                                                                 >  1:310069/1‑149 (MQ=255)
                                                                      aaCCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggc                                                               >  2:249085/1‑149 (MQ=255)
                                                                       accaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcg                                                              <  1:88964/149‑1 (MQ=255)
                                                                       accaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcg                                                              <  2:264600/149‑1 (MQ=255)
                                                                       accaccCTCACGCGGGTAGCCAAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGATGGCCGACTCCCGCCATcggcc                                                              >  6:85999/1‑148 (MQ=255)
                                                                        ccaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                             <  8:130626/149‑1 (MQ=255)
                                                                        ccaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                             <  6:232049/149‑1 (MQ=255)
                                                                        ccaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                             <  8:125234/149‑1 (MQ=255)
                                                                         caccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCAAGTGCAACGCATGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCAATcggcg                                                              >  1:293698/1‑149 (MQ=255)
                                                                          accCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGc                                                           <  8:350261/149‑1 (MQ=255)
                                                                           cccTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                >  2:47309/1‑111 (MQ=255)
                                                                           cccTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                <  1:47309/111‑1 (MQ=255)
                                                                                aCGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAccc                                                     >  5:115306/1‑149 (MQ=255)
                                                                                aCGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAccc                                                     <  1:365433/149‑1 (MQ=255)
                                                                                aCGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAc                                                       <  5:260801/147‑1 (MQ=255)
                                                                                      gTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGcc                                               >  4:247677/1‑149 (MQ=255)
                                                                                         gCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTa                                            >  2:155607/1‑149 (MQ=255)
                                                                                           ccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGc                                                                                                         <  5:367897/86‑1 (MQ=255)
                                                                                           ccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGc                                                                                                         >  6:367897/1‑86 (MQ=255)
                                                                                           ccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTAtt                                          >  3:169542/1‑149 (MQ=255)
                                                                                                       cACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGcc                                               <  7:42597/132‑1 (MQ=255)
                                                                                                       cACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGcc                                               >  8:42597/1‑132 (MQ=255)
                                                                                                        aCGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAc                                                       >  5:8210/1‑123 (MQ=255)
                                                                                                        aCGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAc                                                       <  6:8210/123‑1 (MQ=255)
                                                                                                              cTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTAt                                           <  5:283258/129‑1 (MQ=255)
                                                                                                              cTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTAt                                           >  6:283258/1‑129 (MQ=255)
                                                                                                              cTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgcc                       <  2:293698/149‑1 (MQ=255)
                                                                                                                 cAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgccgc                     >  1:372880/1‑148 (MQ=255)
                                                                                                                 cAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgccgc                     <  4:169542/148‑1 (MQ=255)
                                                                                                                 cAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgccg                      <  2:310069/147‑1 (MQ=255)
                                                                                                                 cAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGAGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgccgc                     >  8:215329/1‑148 (MQ=255)
                                                                                                                   gCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCAc                  <  1:155607/149‑1 (MQ=255)
                                                                                                                           gcgcTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGc          <  6:115306/149‑1 (MQ=255)
                                                                                                                              ccgaactGCAACGCCAGGCAGGCGTTGCA‑‑TGAGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTg                >  8:126061/7‑140 (MQ=37)
                                                                                                                              ccgaactGCAACGCCAGGCAGGCGTTGCA‑‑TGAGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTg                <  7:126061/134‑1 (MQ=37)
                                                                                                                               tCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGAGCAATTGCACCAGGCTGGGCGTGAAGATCGGCACGAAGGCCGACTCCAGCCATCGGCGGGCAAACCCTTCCCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTCGCCTGCAACt      >  6:326020/1‑149 (MQ=255)
                                                                                                                               ggatctGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCaa        <  5:107969/144‑1 (MQ=255)
                                                                                                                                 a‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGc                                                           <  7:190400/94‑1 (MQ=255)
                                                                                                                                 a‑‑TGCAACGCCTGCCTGGCGTTCCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGAACGAAGGCCGACTCCCGCCAACGGCggga                                                           >  8:190400/1‑93 (MQ=255)
                                                                                                                                  actGCAACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc     >  2:42243/3‑149 (MQ=255)
                                                                                                                                    tGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCATGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCa   <  3:45746/149‑1 (MQ=255)
                                                                                                                                    tGCAACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGTGTGACGATCGGCACGAAGACCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc     >  1:352840/1‑147 (MQ=255)
                                                                                                                                    tGCAACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGTGTGACGATCGGCACGAAGACCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc     >  1:347169/1‑147 (MQ=255)
                                                                                                                                     gtaaCGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcg                                                              >  8:298883/3‑89 (MQ=255)
                                                                                                                                     gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCAc                  >  2:160829/1‑133 (MQ=255)
                                                                                                                                     gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCAc                  <  1:160829/133‑1 (MQ=255)
                                                                                                                                     gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAg  >  5:357925/1‑149 (MQ=255)
                                                                                                                                     gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAg  >  8:185992/1‑149 (MQ=255)
                                                                                                                                      cAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGCTAAGCCTTGCCGCCACTGGCCTGCaa        >  6:107969/1‑142 (MQ=255)
                                                                                                                                          gcctgcctGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc     <  8:121647/141‑1 (MQ=255)
                                                                                                                                          gcctgcctGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc     >  7:121647/1‑141 (MQ=255)
                                                                                                                                           cctgcctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTccc                                                                        >  5:77756/1‑73 (MQ=255)
                                                                                                                                           cctgcctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTccc                                                                        <  6:77756/73‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                          
CTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAG  >  NC_002947/4133838‑4134117

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑