Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 476,698 G→T 7.2% Q330H (CAG→CAT) ‡ tsaD → tRNA N6‑adenosine(37)‑threonylcarbamoyltransferase complex transferase subunit TsaD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947476,6980GT 7.2% 213.0 / 4.0 67Q330H (CAG→CAT) ‡tsaDtRNA N6‑adenosine(37)‑threonylcarbamoyltransferase complex transferase subunit TsaD
Reads supporting (aligned to +/- strand):  ref base G (36/26);  new base T (2/3);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGATCGGCAGCAGG  >  NC_002947/476550‑476845
                                                                                                                                                    |                                                                                                                                                   
cgTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAg                                                                                                                                                     <  2:64376/149‑1 (MQ=255)
    ccTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCAc                                                                                                                                                 >  1:16753/1‑149 (MQ=255)
    ccTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACCCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCAc                                                                                                                                                 >  4:112712/1‑149 (MQ=255)
      tGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGc                                                                                                                                               >  3:6551/1‑149 (MQ=255)
       ggaggaCATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACg                                                                                                                                                >  1:312195/1‑147 (MQ=255)
       ggaggaCATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACg                                                                                                                                                <  2:312195/147‑1 (MQ=255)
         aggaCATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTgg                                                                                                                                            <  2:233294/149‑1 (MQ=255)
          ggaCATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTCCCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCCTGCACGCTgat                                                                                                                                           <  3:67997/149‑3 (MQ=255)
          ggaCATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCCTGCACGCTgat                                                                                                                                           <  2:172965/149‑3 (MQ=255)
          ggaCATGCGCGGTAGCATCAAGGGCACCGTGTCCTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCCTGCACGCTgat                                                                                                                                           <  3:63640/149‑3 (MQ=255)
             cATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGa                                                                                                                                        >  5:187112/1‑149 (MQ=255)
              aTGCTCGGTAGCATCAAGGGCACCGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGAt                                                                                                                                       <  4:20391/149‑1 (MQ=255)
              aTGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGAt                                                                                                                                       <  3:295797/149‑1 (MQ=255)
                      tAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGt                                                                                                                               <  2:240571/149‑1 (MQ=255)
                          aTCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTgcc                                                                                                                           >  2:112901/1‑149 (MQ=255)
                           tCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGc                                                                                                                                                   <  4:62919/124‑1 (MQ=255)
                           tCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGc                                                                                                                                                   >  3:62919/1‑124 (MQ=255)
                             aaGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTgccgcc                                                                                                                        >  6:46524/1‑149 (MQ=255)
                                ggCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCt                                                                                                                     >  3:235306/1‑149 (MQ=255)
                                       tgtACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCCTACACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAgg                                                                                                               >  1:109558/1‑148 (MQ=255)
                                       tgtACTACGCGCGGCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGt                                                                                                              <  8:223598/149‑1 (MQ=255)
                                            tacGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAgcgc                                                                                                         >  1:70807/1‑149 (MQ=255)
                                              cgcgcgTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGt                                                                                                       >  5:274226/1‑149 (MQ=255)
                                                   gTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAg                                                                                                  <  5:119752/149‑1 (MQ=255)
                                                      cGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGTTGGCCGATGGAACAGATGCCGCCGCTGTAAGGTAGCGCGTCGAAGGtt                                                                                               >  5:68631/1‑149 (MQ=255)
                                                      cGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGtt                                                                                               >  3:120187/1‑149 (MQ=255)
                                                      cGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGt                                                                                                <  6:86049/148‑1 (MQ=255)
                                                       gCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTc                                                                                              <  8:18218/149‑1 (MQ=255)
                                                         aaTTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAgg                                                                                                               <  7:187209/130‑1 (MQ=255)
                                                         aaTTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAgg                                                                                                               >  8:187209/1‑130 (MQ=255)
                                                                  ccGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGt                                                                                   <  6:187112/149‑1 (MQ=255)
                                                                  ccGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGACGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGt                                                                                   >  2:159078/1‑149 (MQ=255)
                                                                    gACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTc                                                                                 <  7:225812/149‑1 (MQ=255)
                                                                       aaTGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTcgcg                                                                              >  6:221604/1‑149 (MQ=255)
                                                                              cGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGc                                                                        <  6:274226/148‑1 (MQ=255)
                                                                                    tCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAg                                                                 >  2:125967/1‑149 (MQ=255)
                                                                                    tCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGAAAGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAg                                                                 >  1:102283/1‑149 (MQ=255)
                                                                                             ccGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCACGCCCTGCAAACAAGTCGCGCAgg                                                        >  6:255572/1‑149 (MQ=255)
                                                                                                     gtcttgatgctgcCCGGCCAGCAGCAAGACCTGGCGATAAGCGTGCAGGCACGCGGGCGGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTGGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAACTCGCGCAGGTTGTTGCg                                                >  2:184680/8‑149 (MQ=255)
                                                                                                     cAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCCTGCAGGCACGCTGGCCGATTAAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCg                                                >  2:321193/1‑149 (MQ=255)
                                                                                                      aGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAAc                                                                    >  6:287417/1‑128 (MQ=255)
                                                                                                      aGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAAc                                                                    <  5:287417/128‑1 (MQ=255)
                                                                                                      aGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCgg                                               <  2:109558/149‑1 (MQ=255)
                                                                                                          gctgctgCCCGGCCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGgc                                           >  8:117482/2‑149 (MQ=255)
                                                                                                             gctgGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGgcgc                                         >  4:70527/1‑148 (MQ=255)
                                                                                                              ctgGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGTCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGgcgc                                         <  7:2829/147‑1 (MQ=255)
                                                                                                                  ccGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACg                                   <  3:112712/149‑1 (MQ=255)
                                                                                                                    gggCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGat                                 <  3:95773/149‑1 (MQ=255)
                                                                                                                     ggCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGatc                                <  1:159078/149‑1 (MQ=255)
                                                                                                                       cagcagcaAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGatcat                              <  5:46524/149‑1 (MQ=255)
                                                                                                                        agcagcaAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGatcat                              >  2:1357/1‑148 (MQ=255)
                                                                                                                         gcagcaAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGAtcatca                            >  4:228154/1‑149 (MQ=255)
                                                                                                                          cagcaAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCAc                           >  4:59514/1‑149 (MQ=255)
                                                                                                                           agcaAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGatc                                >  7:80558/1‑143 (MQ=255)
                                                                                                                           agcaAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGatc                                <  8:80558/143‑1 (MQ=255)
                                                                                                                           agcaAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGGCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACt                          <  8:243546/149‑1 (MQ=255)
                                                                                                                               aGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCa                      <  3:228154/149‑1 (MQ=255)
                                                                                                                                   cTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGgc                                           <  6:288740/124‑1 (MQ=255)
                                                                                                                                   cTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGgc                                           >  5:288740/1‑124 (MQ=255)
                                                                                                                                      gCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGc                <  2:102283/148‑1 (MQ=255)
                                                                                                                                      gCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGAAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCt               >  4:187775/1‑149 (MQ=255)
                                                                                                                                         gccAGCGTGCCTGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGAt            >  4:281244/3‑149 (MQ=255)
                                                                                                                                           cAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGGCATTACGCTGATCg          >  8:284551/1‑149 (MQ=255)
                                                                                                                                               gTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGATCGgcag      >  7:274760/1‑149 (MQ=255)
                                                                                                                                               gTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGATCGgcag      >  2:199286/1‑149 (MQ=255)
                                                                                                                                                 gCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGATCGgcag      >  8:239323/1‑147 (MQ=255)
                                                                                                                                                  cAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGATCGGcagcag   >  8:184691/1‑149 (MQ=255)
                                                                                                                                                   aGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGATCGGCAGCAgg  >  6:205701/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGATCGGCAGCAGG  >  NC_002947/476550‑476845

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: