Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,082,380 G→A 5.8% T629I (ACC→ATC)  alaS ← alanine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,082,3800GA 5.8% 206.3 / 2.4 69T629I (ACC→ATC) alaSalanine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (33/32);  new base A (2/2);  total (35/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 = 6.06e-01

CAGGCTGATGTCACCGGTGCGCTTGGCGTGGATACCGCCGCACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCGAAGTGGCTGAAGTCGAAACGCAGGCGCTGGCTGTCGACCAGCGAGCCCTTCTGCTGT  >  NC_002947/5082232‑5082526
                                                                                                                                                    |                                                                                                                                                  
cAGGCTGATGTCACCGGTGCGCTTGGCGTGGATACCGCCGCACAATTCCACCGAAAAATCACCGCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGt                                                                                                                                                   >  7:73026/1‑149 (MQ=255)
   gCTGATGTCACCGGTGCGCTTGGCGTGGATACCGCCGCACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGt                                                                                                                                                   >  3:290791/1‑147 (MQ=255)
   gCTGATGTCACCGGTGCGCTTGGCGTGGATACCGCCGCACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGt                                                                                                                                                   >  4:292972/1‑147 (MQ=255)
   gCTGATGTCACCGGTGCGCTTGGCGTGGATACCGCCGCACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGt                                                                                                                                                   >  4:288194/1‑147 (MQ=255)
       aTGTCACCGGTGCGCTTGGCGTGGATACCGCCGCACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTc                                                                                                                                             <  4:290791/149‑1 (MQ=255)
          tCACCGGTGCGCTTGGCGTGGATACCGCCGCACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGt                                                                                                                                          <  1:63793/149‑1 (MQ=255)
             ccGGTGCGCTTGGCGTGGATACCGCCGCACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTc                                                                                                                                       <  7:211987/149‑1 (MQ=255)
             ccGGTGCGCTTGGCGGGGATACCGCCGCACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTc                                                                                                                                       <  4:282919/149‑1 (MQ=255)
                gTGCGCTTGGCGTGGATACCGCCGCACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAAc                                                                                                                                    >  7:294638/1‑149 (MQ=255)
                    gcTTGGCGTGGATACCTCCGCACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGGTTCGATGTCGGTCAGTTCGGTTTCAACCggg                                                                                                                                >  6:77795/1‑149 (MQ=255)
                      ttGGCGTGGATACCGCCGCACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCgggg                                                                                                                               >  8:190213/1‑148 (MQ=255)
                          cGTGGATACCGCCGCACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTc                                                                                                                          >  7:70209/1‑149 (MQ=255)
                            tGGATACCGCCGCACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTAtt                                                                                                                           >  5:261389/1‑146 (MQ=255)
                            tGGATACCGCCGCACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTAtt                                                                                                                           <  6:261389/146‑1 (MQ=255)
                                 aCCGCCGCACAATTCCCCCGAAAAATCACCGCCCATGCTCAGCACCCGCACGGTGTCGCCGTATTTCTCGCCGAACCTCGCCATTGCGCCTGTGGCCTTGGAAGTTTCGATGTCGGTCAGTTCGGTTTCAAACGGGGTATTCTTGCTaa                                                                                                                   >  3:8883/1‑149 (MQ=255)
                                   cgccgcACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCGTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGCAtc                                                                                                                 >  4:87341/1‑149 (MQ=255)
                                    gccgcACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGtt                                                                                                                                         <  4:249890/124‑1 (MQ=255)
                                    gccgcACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGtt                                                                                                                                         >  3:249890/1‑124 (MQ=255)
                                    gccgcACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAAtct                                                                                                                >  6:5411/1‑149 (MQ=255)
                                     ccgcACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGAGGTCGGTCAGTTCGGTTTCAGCCGGGGTATTCTTGCGAAtctc                                                                                                               >  8:146837/1‑149 (MQ=255)
                                       gcACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTcgc                                                                                                             <  7:16871/149‑1 (MQ=255)
                                                cACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCggg                                                                                                                                <  1:142416/121‑1 (MQ=255)
                                                cACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCggg                                                                                                                                >  2:142416/1‑121 (MQ=255)
                                                cACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGa                                                                                                    <  8:294638/149‑1 (MQ=255)
                                                   cGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGt                                                                                                 >  7:272962/1‑149 (MQ=255)
                                                    gAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTc                                                                                                >  2:23571/1‑149 (MQ=255)
                                                    gAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGATCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTc                                                                                                >  6:219414/1‑149 (MQ=255)
                                                    gAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGATCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTc                                                                                                >  6:214747/1‑149 (MQ=255)
                                                     aaaaaTCACCGCCCATTATCAGCACACGCACTGTGTCGCCGTATTTCTCGACGAACAGCGCAAATGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGTTTTCAACAGGCGTATTCTTTCGAATCTCGCTGTTGCCGATGTCt                                                                                               >  4:176039/1‑149 (MQ=255)
                                                     aaaaaTCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCggg                                                                                                                                <  7:45348/116‑1 (MQ=255)
                                                     aaaaaTCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCggg                                                                                                                                >  8:45348/1‑116 (MQ=255)
                                                      aaaaTCACCGCCCATGCTCAGAACACGCACTGTGTCGCCGTATTTCTCGCCGAACAGCGGCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCtt                                                                                              <  3:176039/149‑1 (MQ=255)
                                                        aaTCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTcc                                                                                            <  8:272962/149‑1 (MQ=255)
                                                        aaTCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTcc                                                                                            <  3:122460/149‑1 (MQ=255)
                                                              ccgccATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAgctgc                                                                                      <  8:73026/146‑1 (MQ=255)
                                                                cccATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCtt                                                                                                                        >  8:85453/1‑113 (MQ=255)
                                                                cccATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCtt                                                                                                                        <  7:85453/113‑1 (MQ=255)
                                                                         aGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAgctgc                                                                                      >  5:33483/1‑138 (MQ=255)
                                                                         aGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAgctgc                                                                                      <  6:33483/138‑1 (MQ=255)
                                                                           cacaCGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTAtt                                                                                                                           <  7:113841/99‑1 (MQ=255)
                                                                           cacaCGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTAtt                                                                                                                           >  8:113841/1‑99 (MQ=255)
                                                                             cacgcacgTTGTCGCCGTATTTCTCGCCGAACCGCGCCATTGCGCCTTTTGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGTTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGc                                                                       <  4:8883/149‑1 (MQ=255)
                                                                             cacgcacgGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGc                                                                       <  5:5411/149‑1 (MQ=255)
                                                                              acgcacgGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGACTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCt                                                                      >  1:61765/1‑149 (MQ=255)
                                                                                  acgGTGTCGCCGTATTTCTAGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTGGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATGTCGCGGTTGAGGATGTCTTGCAGGTGCTTGATGTGCTCCGGCTGCAc                                                                  >  8:43742/1‑149 (MQ=255)
                                                                                          gCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCg                                                           >  8:178602/1‑148 (MQ=255)
                                                                                           ccGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTc                                                                                             <  5:305784/113‑1 (MQ=255)
                                                                                           ccGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTc                                                                                             >  6:305784/1‑113 (MQ=255)
                                                                                           ccGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCg                                                           >  2:256613/1‑147 (MQ=255)
                                                                                             gTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCGAAGt                                                       <  8:232406/149‑1 (MQ=255)
                                                                                                    tcGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCGAAGTgg                                                     <  3:15856/144‑1 (MQ=255)
                                                                                                    tcGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCGAAGTgg                                                     >  4:15856/1‑144 (MQ=255)
                                                                                                      gCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCGAAGTGGCTGAAGt                                              >  4:58711/1‑149 (MQ=255)
                                                                                                           aCAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTGGGTCAGTTCGGTTTCAACCGGGGTTTTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCGAAGTGGCTGAAGTCGaaa                                         <  7:43742/149‑1 (MQ=255)
                                                                                                                       cgcCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCGAAGTGGCTGAAGTCGAAACGCAGGCgctgg                             >  5:97195/1‑149 (MQ=255)
                                                                                                                       cgcCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCGAAGTGGCTGAAGTCGAAACGCAGGCgctgg                             <  3:288194/149‑1 (MQ=255)
                                                                                                                       cgcCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCGAAGTGGCTGAAGTCGAAACGCAGGCgctgg                             <  3:67954/149‑1 (MQ=255)
                                                                                                                       cgcCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTAGATCTGCTCCGGCTTCACCGCTTCGAAGTGGCTGAAGTCGAAACGCAGGCgctgg                             <  1:283486/149‑1 (MQ=255)
                                                                                                                        gcCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGATTCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCGAAGTGGCTGAAGTCGAAACGCAGGCgctggc                            <  2:61765/149‑1 (MQ=255)
                                                                                                                           tttGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCggg                                                                                                                                <  3:44217/46‑1 (MQ=255)
                                                                                                                           tttGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCggg                                                                                                                                >  4:44217/1‑46 (MQ=255)
                                                                                                                                    ggCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCGAAGTGGCTGAAGTCGAAACGCAGGCGCTGGCTGTCGACCAGc                 <  7:178602/148‑1 (MQ=255)
                                                                                                                                     gCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCGAAGTGGCTGAAGTCGAAACGCAGGCGCTGGCTGTCGACCAGCGa               <  3:292972/149‑1 (MQ=255)
                                                                                                                                      cAGTTTCGATGTCGGTCAGGTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCGAAGTGGCTGAAGTCGAAACGCAGGCGCTGGCTGTCGACCAGCGAg              <  4:146204/149‑1 (MQ=255)
                                                                                                                                            cGATGTCGATCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCGAAGTGGCTGAAGTCGAAACGCAGGCGCTGGCTGTCGACCAGCGAGCCCtt         <  5:214747/148‑1 (MQ=255)
                                                                                                                                            cGATGTCGATCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCGAAGTGGCTGAAGTCGAAACGCAGGCGCTGGCTGTCGACCAGCGAGCCCtt         <  5:219414/148‑1 (MQ=255)
                                                                                                                                                 tCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCGAAGTGGCTGAAGTCGAAACGCAGGCGCTGGCTGTCGACCAGCGAGCCCTTctgctg   >  6:229980/1‑149 (MQ=255)
                                                                                                                                                  cGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGtt                                                                                                         <  2:204305/46‑1 (MQ=255)
                                                                                                                                                  cGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGtt                                                                                                         >  1:204305/1‑46 (MQ=255)
                                                                                                                                                   ggTCAGTTCTGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCGAAGTGGCTGAAGTCGAAACGCAGGCGCTGGCTGTCGACCAGCGAGCCCTTCTGCTGt  <  8:70209/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAGGCTGATGTCACCGGTGCGCTTGGCGTGGATACCGCCGCACAATTCCACCGAAAAATCACCGCCCATGCTCAGCACACGCACGGTGTCGCCGTATTTCTCGCCGAACAGCGCCATTGCGCCTTTGGCCTTGGCAGTTTCGATGTCGGTCAGTTCGGTTTCAACCGGGGTATTCTTGCGAATCTCGCGGTTGACGATGTCTTCCAGCTGCTTGATCTGCTCCGGCTTCACCGCTTCGAAGTGGCTGAAGTCGAAACGCAGGCGCTGGCTGTCGACCAGCGAGCCCTTCTGCTGT  >  NC_002947/5082232‑5082526

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: