Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 664,889 C→T 5.4% A584A (GCC→GCT)  yfhM → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947664,8890CT 5.4% 211.6 / 2.3 74A584A (GCC→GCT) yfhMlipoprotein
Reads supporting (aligned to +/- strand):  ref base C (34/36);  new base T (2/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

TACCAATTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCGA  >  NC_002947/664743‑665038
                                                                                                                                                  |                                                                                                                                                     
tACCAATTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGc                                                                                                                                                    <  2:234208/150‑1 (MQ=255)
caCCGATTCGGCTCGGGCACCGAAACCGAGCTGGACCAGGACCTGGAACGGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATAGCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGc                                                                                                                                                    <  1:742461/149‑1 (MQ=255)
       tCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGCCCCGGCCGGATGGGCGGCGGTCGATAGCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCTGGCCTGCacg                                                                                                                                             <  1:1671106/150‑3 (MQ=255)
        tggCTCGTTCACTGACATCTTTCTCCACCCGGCCCTGGAATTGTTCTCAGTCTCCCTTTTCCAGGCGGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGa                                                                                                                                            <  1:841424/149‑1 (MQ=255)
           cTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTc                                                                                                                                         >  2:836022/1‑150 (MQ=255)
             cGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCATGAGTcag                                                                                                                                       <  2:806145/150‑3 (MQ=255)
               gtaaCCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGACCAGCCTGCAGGAGTCGGGt                                                                                                                                     >  1:212958/4‑150 (MQ=255)
                    cGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCGGCCTGCAGGAGTCGGGTggccg                                                                                                                                >  1:123526/1‑150 (MQ=255)
                     ggaaCCGGTCAGAGCCGGGACGCGGAACTGGACGAAGTCACCCTTGCCCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTggccgg                                                                                                                               <  1:1238399/148‑1 (MQ=255)
                          cgaagtGACCCAGGACCTGGACCTGGCCGATCTCCATGATGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGACGGCCGAt                                                                                                                          <  1:1841460/145‑1 (MQ=255)
                           gATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATc                                                                                                                         >  2:1844806/1‑150 (MQ=255)
                           gATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATc                                                                                                                         >  2:1103918/1‑150 (MQ=255)
                               ctacccAGGCCCTGGAACTTGTCGAGGTCGCCGTCGACCGGGCCGGGAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCAccc                                                                                                                     <  1:674736/146‑1 (MQ=255)
                                 gaCCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCgg                                                                                                                   >  1:1749423/2‑150 (MQ=255)
                                      ggACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGCCCGGCAATTAACTCTGCAACTGACCGTGCAGGCCAGCCCGCAGGAGTCGGGTGGCCGGCCGAGCACCCGGCGCCt                                                                                                              >  1:329627/1‑150 (MQ=255)
                                         ccTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGa                                                                                                           >  1:1198369/1‑150 (MQ=255)
                                          cTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGag                                                                                                          >  2:1840153/1‑150 (MQ=255)
                                          cTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGag                                                                                                          >  2:1780387/1‑150 (MQ=255)
                                           tGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGagc                                                                                                         <  1:1372396/150‑1 (MQ=255)
                                           tGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGagc                                                                                                         >  2:764053/1‑150 (MQ=255)
                                            ggAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGagca                                                                                                        >  1:1463580/1‑150 (MQ=255)
                                              aaCTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGagcagc                                                                                                      <  2:289653/150‑1 (MQ=255)
                                               aCTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGGAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGCTCAACCACCCCCTAGTGCAGcc                                                                                                     >  2:1050195/1‑150 (MQ=255)
                                                 tGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTa                                                                                                   >  1:2756/1‑150 (MQ=255)
                                                  ggaggaCGTCACCCTTGACCAGGCGGGAAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGagcaccatag                                                                                                  <  2:442186/150‑7 (MQ=255)
                                                   gACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATc                                                                                                 >  1:488549/1‑150 (MQ=255)
                                                   gACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGCAGACGTGAGCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATc                                                                                                 <  2:1589901/150‑1 (MQ=255)
                                                          tCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCAGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCgg                                                                                          >  1:1760476/1‑150 (MQ=255)
                                                              ccTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGa                                                                                      >  1:1034304/1‑150 (MQ=255)
                                                                ttcaCCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAgc                                                                                    <  1:1161491/147‑1 (MQ=255)
                                                                ttGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAgc                                                                                    >  1:1121629/1‑150 (MQ=255)
                                                                      aGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGCGCGCCTGc                                                                              >  1:790012/1‑150 (MQ=255)
                                                                       ggCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGcc                                                                             >  2:1867265/1‑150 (MQ=255)
                                                                         ccGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCgg                                                                           <  2:2756/150‑1 (MQ=255)
                                                                           ggCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGc                                                                         <  2:1640189/150‑1 (MQ=255)
                                                                           ggCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGc                                                                         <  2:488549/150‑1 (MQ=255)
                                                                             cAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTTAGCGCCTGGCGGGCCt                                                                       >  2:617886/1‑150 (MQ=255)
                                                                              aaGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTg                                                                      <  1:334647/150‑1 (MQ=255)
                                                                                 gCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGa                                                                   <  1:1844806/150‑1 (MQ=255)
                                                                                 gCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGa                                                                   <  1:1333366/150‑1 (MQ=255)
                                                                                       gtCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTg                                                             >  1:1798695/1‑150 (MQ=255)
                                                                                        tCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGt                                                            <  1:764053/150‑1 (MQ=255)
                                                                                             atCGAAAGCCGTTGGGCCGAAGCCCGCGCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGaa                                                       <  1:392420/150‑1 (MQ=255)
                                                                                                  aaGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGCGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTga                                                  <  2:1503396/150‑1 (MQ=255)
                                                                                                   aGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTgag                                                 <  1:1050195/150‑1 (MQ=255)
                                                                                                     ccGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTgag                                                 <  1:1105129/148‑1 (MQ=255)
                                                                                                     ccGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTgag                                                 >  2:1105129/1‑148 (MQ=255)
                                                                                                     ccGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTCGCCACCGGCCGTGCCTGCGAGGCCTGTTCGAAGGTgagga                                               >  1:1040010/1‑150 (MQ=255)
                                                                                                      cGTTGGGCCGAAGCCCGCTCTCCGCGGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTgaggag                                              <  2:348908/150‑1 (MQ=255)
                                                                                                             ccGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCAGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATa                                       >  1:1314608/1‑150 (MQ=255)
                                                                                                                 ggCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGa                                   <  1:1747409/149‑1 (MQ=255)
                                                                                                                   cccGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACg                                 >  2:1863967/1‑150 (MQ=255)
                                                                                                                    ccGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGAc                                            >  2:1272820/1‑138 (MQ=255)
                                                                                                                    ccGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGAc                                            <  1:1272820/138‑1 (MQ=255)
                                                                                                                    ccGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACgg                                >  2:271694/1‑150 (MQ=255)
                                                                                                                     cGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCTCGTGTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACggg                               >  1:797638/1‑150 (MQ=255)
                                                                                                                          tcCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGt                          <  1:1840153/150‑1 (MQ=255)
                                                                                                                           ccGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTg                         <  1:1780387/150‑1 (MQ=255)
                                                                                                                           ccGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCAGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTg                         <  2:1760476/150‑1 (MQ=255)
                                                                                                                            cGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTgg                        >  1:1830229/1‑150 (MQ=255)
                                                                                                                              cTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGaa                      <  2:80658/150‑1 (MQ=255)
                                                                                                                               tGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAAt                     >  1:937490/1‑150 (MQ=255)
                                                                                                                                 cAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCTATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTc                   <  2:123526/150‑1 (MQ=255)
                                                                                                                                  aaCTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCg                  >  1:194735/1‑150 (MQ=255)
                                                                                                                                   aCTGCCCGGGCAGGCCAGCCTGCAGGAGTCGGGGGGCCGGCCGATCACCCGGCGCCGGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGa                 <  1:286879/150‑1 (MQ=255)
                                                                                                                                   aCTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGa                 <  2:790012/150‑1 (MQ=255)
                                                                                                                                     tGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCgg                                                                           >  1:835956/1‑90 (MQ=255)
                                                                                                                                     tGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCgg                                                                           <  2:835956/90‑1 (MQ=255)
                                                                                                                                        tcctGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCGGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTcc            >  1:525552/4‑150 (MQ=255)
                                                                                                                                           tGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCgg                                                                           <  2:419035/84‑1 (MQ=255)
                                                                                                                                           tGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCgg                                                                           >  1:419035/1‑84 (MQ=255)
                                                                                                                                           tGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCACGCCGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCtgg         <  1:1367468/150‑1 (MQ=255)
                                                                                                                                             cAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCtggtt       <  2:1749423/150‑2 (MQ=255)
                                                                                                                                               ggCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGCGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGc     >  1:1009625/1‑150 (MQ=255)
                                                                                                                                                  cAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCGa  >  2:62328/1‑150 (MQ=255)
                                                                                                                                                  cAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCATGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCAGTGGCGTACGAGTTCCTGGTGGCCGa  >  2:585641/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TACCAATTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCGA  >  NC_002947/664743‑665038

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