Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,384,157 T→G 14.3% I292L (ATT→CTT)  proS ← proline‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,384,1570TG14.3% 159.7 / 20.6 63I292L (ATT→CTT) proSproline‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (26/28);  new base G (4/5);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

GGTTGAGCGGGCCGAGCGAGCCGGCGCCGGCGCCAATGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCAATGGTCTTGGCATTCGGGGTGTCGACCAGGCGCAGCTCCTCGGTAGGGGCAGGGCGCAC  >  NC_002947/1384010‑1384301
                                                                                                                                                   |                                                                                                                                                  
ggTTGAGCGGGCCGAGCGAGCCGGCGCCGGCGCCAATGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGAGTTCGTTCAGTTCGTGGTCGCCACGGCCGATCAGGGCAAt                                                                                                                                                    >  2:315203/1‑148 (MQ=255)
     aGCGGGCCGAGCGAGCCGGCGCCGGCGCCAATGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCtt                                                                                                                                              <  1:374631/149‑1 (MQ=255)
      gCGGGCCGAGCGAGCCGGCGCCGGCGCCAATGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTg                                                                                                                                             <  6:228433/149‑1 (MQ=255)
       cGGGCCGAGCGAGCCGGCGCCGGCGCCAATGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGc                                                                                                                                            <  4:124431/149‑1 (MQ=255)
        gggCCGAGCGAGCCGGCGCCGGCGCCAATGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGcc                                                                                                                                           <  2:146972/149‑1 (MQ=255)
        gggCCGAGCGAGCCGGCGCCGGCGCCAATGGCGTCGCGCAGTTCGGCTTCGGTGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGccc                                                                                                                                          <  3:352558/149‑1 (MQ=255)
          gCCGAGCGAGCCGGCGCCGGCGCCAATGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCct                                                                                                                                         >  4:156085/1‑149 (MQ=255)
          gCCGAGCGAGCCGGCGCCGGCGCCAATGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCct                                                                                                                                         >  8:336163/1‑149 (MQ=255)
           ccGAGCGAGCCGGCGCCGGCGCCAATGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTCGTGGTCGCCACGGACGATCAGGGCAAGCTGATTGCCctga                                                                                                                                       <  8:207387/149‑3 (MQ=255)
             gagcgagcCGGCGCCGGCGCCAATGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAAGCTGATTGCCctgat                                                                                                                                      <  1:80411/149‑4 (MQ=255)
              agcgagcCGGCGCCGGCGCCAATGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAAGCTGATTGCCctgatc                                                                                                                                     <  6:311964/149‑5 (MQ=255)
               gcgagcCGGCGCCGACGCCAATGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCg                                                                                                                                    <  5:1970/149‑1 (MQ=255)
                cgagcCGGCGCCGGCGCCAATGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCg                                                                                                                                    <  3:253824/148‑1 (MQ=255)
                 gagcCGGCGCCGGCGCCAATGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTg                                                                                                                                             <  3:294238/138‑1 (MQ=255)
                    ccggcgccggcgCCAATGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCAcc                                                                                                                               >  4:128522/1‑149 (MQ=255)
                      ggcgccggcgccAATGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACcgc                                                                                                                             <  8:19410/149‑1 (MQ=255)
                                 cAATGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGTACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCa                                                                                                                    <  8:158243/147‑1 (MQ=255)
                                   aTGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGt                                                                                                                <  6:111488/149‑1 (MQ=255)
                                        gTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGCCGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTgc                                                                                                           >  5:361289/1‑148 (MQ=255)
                                        gTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGa                                                                                                           >  4:153455/1‑149 (MQ=255)
                                        gTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGa                                                                                                           <  1:112783/149‑1 (MQ=255)
                                           gcgcAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACgg                                                                                                        >  6:131981/1‑149 (MQ=255)
                                              cAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCt                                                                                                    >  8:375665/1‑149 (MQ=255)
                                                 ttCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGTTCtttt                                                                                                 >  3:327515/1‑149 (MQ=255)
                                                 ttCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCtttt                                                                                                 <  4:398664/149‑1 (MQ=255)
                                                  tCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACGCGTTCGAGGTTGGTGGCCTTGATTTAGTTCAGTTCGTGGTCGCGAGGGACGATCAGGGCAATCATCTTGCTGTCTTCGGCACCGCTTACGATCAGGGTCTTGACGGTCtttt                                                                                                 <  4:243522/148‑1 (MQ=255)
                                                             ggCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGc                                                                                     >  6:141036/1‑149 (MQ=255)
                                                                cATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGt                                                                                  >  4:334281/1‑149 (MQ=255)
                                                                cATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCAACAGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGt                                                                                  >  7:76337/1‑149 (MQ=255)
                                                                 aTGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTg                                                                                 <  1:21356/149‑1 (MQ=255)
                                                                 aTGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTg                                                                                 <  7:375665/149‑1 (MQ=255)
                                                                      cAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGtttt                                                                            >  8:16879/1‑149 (MQ=255)
                                                                      cAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGCCCGTGGtttt                                                                            >  7:172235/1‑149 (MQ=255)
                                                                               ggCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACt                                                                   >  2:41038/1‑149 (MQ=255)
                                                                                gCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGACGTGGTTTTCCACCAACTg                                                                  >  3:392787/1‑149 (MQ=255)
                                                                                 cAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTgg                                                                 <  5:131981/149‑1 (MQ=255)
                                                                                     cTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCaa                                                             >  6:92659/1‑149 (MQ=255)
                                                                                          ccAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTc                                                                                                     <  4:151999/104‑1 (MQ=255)
                                                                                          ccAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTc                                                                                                     >  3:151999/1‑104 (MQ=255)
                                                                                             gCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCAATGGTCTTg                                                     <  1:41038/149‑1 (MQ=255)
                                                                                                      ccTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTg                                                                                                            <  6:382394/86‑1 (MQ=255)
                                                                                                      ccTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTg                                                                                                            >  5:376303/1‑86 (MQ=255)
                                                                                                      ccTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTg                                                                                                            >  5:382394/1‑86 (MQ=255)
                                                                                                      ccTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTg                                                                                                            <  6:376303/86‑1 (MQ=255)
                                                                                                      ccTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCAATGGTCTTGGCATTCggg                                           <  3:153455/149‑1 (MQ=255)
                                                                                                      ccTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCAATGGTCTTGGCATTCggg                                           <  3:334281/149‑1 (MQ=255)
                                                                                                       cTTGATTTCGTTCAGTTCGGGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGGGTACGATCAGGGGCTTGACGGTATTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCAATGGTCTTGGCATTCggg                                           <  1:315203/149‑1 (MQ=255)
                                                                                                        ttGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCAATGGTCTTGGCATTCGGGgt                                         >  5:277715/1‑149 (MQ=255)
                                                                                                          gATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGt                                                                                  <  8:219035/107‑1 (MQ=255)
                                                                                                          gATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGt                                                                                  >  7:219035/1‑107 (MQ=255)
                                                                                                             ttCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCAATGGTCTTGGCATTCGGGGTGTCGa                                    >  4:14953/1‑149 (MQ=255)
                                                                                                                gTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCAATGGTCTTGGCATTCGGGGTGTCGACCa                                 <  3:128522/149‑1 (MQ=255)
                                                                                                                     gTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCAATGGTCTTGGCATTCGGGGTGTCGACCAGgcg                             >  1:318642/1‑148 (MQ=255)
                                                                                                                          tGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACAAACTGGGCAATGGTCTTGGCATTCGGGGTGTCGACCAGGCGCAGctc                       >  1:327555/1‑149 (MQ=255)
                                                                                                                               gCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCAATGGTCTTGGCATTCGGGGTGTCGACCAGGCGCAGCTCCTCgg                  >  5:9498/1‑149 (MQ=255)
                                                                                                                                    ggACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCAATGGTCTTGGCATTCGGGGTGTCGACCAGGCGCAGCTCCTCGGTAggg             <  2:327555/149‑1 (MQ=255)
                                                                                                                                        gaagaGGGCAAGCTGATTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCGACCAACAGGGCAATGGTCTTGGCATTCGGGGTGTCGACCAGGCGCAGCTCCTCGGTAGGGGCAg         >  5:1326/5‑149 (MQ=255)
                                                                                                                                        gaagaGGGCAAGCTGATTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCAATGGTCtt                                                      >  8:199407/5‑105 (MQ=255)
                                                                                                                                        gaagaGGGCAAGCTGATTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCAATGGTCTTGGCCATTCGGGGTGTCGACCAGGCGCAGc                         <  1:209974/130‑1 (MQ=255)
                                                                                                                                        gaagaGGGCAAGCTGATTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCAATGGTCTTGGCATTCGGGGTGTCGACCAGGCGCAGc                         >  2:209974/5‑133 (MQ=255)
                                                                                                                                          agaGGGCAAGCTGATTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCAATGGTCTTGGCATTCGGGGTGTc                                      >  4:383898/3‑118 (MQ=255)
                                                                                                                                          agaGGGCAAGCTGATTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCAATGGTCTTGGCATTCGGGGTGTc                                      <  3:383898/116‑1 (MQ=255)
                                                                                                                                              ggCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCAATGGTCTTGGCATTCGGGGTGTCGACCAGGCGCAGCTCCTCGGTAGGGGCAGGGCGCa   >  6:148466/1‑149 (MQ=255)
                                                                                                                                               gCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTCGTTTTCCACCAACTGGGCAATGGTCTTGGCATTCGGGGTGTCGACCAGGCGCAGCTCCTCGGTAGGGGCAGGGCGCAc  <  5:92659/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGTTGAGCGGGCCGAGCGAGCCGGCGCCGGCGCCAATGGCGTCGCGCAGTTCGGCTTCGGTGGCCATGACCAGCGGGTCGGCAACCTGTTCCAGCTTGGTGGCCTTGATTTCGTTCAGTTCGTGGTCGCCACGGACGATCAGGGCAATCAGCTTGCCCTCTTCGGCACCGCGTACGATCAGGGTCTTGACGGTCTTTTCGATCGGCAGGCCGTGGTTTTCCACCAACTGGGCAATGGTCTTGGCATTCGGGGTGTCGACCAGGCGCAGCTCCTCGGTAGGGGCAGGGCGCAC  >  NC_002947/1384010‑1384301

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