Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 476,701 G→T 6.8% 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,7010GT 6.8% 183.2 / 2.5 59Q330H (CAG→CAT) ‡tsaDtRNA N6‑adenosine(37)‑threonylcarbamoyltransferase complex transferase subunit TsaD
Reads supporting (aligned to +/- strand):  ref base G (31/24);  new base T (2/2);  total (33/26)
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.19e-01

TGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGAT  >  NC_002947/476559‑476838
                                                                                                                                              |                                                                                                                                          
tGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAg                                                                                                                                            <  1:463945/143‑1 (MQ=255)
tGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAg                                                                                                                                            >  2:463945/1‑143 (MQ=255)
   aggaCATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGc                                                                                                                                  <  1:478417/150‑1 (MQ=255)
        aTGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATg                                                                                                                             >  2:98719/1‑150 (MQ=255)
        aTGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATg                                                                                                                             >  1:1157918/1‑150 (MQ=255)
        aTGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATg                                                                                                                             <  1:616798/150‑1 (MQ=255)
        aTGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATg                                                                                                                             >  2:535876/1‑150 (MQ=255)
           cTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGaa                                                                                                                          <  1:144221/150‑1 (MQ=255)
            tCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAAc                                                                                                                         >  1:354333/1‑150 (MQ=255)
               gTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGt                                                                                                                      >  1:731317/1‑150 (MQ=255)
                         gggCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCt                                                                                                            >  2:1057087/1‑150 (MQ=255)
                            cAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGCGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTa                                                                                                         >  2:1377114/1‑150 (MQ=255)
                                  gtACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAg                                                                                                   >  2:293340/1‑150 (MQ=255)
                                     ctacGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAgcgc                                                                                                >  2:42622/1‑150 (MQ=255)
                                     ctacGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATACGTCGGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAgcgc                                                                                                >  1:448496/1‑150 (MQ=255)
                                      tacGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTACCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGcgcg                                                                                               >  2:1262819/1‑150 (MQ=255)
                                         gcgcgTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCg                                                                                            >  1:1004361/1‑150 (MQ=255)
                                            cgTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAg                                                                                         <  2:677758/150‑1 (MQ=255)
                                            cgTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAg                                                                                         <  2:1028852/150‑1 (MQ=255)
                                              tCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGt                                                                                       >  1:143205/1‑150 (MQ=255)
                                               ccGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGtt                                                                                      <  1:42622/150‑1 (MQ=255)
                                                cGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGGGCTTGCTCGCCGGGTGGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCGCGCCGATGGAACGGTGGCCACCGCCGTTAGGTAGCGCGTCGTAgggga                                                                                     >  2:555329/1‑147 (MQ=255)
                                                cGCAATTCTGTACCGACAATGGTGCGATCATCGCCTACGCCGGTGGGCAGCGCCTTCTCGTGCCGCAGCAGCCCGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACCGTTGCCGACGCTGGAAGGTAGCGCGTCGAAGGTTc                                                                                     >  1:535304/1‑150 (MQ=255)
                                                 gCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATgg                                                                                                                            <  1:963366/110‑1 (MQ=255)
                                                 gCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATgg                                                                                                                            >  2:963366/1‑110 (MQ=255)
                                                      tCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGt                                                                               <  1:273832/150‑1 (MQ=255)
                                                         gTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATgg                                                                                                                            <  2:796642/102‑1 (MQ=255)
                                                         gTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATgg                                                                                                                            >  1:796642/1‑102 (MQ=255)
                                                                 aaTGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCgcgc                                                                   >  1:52105/1‑150 (MQ=255)
                                                                         gatgatCGCCTACGCCGGTTGGCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCaaa                                                           <  2:731317/150‑1 (MQ=255)
                                                                         gatgatCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCaaa                                                           <  1:1057087/150‑1 (MQ=255)
                                                                            gatCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGtt                                                                                      >  1:1377459/1‑121 (MQ=255)
                                                                            gatCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGtt                                                                                      <  2:1377459/121‑1 (MQ=255)
                                                                            gatCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACaa                                                        <  1:1377114/150‑1 (MQ=255)
                                                                            gatCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACaa                                                        <  1:790914/150‑1 (MQ=255)
                                                                            gatCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACaa                                                        <  1:535876/150‑1 (MQ=255)
                                                                              tCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGGTCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTACGGTAGGGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGt                                                      >  2:855652/1‑150 (MQ=255)
                                                                                    aCGCCGGTTGCCAGCGCCTACTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGTCGTTCGCGCCCTGCAAACAAGTCGCGCa                                                <  2:939993/150‑1 (MQ=255)
                                                                                      gCCGGTTGCCAGCGCCTGCTGGTCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGAAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAgg                                              >  1:807546/1‑150 (MQ=255)
                                                                                      gCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAgg                                              <  2:293396/150‑1 (MQ=255)
                                                                                      gCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAgg                                              <  1:990979/150‑1 (MQ=255)
                                                                                       ccGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCAGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGgt                                             <  1:293340/150‑1 (MQ=255)
                                                                                              gcAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGCACAGTTGCCGCCGCAGTAAGGTAGCGCGTCAAAGGGTGCGAAGTTCCGTAGATCTCCGTGAAGCAAGTCGCGCTGGTTGTTGCg                                      >  2:514659/2‑150 (MQ=255)
                                                                                               cAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCgg                                     >  1:1124528/1‑150 (MQ=255)
                                                                                               cAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCgg                                     >  1:227927/1‑150 (MQ=255)
                                                                                                     ctgctgGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGgcgc                               >  2:1104470/1‑150 (MQ=255)
                                                                                                                   gcagcaAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAAGAAGCGCGTCGAAGAACCGTAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTCCGGTGGCGACACACGATCATCAc                 >  1:1156551/1‑150 (MQ=255)
                                                                                                                    cagcaAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACt                <  2:143205/150‑1 (MQ=255)
                                                                                                                     agcaAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGt              >  1:1088778/1‑150 (MQ=255)
                                                                                                                      gcaAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGt              >  1:1364939/1‑150 (MQ=255)
                                                                                                                       caAGACCTGGCGATCAGCGTGCAGGCACGCTGGTCGATGGAATAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTc             <  2:662997/150‑1 (MQ=255)
                                                                                                                             cTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACg       >  1:1063701/1‑150 (MQ=255)
                                                                                                                              tGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGc      >  1:961965/1‑150 (MQ=255)
                                                                                                                                 cGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGAt  <  2:1088778/150‑1 (MQ=255)
                                                                                                                                  gATCAGCGTGCAGGCACGCTGGCCGCTGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGAt  <  1:1161817/150‑1 (MQ=255)
                                                                                                                                   gccAGCGTGCCTGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGaa                                                                                          <  2:621991/60‑1 (MQ=38)
                                                                                                                                   gccAGCGTGCCTGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGaa                                                                                          >  1:621991/3‑62 (MQ=38)
                                                                                                                                   gccAGCGTGCCTGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGaa                                                                                          >  1:975664/3‑62 (MQ=38)
                                                                                                                                   gccAGCGTGCCTGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGaa                                                                                          <  2:975664/60‑1 (MQ=38)
                                                                                                                                             aGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCa            <  1:1334288/129‑1 (MQ=255)
                                                                                                                                             aGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCa            >  2:1334288/1‑129 (MQ=255)
                                                                                                                                              |                                                                                                                                          
TGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGAT  >  NC_002947/476559‑476838

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