Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,693,502 T→C 5.7% V444A (GTG→GCG)  PP_4154 → class 3 aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,693,5020TC 5.7% 213.5 / 2.1 70V444A (GTG→GCG) PP_4154class 3 aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (37/29);  new base C (2/2);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

CAACTGGGAGCAACCCGGTTACGACCAGCATCTGCTGACTTCGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCGTGGAAAAAGGCC  >  NC_002947/4693353‑4693651
                                                                                                                                                     |                                                                                                                                                     
cAACTGGGAGCAACCCGGTTACGACCAGCATCTGCTGACTTCGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCAtgt                                                                                                                                                       <  1:1166795/150‑1 (MQ=255)
cAACTGGGAGCAACCCGGTTACGACCAGCATCTGCTGACTTCGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCAtgt                                                                                                                                                       <  1:1378255/150‑1 (MQ=255)
aaaCTGGGAGCAACCCGGTTACGACCAGCATCTGCTGACTTCGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCAtgt                                                                                                                                                       <  2:960287/149‑1 (MQ=255)
           aaCCCGGTTACGACCAGCATCTGCTGACTTCGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCa                                                                                                                                            <  2:928101/150‑1 (MQ=255)
             cccGGTTACGACCAGCATCTGCTGACTTCGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCACACCCTCGACGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGAATGTGCCCAACGGCAGc                                                                                                                                          >  1:798636/1‑150 (MQ=255)
               cGGTTACGACCAGCATCTGCTGACTTCGCAAGCCGCCCGCTCGTGCAGCCTGCACGGGCGGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGccc                                                                                                                                        <  1:1085219/150‑1 (MQ=255)
                         cAGCATCTGCTGACTTCGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAgcg                                                                                                                              >  2:575385/1‑150 (MQ=255)
                         cAGCATCTGCTGACTTCGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAgcg                                                                                                                              >  1:126637/1‑150 (MQ=255)
                          aGCATCTGCTGACTTCGCAAGCCGCCGGGTGCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAgcgc                                                                                                                             >  1:398304/1‑150 (MQ=255)
                                 gctgACTTCGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCTAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCgccggcagc                                                                                                                      >  1:82072/1‑150 (MQ=255)
                                   tgACTTCGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCgccggcagccg                                                                                                                    <  1:444742/150‑1 (MQ=255)
                                    gACTTCGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCgccggcagccgg                                                                                                                   >  1:1249495/1‑150 (MQ=255)
                                       ttCGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCgccggctgccggcgc                                                                                                                <  2:845061/150‑3 (MQ=255)
                                       ttCGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCAc                                                                                                                >  1:671420/1‑150 (MQ=255)
                                         cGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCgccggctgccggcgctt                                                                                                              <  1:575385/150‑5 (MQ=255)
                                          gCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTg                                                                                                             >  2:492249/1‑150 (MQ=255)
                                          gCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTg                                                                                                             >  2:1334202/1‑150 (MQ=255)
                                          gCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCATTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTg                                                                                                             >  1:19163/1‑150 (MQ=255)
                                             aGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGgca                                                                                                          <  2:495774/150‑1 (MQ=255)
                                                  ccccGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCAGgcagcccc                                                                                                     >  1:102498/1‑147 (MQ=255)
                                                  ccccGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGCGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGcaca                                                                                                     >  2:569003/1‑150 (MQ=255)
                                                  ccccGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGCGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGcaca                                                                                                     >  2:569197/1‑150 (MQ=255)
                                                       gcctgcagcctgcaCGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCTAGCTGCATGTGCCCAACGGCAGCCCGTTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAgcggc                                                                                                <  1:364585/150‑1 (MQ=255)
                                                       gcctgcagcctgcaCGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAtcggc                                                                                                <  2:798636/150‑1 (MQ=255)
                                                       gcctgcagcctgcaCGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGCGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAgcggc                                                                                                <  1:569197/150‑1 (MQ=255)
                                                       gcctgcagcctgcaCGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGCGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAgcggc                                                                                                <  1:569003/150‑1 (MQ=255)
                                                              gcctgcaCGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGcc                                                                                         <  2:1249495/150‑1 (MQ=255)
                                                                ctgcaCGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACg                                                                                                                                               <  2:826358/94‑1 (MQ=255)
                                                                ctgcaCGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACg                                                                                                                                               >  1:826358/1‑94 (MQ=255)
                                                                 tgcaCGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCggg                                                                                      <  2:82072/150‑1 (MQ=255)
                                                                 tgcaCGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCggg                                                                                      <  1:1126372/150‑1 (MQ=255)
                                                                                     tCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAg                                                                                                                                <  1:34823/88‑1 (MQ=255)
                                                                                     tCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAg                                                                                                                                >  2:34823/1‑88 (MQ=255)
                                                                                     tCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCgccggca                                                                                                                        <  2:1387996/96‑1 (MQ=255)
                                                                                     tCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCgccggca                                                                                                                        >  1:1387996/1‑96 (MQ=255)
                                                                                      cGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGGCTGCGCACCCCGCAtt                                                                 <  2:847634/150‑1 (MQ=255)
                                                                                        cAAACCCATATCGATCGCCCCGAGGAGCACGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGACGGACGGCTGGACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTgg                                                               <  2:126637/150‑1 (MQ=255)
                                                                                         aaaCCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTggg                                                              <  2:1157093/150‑1 (MQ=255)
                                                                                             ccATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCttt                                                          >  1:213333/1‑150 (MQ=255)
                                                                                             ccATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCttt                                                          >  1:1042175/1‑150 (MQ=255)
                                                                                                 atcgatcgCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTgg                                                      >  2:362777/1‑150 (MQ=255)
                                                                                                 atcgatcgCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTgg                                                      <  2:213333/150‑1 (MQ=255)
                                                                                                    gatcgCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTg                                                   >  2:977160/1‑150 (MQ=255)
                                                                                                      tcgCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCa                                                 >  1:1111741/1‑150 (MQ=255)
                                                                                                       cgCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCaa                                                <  2:671420/150‑1 (MQ=255)
                                                                                                       cgCCCCGAGGAGCCCGCCACCTGCGCCCTTTGCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCaa                                                <  1:1374232/150‑1 (MQ=255)
                                                                                                         ccccGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAgg                                              >  1:896912/1‑150 (MQ=255)
                                                                                                            cGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCt                                           <  2:1042175/150‑1 (MQ=255)
                                                                                                              aggagCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGccatt                                         <  1:1119705/150‑4 (MQ=255)
                                                                                                                   cccGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGAc                                    >  1:245118/1‑150 (MQ=255)
                                                                                                                   cccGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGCGCAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGAc                                    >  1:1244556/1‑150 (MQ=255)
                                                                                                                    ccGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACg                                   >  1:587481/1‑150 (MQ=255)
                                                                                                                     cGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGc                                  >  2:575380/1‑150 (MQ=255)
                                                                                                                           cTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGc                            >  2:672607/1‑150 (MQ=255)
                                                                                                                            tGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTAGCTGCAAGGCCTGGAAGACGCACCGGCa                           >  2:998367/1‑150 (MQ=255)
                                                                                                                             gcgcCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAg                          >  2:160332/1‑150 (MQ=255)
                                                                                                                              cgcCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGa                         <  1:977160/150‑1 (MQ=255)
                                                                                                                                  cTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCgccg                                                                                                                           <  2:647810/48‑1 (MQ=255)
                                                                                                                                  cTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCgccg                                                                                                                           >  1:647810/1‑48 (MQ=255)
                                                                                                                                   ttttCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGAAGGCACACGTCGGCTGCGCACCCTGCATTTGGCTTTGTGGCTGCAAGGCTTGGAAGACGCATCGGCAGACGGcc                    >  1:386288/1‑150 (MQ=255)
                                                                                                                                        gtcgAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGcc               >  1:1154790/1‑150 (MQ=255)
                                                                                                                                        gtcgAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGcc               >  1:730121/1‑150 (MQ=255)
                                                                                                                                         tcgAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCgggg                                                                                     >  1:174864/1‑79 (MQ=255)
                                                                                                                                         tcgAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCgggg                                                                                     <  2:174864/79‑1 (MQ=255)
                                                                                                                                         tcgAGCTGCATGTGCACAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCg              >  2:1195716/1‑150 (MQ=255)
                                                                                                                                          cgAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCGt             <  1:1334202/150‑1 (MQ=255)
                                                                                                                                           gAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTg                                          <  1:1145890/120‑1 (MQ=255)
                                                                                                                                           gAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTg                                          >  2:1145890/1‑120 (MQ=255)
                                                                                                                                                 cATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCGTGGaaaaa      <  1:216895/150‑1 (MQ=255)
                                                                                                                                                  aTGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTtgt                                                        >  1:474857/1‑99 (MQ=255)
                                                                                                                                                  aTGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTtgt                                                        <  2:474857/99‑1 (MQ=255)
                                                                                                                                                  aTGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCGTGGAAAAag     >  1:976485/1‑150 (MQ=255)
                                                                                                                                                  aTGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCGTGGAAAAag     >  1:978591/1‑150 (MQ=255)
                                                                                                                                                     tgCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGaa                                       >  1:1359492/1‑113 (MQ=255)
                                                                                                                                                     tgCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGaa                                       <  2:1359492/113‑1 (MQ=255)
                                                                                                                                                     tgCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGAACCGGAAGAAGGCCAGGCGCTGGaagcagtcc  >  2:305552/1‑143 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAACTGGGAGCAACCCGGTTACGACCAGCATCTGCTGACTTCGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCGTGGAAAAAGGCC  >  NC_002947/4693353‑4693651

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