Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,737,658 G→C 16.8% intergenic (‑157/+45) hutU ← / ← PP_5034 urocanate hydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,737,6580GC16.8% 238.8 / 39.8 95intergenic (‑157/+45)hutU/PP_5034urocanate hydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (39/40);  new base C (9/7);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 102 total aligned reads displayed.

GTCCTCTGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAACCCCAA  >  NC_002947/5737509‑5737801
                                                                                                                                                     |                                                                                                                                                
gTCCTCTGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCgcgg                                                                                                                                                  <  1:1997174/150‑1 (MQ=255)
 tCCTCTGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCgcggg                                                                                                                                                 <  1:51483/150‑1 (MQ=255)
 tCCTCTGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCgcggg                                                                                                                                                 >  1:2155561/1‑150 (MQ=255)
 tCCTCTGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCgcggg                                                                                                                                                 <  1:945697/150‑1 (MQ=255)
   ctctGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCgcgggcg                                                                                                                                               <  1:1419415/150‑1 (MQ=255)
    tctGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCgcgggcgg                                                                                                                                              >  1:2750096/1‑150 (MQ=255)
     ctGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGggcggc                                                                                                                                             <  2:2750096/150‑1 (MQ=255)
       gCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTCGGGGTATAGATCGAGCGCCGCCCGCGCGgcgc                                                                                                                                           <  1:2106070/150‑1 (MQ=255)
         ggTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACGAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                         >  1:1510440/1‑150 (MQ=255)
         ggTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                         >  2:260354/1‑150 (MQ=255)
         ggTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                         <  1:821783/150‑1 (MQ=255)
            cgtcCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                      <  1:260354/150‑1 (MQ=255)
                 gcgggcggGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGCACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                                                                <  2:2155561/150‑1 (MQ=255)
                              tCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAAt                                                                                                                   >  2:1655418/1‑150 (MQ=255)
                              tCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAAt                                                                                                                   >  1:986756/1‑150 (MQ=255)
                                  ggTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTAc                                                                                                               <  1:582303/150‑1 (MQ=255)
                                  ggTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTAc                                                                                                               <  2:2373563/150‑1 (MQ=255)
                                         ataCAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGaa                                                                                                        >  1:2450212/1‑150 (MQ=255)
                                             aaaCACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCtt                                                                                                    >  1:1413400/1‑150 (MQ=255)
                                             aaaCACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCtt                                                                                                    >  2:2177220/1‑150 (MQ=255)
                                             aaaCACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCtt                                                                                                    >  1:2965404/1‑150 (MQ=255)
                                                 acacCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAgc                                                                                                <  2:1510440/150‑1 (MQ=255)
                                                 acacCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAgc                                                                                                >  1:1416706/1‑150 (MQ=255)
                                                   cccTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAgcgc                                                                                              <  1:1421051/149‑1 (MQ=255)
                                                    ccTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGcc                                                                                             <  1:258180/150‑1 (MQ=255)
                                                     cTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGccc                                                                                            <  2:2450212/150‑1 (MQ=255)
                                                     cTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGccc                                                                                            >  1:2037329/1‑150 (MQ=255)
                                                        tACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATc                                                                                         <  1:1759621/150‑1 (MQ=255)
                                                          cTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCaa                                                                                       <  1:877326/150‑1 (MQ=255)
                                                           ttGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAg                                                                                      >  1:60050/1‑150 (MQ=255)
                                                             gtatgtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGtt                                                                                    <  2:2965404/150‑1 (MQ=255)
                                                             gtatgtaCAAGCATGGACAACCAAAATAATCGGACCTAAGCGTTCCTGGCGGGAGATTTCTGTTGGGGTATAGAACCACCTGTTCGAACACCACGCGCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTACAGACCCCATCAAGtt                                                                                    >  1:2994914/1‑150 (MQ=255)
                                                              tatctaCAAGCATCGCCAACCAAAATCATCGGAACTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTc                                                                                   <  2:1413400/150‑1 (MQ=39)
                                                               atgtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGGCCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCg                                                                                  >  2:645022/1‑150 (MQ=255)
                                                                tgtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCCCGCGCGGCGCTCg                                                                                                                                        <  2:1514211/96‑1 (MQ=255)
                                                                tgtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCCCGCGCGGCGCTCg                                                                                                                                        >  1:1514211/1‑96 (MQ=255)
                                                                 gtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGAt                                                                                >  2:1411080/1‑150 (MQ=255)
                                                                 gtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGAt                                                                                <  1:1655418/150‑1 (MQ=255)
                                                                   aCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCa                                                                              >  1:23291/1‑150 (MQ=255)
                                                                   aCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCa                                                                              <  1:1220375/150‑1 (MQ=255)
                                                                   aCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCa                                                                              >  1:1653280/1‑150 (MQ=255)
                                                                     aaGCATGGACAACCAAAATCATCGGACCTTATAGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAgg                                                                            <  2:2037329/150‑1 (MQ=255)
                                                                     aaGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATc                                                                                         >  2:2224449/1‑137 (MQ=255)
                                                                     aaGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATc                                                                                         <  1:2224449/137‑1 (MQ=255)
                                                                     aaGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAgg                                                                            <  2:60050/150‑1 (MQ=255)
                                                                     aaGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCCCGCGCGgcgc                                                                                                                                           <  2:316001/88‑1 (MQ=255)
                                                                     aaGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCCCGCGCGgcgc                                                                                                                                           >  1:316001/1‑88 (MQ=255)
                                                                       gCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATc                                                                               >  1:510229/1‑145 (MQ=255)
                                                                       gCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATc                                                                               <  2:510229/145‑1 (MQ=255)
                                                                               aaCCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGcc                                                                  >  2:642937/1‑150 (MQ=255)
                                                                                  cAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGAt                                                                                >  2:2469664/1‑133 (MQ=255)
                                                                                  cAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGAt                                                                                <  1:2469664/133‑1 (MQ=255)
                                                                                  cAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGgc                                                                           <  2:2622839/138‑1 (MQ=255)
                                                                                  cAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGgc                                                                           >  1:2622839/1‑138 (MQ=255)
                                                                                  cAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCt                                                               <  1:2177220/150‑1 (MQ=255)
                                                                                       taatCGGACCTTACCGTTCCTGGCTTTAGATTTATGTTGGGCTATAGAGCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGACGAACCTTCAGCGCCCCTCAAGTTCGATCAGGCAGCACCGCCCCTGcaca                                                          >  1:716021/1‑150 (MQ=39)
                                                                                                cTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCa                                                                                                  <  2:387568/101‑1 (MQ=255)
                                                                                                cTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCa                                                                                                  >  1:387568/1‑101 (MQ=255)
                                                                                                cTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAgg                                                                            <  2:2847661/123‑1 (MQ=255)
                                                                                                cTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAgg                                                                            >  1:2847661/1‑123 (MQ=255)
                                                                                                  tACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAg                                               >  1:891294/1‑150 (MQ=255)
                                                                                                     cGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCaa                                            >  2:2944628/1‑150 (MQ=255)
                                                                                                          cTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGcccc                                                                <  2:1211314/125‑1 (MQ=255)
                                                                                                          cTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGcccc                                                                >  1:1211314/1‑125 (MQ=255)
                                                                                                            ggCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATtcatca                                     <  1:1411080/150‑1 (MQ=255)
                                                                                                              cTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATtcatcact                                   <  2:1653280/150‑3 (MQ=255)
                                                                                                               ttGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCa                                                                                                  >  1:1557638/1‑86 (MQ=255)
                                                                                                               ttGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCa                                                                                                  <  2:1557638/86‑1 (MQ=255)
                                                                                                               ttGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCa                                                                              <  1:1102158/106‑1 (MQ=255)
                                                                                                               ttGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCa                                                                              >  2:1102158/1‑106 (MQ=255)
                                                                                                               ttGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTg                                  >  1:1036021/1‑150 (MQ=255)
                                                                                                                 gagaTTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCg                                                                                  >  1:133267/1‑100 (MQ=255)
                                                                                                                 gagaTTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCg                                                                                  <  2:133267/100‑1 (MQ=255)
                                                                                                                 gagaTTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGcccc                                                 >  1:2797494/1‑133 (MQ=255)
                                                                                                                 gagaTTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGcccc                                                 <  2:2797494/133‑1 (MQ=255)
                                                                                                                   gaTTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCtt                              <  2:891294/150‑1 (MQ=255)
                                                                                                                    aTTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCGGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTc                             >  1:878892/1‑150 (MQ=255)
                                                                                                                    aTTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTc                             >  2:2977553/1‑150 (MQ=255)
                                                                                                                    aTTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTc                             >  2:1936273/1‑150 (MQ=255)
                                                                                                                           ttGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGc                                              >  1:1562867/1‑126 (MQ=255)
                                                                                                                           ttGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGc                                              <  2:1562867/126‑1 (MQ=255)
                                                                                                                               ctgtgaGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGc                  >  2:527821/6‑150 (MQ=255)
                                                                                                                                 tataGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAg                <  2:2994914/150‑1 (MQ=255)
                                                                                                                                   taGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATc                                                                               <  2:160436/85‑1 (MQ=255)
                                                                                                                                   taGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATc                                                                               >  1:160436/1‑85 (MQ=255)
                                                                                                                                     gATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGa            >  2:1215088/1‑150 (MQ=255)
                                                                                                                                      aTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGc                                              >  1:742413/1‑115 (MQ=255)
                                                                                                                                      aTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGc                                              <  2:742413/115‑1 (MQ=255)
                                                                                                                                      aTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAg                <  1:3005361/145‑1 (MQ=255)
                                                                                                                                      aTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAg                >  2:3005361/1‑145 (MQ=255)
                                                                                                                                       tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGccc                                                                 >  1:2544133/1‑95 (MQ=255)
                                                                                                                                       tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGccc                                                                 <  2:2544133/95‑1 (MQ=255)
                                                                                                                                       tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGaac          <  1:1936273/150‑1 (MQ=255)
                                                                                                                                          aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAAt                                           <  1:136700/114‑1 (MQ=255)
                                                                                                                                          aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAAt                                           >  2:136700/1‑114 (MQ=255)
                                                                                                                                          aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCAAAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTc              <  2:1763435/143‑1 (MQ=255)
                                                                                                                                          aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCAAAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTc              >  1:1763435/1‑143 (MQ=255)
                                                                                                                                           gcgcCGCGCGGGCGGCGCCTATCTACACGGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAAcc      <  2:716021/150‑1 (MQ=255)
                                                                                                                                             gccgcGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAAcccc    >  1:2785176/1‑150 (MQ=255)
                                                                                                                                               cgcgcgGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAACCCCaa  >  1:441013/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GTCCTCTGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAACCCCAA  >  NC_002947/5737509‑5737801

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: