Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,963,902 T→G 9.6% Y216D (TAT→GAT)  PP_5229 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,963,9020TG 9.6% 208.1 / 12.8 83Y216D (TAT→GAT) PP_5229hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/36);  new base G (3/5);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

GGCGTTCCTGTTCGGTGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGCTACCG  >  NC_002947/5963754‑5964046
                                                                                                                                                    |                                                                                                                                                
ggCGTTCCTGTTCGGTGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTAAAATGGCTGCGTCGGTACCCTTTTCCTCGGCta                                                                                                                                                 <  1:815533/150‑1 (MQ=255)
   gTTCCTGTTCGGTGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTatat                                                                                                                                              >  2:238190/1‑150 (MQ=255)
       cTGTTCGGTGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGcc                                                                                                                                          <  2:1008030/150‑1 (MQ=255)
         gTTCGGTGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGa                                                                                                                                        >  2:85018/1‑150 (MQ=255)
             ggTGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTg                                                                                                                                    <  1:756900/150‑1 (MQ=255)
             ggTGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCGATATAGCCGAggaa                                                                                                                                    >  2:270995/1‑148 (MQ=255)
               tGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCt                                                                                                                                  >  1:1119147/1‑150 (MQ=255)
               tGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCGATATAGCCGAggaaaa                                                                                                                                  <  1:1546108/150‑5 (MQ=255)
                gaagaaCAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGGAGCTCTTGAGTACCAACGGCTAAATGGGGGTCTGGACACCGGCAGCCGCGACACGCTACACCACACGATCAGCATCGGTACCCTTTTCCTCACATACATCGACGAgaccctt                                                                                                                                 >  2:679599/1‑144 (MQ=255)
                gaagaaCAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCGATATAGCCGAggaaaag                                                                                                                                 <  1:270995/150‑6 (MQ=255)
                gaagaaCAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCGATATAGCCGAggaaaag                                                                                                                                 <  1:1494092/150‑6 (MQ=255)
                 aagaaCAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTg                                                                                                                                >  1:587774/1‑150 (MQ=255)
                    aaCAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGcc                                                                                                                             <  1:1734284/150‑1 (MQ=255)
                        gcgcCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGt                                                                                                                         >  2:23518/1‑150 (MQ=255)
                        gcgcCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGt                                                                                                                         <  1:1092167/150‑1 (MQ=255)
                           ccGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTtgt                                                                                                                      <  2:587774/150‑1 (MQ=255)
                                    aGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGt                                                                                                             <  1:822484/150‑1 (MQ=255)
                                     ggCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTc                                                                                                            >  2:1504269/1‑150 (MQ=255)
                                     ggCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTc                                                                                                            >  2:1169435/1‑150 (MQ=255)
                                        tCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGATATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCAcc                                                                                                         <  2:716854/150‑1 (MQ=255)
                                            gcgcCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGCGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGa                                                                                                     >  1:1193317/1‑150 (MQ=255)
                                               ccGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTGTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGAccc                                                                                                  <  1:372678/150‑1 (MQ=255)
                                                 gTGGCAGGTCTTGAGTATCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTTGCCGCGTTGTGCCACGCGTCACCCAGACCCTg                                                                                                >  2:776843/1‑150 (MQ=255)
                                                    gCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGACACGCGTGAAACAGACGCTgcgc                                                                                             >  2:585662/1‑150 (MQ=255)
                                                    gCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCGATATAGCCGAg                                                                                                                                       >  2:778983/1‑108 (MQ=255)
                                                    gCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCGATATAGCCGAg                                                                                                                                       <  1:778983/108‑1 (MQ=255)
                                                      aGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTGCCTCGGTTAGAGCGCCTAGGGGCGTGGGCGGTGTGTTCCACGCGTCACCCAGACCCTGCGgcc                                                                                           >  2:398549/1‑150 (MQ=255)
                                                           ttGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTAc                                                                                      >  2:857592/1‑150 (MQ=255)
                                                             gAGTACCTGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGc                                                                                    >  2:297843/1‑150 (MQ=255)
                                                             gAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACCAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCCGCGCCCGGTACGc                                                                                    <  1:711545/150‑1 (MQ=255)
                                                             gAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCCGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGc                                                                                    <  2:380026/150‑1 (MQ=255)
                                                              ggTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGCGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCt                                                                                   <  2:1056033/149‑1 (MQ=255)
                                                                tACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTc                                                                                                            <  2:1806670/123‑1 (MQ=255)
                                                                tACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTc                                                                                                            >  1:1806670/1‑123 (MQ=255)
                                                                tACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCtga                                                                                 <  2:766289/150‑1 (MQ=255)
                                                                tACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCtga                                                                                 <  1:85018/150‑1 (MQ=255)
                                                                 aCCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTtgt                                                                                                                      >  2:197401/1‑112 (MQ=255)
                                                                 aCCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTtgt                                                                                                                      <  1:197401/112‑1 (MQ=255)
                                                                 aCCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGCGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCtgat                                                                                <  2:1821728/150‑1 (MQ=255)
                                                                 aCCAGGGGCTGCATGGGGTTCTGGCCATCGGAAGACGAGAACAAAAACACTACGATAGCAGCGTCGGTACACTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCtgat                                                                                <  1:425066/150‑1 (MQ=255)
                                                                     ggggCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCAAGACCCTGCGCCCGGTACGCTGATGaaa                                                                            >  1:475544/1‑150 (MQ=255)
                                                                        gCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGt                                                                                                             <  1:1320520/114‑1 (MQ=255)
                                                                        gCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGt                                                                                                             >  2:1320520/1‑114 (MQ=255)
                                                                             aTGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCt                                                                    <  2:23944/150‑1 (MQ=255)
                                                                             aTGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCt                                                                    >  1:23944/1‑150 (MQ=255)
                                                                                  gttttGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCAc                                                                                                                 <  1:503668/96‑1 (MQ=255)
                                                                                  gTTTTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCAc                                                                                                                 >  2:503668/1‑100 (MQ=255)
                                                                                      tGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCg                                                                                          <  2:47487/119‑1 (MQ=255)
                                                                                      tGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCg                                                                                          >  1:47487/1‑119 (MQ=255)
                                                                                      tGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAAc                                                                           >  1:1584501/1‑134 (MQ=255)
                                                                                      tGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAAc                                                                           <  2:1584501/134‑1 (MQ=255)
                                                                                       ggCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAg                                                                      <  1:1151267/138‑1 (MQ=255)
                                                                                       ggCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAg                                                                      >  2:1151267/1‑138 (MQ=255)
                                                                                        gCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTAtt                                                         >  2:1380883/1‑150 (MQ=255)
                                                                                        gCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTAtt                                                         <  1:459619/150‑1 (MQ=255)
                                                                                          cATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTgc                                                       >  2:1650303/1‑150 (MQ=255)
                                                                                                   ccGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCAc                                                                                                                 >  1:494893/1‑83 (MQ=255)
                                                                                                   ccGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCAc                                                                                                                 <  2:494893/83‑1 (MQ=255)
                                                                                                    cgcgACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTgc                                             >  1:75089/1‑150 (MQ=255)
                                                                                                            gCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTgg                                                                  <  2:386593/121‑1 (MQ=255)
                                                                                                            gCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTgg                                                                  >  1:386593/1‑121 (MQ=255)
                                                                                                            gCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAg                                     >  1:1267153/1‑150 (MQ=255)
                                                                                                             cAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGTTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAgc                                    <  1:398549/150‑1 (MQ=255)
                                                                                                             cAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCAAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAgc                                    <  2:475544/150‑1 (MQ=255)
                                                                                                               acacTACAAGAGACGAGGCGTTACCCTTTTCCTCGGCTATATCGCCTAGTTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAgcgc                                  >  1:1499551/1‑150 (MQ=255)
                                                                                                                 acTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCa                                >  2:1729395/1‑150 (MQ=255)
                                                                                                                    aCAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGgt                             >  2:1474231/1‑150 (MQ=255)
                                                                                                                    aCAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCCGTCCGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGgt                             >  2:57164/1‑150 (MQ=255)
                                                                                                                     cAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATAACGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGgtg                            >  1:627319/1‑150 (MQ=255)
                                                                                                                     cAAGAGCAGCGTCGGTACCCTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGAAGAAACGTCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGgtgt                           <  1:987485/150‑1 (MQ=255)
                                                                                                                      aaGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGgtgt                           >  1:823998/1‑150 (MQ=255)
                                                                                                                       agagCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTc                          <  1:1380883/150‑1 (MQ=255)
                                                                                                                           cagcGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAg                      <  1:23518/150‑1 (MQ=255)
                                                                                                                              cGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGc                                                                         <  2:308883/96‑1 (MQ=255)
                                                                                                                              cGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGc                                                                         >  1:308883/1‑96 (MQ=255)
                                                                                                                                tCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCacac                 <  2:823998/150‑1 (MQ=255)
                                                                                                                                        aagcaccTCGGCGATATAGCCGAGGTGCTTGGCCGCGTTGTGCCACGTGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAgc                                    >  1:1249497/6‑123 (MQ=255)
                                                                                                                                        aagcaccTCGGCGATATAGCCGAGGTGCTTGGCCGCGTTGTGCCACGTGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAgc                                    <  2:1249497/118‑1 (MQ=255)
                                                                                                                                         ttttCCTCGGCTATATCGCCGTGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAggg        <  1:81967/150‑1 (MQ=255)
                                                                                                                                             ccTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGgtg                            <  2:1490339/126‑1 (MQ=255)
                                                                                                                                             ccTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGgtg                            >  1:1490339/1‑126 (MQ=255)
                                                                                                                                               tCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAgcgc                                  <  1:204686/118‑1 (MQ=255)
                                                                                                                                               tCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAgcgc                                  >  2:204686/1‑118 (MQ=255)
                                                                                                                                               tCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGCTACCg  >  1:1601786/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GGCGTTCCTGTTCGGTGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGCTACCG  >  NC_002947/5963754‑5964046

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: