Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 476,701 G→T 8.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,7010GT 8.2% 259.8 / 9.4 85Q330H (CAG→CAT) ‡tsaDtRNA N6‑adenosine(37)‑threonylcarbamoyltransferase complex transferase subunit TsaD
Reads supporting (aligned to +/- strand):  ref base G (39/39);  new base T (2/5);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

GCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGATCGGCAGCAG  >  NC_002947/476552‑476847
                                                                                                                                                     |                                                                                                                                                  
gcgTGCCGGGAGGACATGTTCGGTAGTATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAg                                                                                                                                                    <  2:424414/150‑1 (MQ=255)
gcgTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAg                                                                                                                                                    <  2:1941172/150‑1 (MQ=255)
gcgTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAg                                                                                                                                                    <  1:1645198/150‑1 (MQ=255)
gcgTCCCTGGAGGACATGCTCAGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAg                                                                                                                                                    <  1:82046/150‑1 (MQ=255)
 cgTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAgg                                                                                                                                                   <  1:1728674/150‑1 (MQ=255)
   tCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCCGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCa                                                                                                                                                 >  2:1564030/1‑150 (MQ=255)
      cTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCCTGCACGc                                                                                                                                              >  1:958750/1‑150 (MQ=255)
      cTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGc                                                                                                                                              <  2:1312338/150‑1 (MQ=255)
      cTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGc                                                                                                                                              >  2:312675/1‑150 (MQ=255)
       tGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCt                                                                                                                                             >  1:1863255/1‑150 (MQ=255)
         gaggaCATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGACCAGCGTGCCTGCACGCTga                                                                                                                                           <  1:1006345/150‑2 (MQ=255)
          aggaCATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCCTGCACGCTgat                                                                                                                                          <  2:958750/150‑3 (MQ=255)
          aggaCATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCCTGCACGCTgat                                                                                                                                          <  1:970779/150‑3 (MQ=255)
          aggaCATGCTCGGTAGCATCAAGGGAAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAACCCGACCTGGCGATCAGCGTGCAGGCACGCTGGc                                                                                                                                          <  1:2177393/150‑1 (MQ=255)
           ggaCATGCTCGGTAGCATCAAGCGCAATGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGcc                                                                                                                                         >  2:2020337/1‑150 (MQ=255)
            gaCATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCCTGCACGCTgatcg                                                                                                                                        <  1:1042092/150‑5 (MQ=255)
               aTGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATg                                                                                                                                     <  2:1863255/150‑1 (MQ=255)
                  cTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGaa                                                                                                                                  <  2:258666/150‑1 (MQ=255)
                   tCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAAc                                                                                                                                 >  2:1369352/1‑150 (MQ=255)
                   tCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTAACGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCGGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCGGGCACGCTGGCCGATGGaag                                                                                                                                 >  2:861405/1‑149 (MQ=255)
                      gTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGt                                                                                                                              <  1:1393995/150‑1 (MQ=255)
                      gTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGt                                                                                                                              <  1:312675/150‑1 (MQ=255)
                       tAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGtt                                                                                                                             <  1:1369352/150‑1 (MQ=255)
                       tAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCCGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGtt                                                                                                                             <  1:1564030/150‑1 (MQ=255)
                         gCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTgc                                                                                                                           <  2:1279232/150‑1 (MQ=255)
                         gCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTgc                                                                                                                           <  2:479866/150‑1 (MQ=255)
                          cATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTgcc                                                                                                                          >  1:904080/1‑150 (MQ=255)
                           aTCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTAGCGATCGGCGTGCAGGCACGCTGGCCGATGGAACAGTTgccg                                                                                                                         >  1:538945/1‑150 (MQ=255)
                            tCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGGAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTgccgc                                                                                                                        >  2:302415/1‑150 (MQ=255)
                            tCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCCTGCACGCTgatc                                                                                                                                         <  1:1742248/133‑4 (MQ=255)
                            tCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCCTGCACGCTgatc                                                                                                                                         >  2:1742248/1‑130 (MQ=255)
                               aGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCCGGCTTGCCGATGGAACAGTTGCcgccgc                                                                                                                     >  1:552/1‑150 (MQ=255)
                                cagCAACGTGTACTACGCGCGTCCGCAACTATGTACCGAGAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGCTCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCcaccgca                                                                                                                    >  2:538482/3‑149 (MQ=255)
                                        tgtACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTa                                                                                                            <  1:1726282/150‑1 (MQ=255)
                                                  gcgTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGaa                                                                                                  <  2:951203/150‑1 (MQ=255)
                                                   cgTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAg                                                                                                 <  1:1138069/150‑1 (MQ=255)
                                                     tCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGGAAGACCTAGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGt                                                                                               >  1:299120/1‑150 (MQ=255)
                                                        gCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCa                                                                                            >  2:678188/1‑150 (MQ=255)
                                                            ttCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAg                                                                                        >  2:1997492/1‑150 (MQ=255)
                                                                gTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGcc                                                                                    <  1:1117510/150‑1 (MQ=255)
                                                                 tACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCg                                                                                   >  1:1018694/1‑150 (MQ=255)
                                                                  aCCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGt                                                                                  >  1:1207871/1‑150 (MQ=255)
                                                                   ccGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGtt                                                                                 >  2:1968200/1‑150 (MQ=255)
                                                                     gACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTcg                                                                               >  2:526506/1‑150 (MQ=255)
                                                                      aCAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTcgc                                                                              >  2:1234980/1‑150 (MQ=255)
                                                                         aTGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCTGTTcgccca                                                                           <  2:31468/150‑4 (MQ=255)
                                                                          tGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGATGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGcgcg                                                                                                       >  2:237836/1‑121 (MQ=255)
                                                                          tGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGATGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGcgcg                                                                                                       <  1:237836/121‑1 (MQ=255)
                                                                              gCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCa                                                                      >  2:1632573/1‑150 (MQ=255)
                                                                                gatgatCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCaaa                                                                    <  2:2044804/150‑1 (MQ=255)
                                                                                     tCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGt                                                               >  1:1605220/1‑150 (MQ=255)
                                                                                      cGCCTACGCCGGTTGCCAGCGCCTGCTGGTCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTc                                                              <  1:208020/150‑1 (MQ=255)
                                                                                      cGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTc                                                              <  1:707584/150‑1 (MQ=255)
                                                                                       gCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGGAAGACCTAGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTcg                                                             <  2:299120/150‑1 (MQ=255)
                                                                                        ccTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTcgc                                                            <  1:1767273/150‑1 (MQ=255)
                                                                                        ccTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTcgc                                                            >  2:1520673/1‑150 (MQ=255)
                                                                                            cGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAg                                                        <  2:728389/150‑1 (MQ=255)
                                                                                                ggTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGc                                                                                                                                                  <  2:2191333/56‑1 (MQ=255)
                                                                                                ggTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGc                                                                                                                                                  >  1:2191333/1‑56 (MQ=255)
                                                                                                     ccAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCg                                               >  1:1549376/1‑150 (MQ=255)
                                                                                                      cAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCgg                                              <  1:1674742/150‑1 (MQ=255)
                                                                                                      cAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTAAGAAGTGTCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCgg                                              >  1:161748/1‑150 (MQ=255)
                                                                                                       aGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGt                                             <  1:1520673/150‑1 (MQ=255)
                                                                                                       aGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGt                                             >  2:979826/1‑150 (MQ=255)
                                                                                                        gcgcCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTg                                            >  2:1841677/1‑150 (MQ=255)
                                                                                                          gcCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGgc                                          <  1:1841677/150‑1 (MQ=255)
                                                                                                            ctgctgGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTaa                                                                                                                <  1:2032732/78‑1 (MQ=255)
                                                                                                            ctgctgGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTaa                                                                                                                >  2:2032732/1‑78 (MQ=255)
                                                                                                            ctgctgGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGgcgc                                        >  2:1506295/1‑150 (MQ=255)
                                                                                                            ctgctgGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGgcgc                                        <  1:1506295/150‑1 (MQ=255)
                                                                                                              gctgGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGATGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCca                                      >  1:1688112/1‑150 (MQ=255)
                                                                                                                tgGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCcaca                                    >  1:2009587/1‑150 (MQ=255)
                                                                                                                 ggCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGGCCAGCAAACAAGTCGCGCAGGTTGTTGCGGTGGAGCCacac                                   >  1:1832608/1‑150 (MQ=255)
                                                                                                                     gggCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGatc                               <  1:1632573/150‑1 (MQ=255)
                                                                                                                       gcagcagcaAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGatcat                             >  2:107237/1‑150 (MQ=255)
                                                                                                                        cagcagcaAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGatcatc                            <  2:1605220/150‑1 (MQ=255)
                                                                                                                           cagcaAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACt                         <  2:1018694/150‑1 (MQ=255)
                                                                                                                             gcaAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGt                       <  1:526506/150‑1 (MQ=255)
                                                                                                                                     tGGCGATCAGCGTGCAGGCCCGCTGGCCGATGGAACATTTGCCGCCGCTGTAAGGTATCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGc               <  1:861405/150‑1 (MQ=255)
                                                                                                                                     tGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGc               <  2:2009587/150‑1 (MQ=255)
                                                                                                                                     tGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCAAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGACGCCACACTGTCATCACTGTCAGTACGc               >  2:127165/1‑150 (MQ=255)
                                                                                                                                         gATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGAt           <  1:678188/150‑1 (MQ=255)
                                                                                                                                         gATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGAt           >  2:984250/1‑150 (MQ=255)
                                                                                                                                          aTCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTCTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGTTCATCACGGTCAGTACGCTGATc          <  2:1735706/150‑1 (MQ=255)
                                                                                                                                           tCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACgatgat           >  2:135434/1‑148 (MQ=255)
                                                                                                                                           tCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGAt           <  1:135434/148‑1 (MQ=255)
                                                                                                                                             aGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGATCGgc       >  2:1134157/1‑150 (MQ=255)
                                                                                                                                                gTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGATCGgcagc    >  1:280241/1‑150 (MQ=255)
                                                                                                                                                 tGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAATGTTCAGAAGTGCCGTTCGAGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGATCGgcagca   <  1:1968200/150‑1 (MQ=255)
                                                                                                                                                  gCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGCTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGATCGGcagcag  <  1:1234980/150‑1 (MQ=255)
                                                                                                                                                     ggCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGatcatc                            >  1:2102452/1‑121 (MQ=255)
                                                                                                                                                     ggCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGatcatc                            <  2:2102452/121‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCCACACGATCATCACTGTCAGTACGCTGATCGGCAGCAG  >  NC_002947/476552‑476847

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