Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,450,173 C→G 5.8% A68G (GCC→GGC)  nicA → nicotinate dehydrogenase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,450,1730CG 5.8% 219.1 / 4.0 69A68G (GCC→GGC) nicAnicotinate dehydrogenase subunit A
Reads supporting (aligned to +/- strand):  ref base C (30/35);  new base G (2/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

CATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATC  >  NC_002947/4450024‑4450319
                                                                                                                                                     |                                                                                                                                                  
cATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGgcgc                                                                                                                                                    >  2:351091/1‑150 (MQ=255)
     cTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGccgccg                                                                                                                                               >  2:503708/1‑150 (MQ=255)
      tgccacCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCgccgg                                                                                                                                              <  1:1249512/147‑1 (MQ=255)
      tGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCgccgg                                                                                                                                              >  1:284796/1‑150 (MQ=255)
              cGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAata                                                                                                                                      <  1:895331/150‑1 (MQ=255)
               gctgctgctGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGGTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAatat                                                                                                                                     <  2:1402938/150‑1 (MQ=255)
                 tgctgctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATca                                                                                                                                   <  2:1174484/150‑1 (MQ=255)
                 tgctgctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATca                                                                                                                                   >  2:1325659/1‑150 (MQ=255)
                  gctgctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATcac                                                                                                                                  <  2:1369065/150‑1 (MQ=255)
                                  aaCGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaag                                                                                                                                         <  1:610291/127‑4 (MQ=255)
                                  aaCGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaag                                                                                                                                         >  2:610291/1‑124 (MQ=255)
                                  aaCGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTg                                                                                                                  >  1:766443/1‑150 (MQ=255)
                                  aaCGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTg                                                                                                                  <  1:351091/150‑1 (MQ=255)
                                  aaCGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTg                                                                                                                  >  1:1393841/1‑150 (MQ=255)
                                  aaCGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCACCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTg                                                                                                                  <  1:801802/150‑1 (MQ=255)
                                         tgtgCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgcagca                                                                                                           >  2:653614/1‑150 (MQ=255)
                                          gtgCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGaa                                                                                                                        <  1:758985/136‑1 (MQ=255)
                                          gtgCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGaa                                                                                                                        >  2:758985/1‑136 (MQ=255)
                                             ccTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTCGGCCGGCGGCgccggccaaggg                                                                                                                                       >  2:46469/1‑113 (MQ=255)
                                             ccTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTCGGCCGGCGGCgccggccaaggg                                                                                                                                       <  1:46469/118‑6 (MQ=255)
                                              cTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAg                                                                                                                       >  2:1056999/1‑133 (MQ=255)
                                              cTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAg                                                                                                                       <  1:1056999/133‑1 (MQ=255)
                                              cTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGcc                                                                                                      >  2:174971/1‑150 (MQ=255)
                                              cTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGcc                                                                                                      >  1:365610/1‑150 (MQ=255)
                                               tGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGgcgccg                                                                                                                                               <  2:550203/108‑1 (MQ=255)
                                               tGAACGGCCCCAAGTACGGCTGAGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGgcgccg                                                                                                                                               >  1:550203/1‑108 (MQ=255)
                                                  aCGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGGGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGccgcac                                                                                                  <  2:93220/150‑1 (MQ=255)
                                                     gCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgca                                                                                                              <  1:991186/135‑1 (MQ=255)
                                                     gCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgca                                                                                                              >  2:991186/1‑135 (MQ=255)
                                                       cccAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCcgcaccgcac                                                                                             >  2:1165533/1‑150 (MQ=255)
                                                        ccAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCgcaccgcacc                                                                                            >  1:202272/1‑150 (MQ=255)
                                                         cAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCAccc                                                                                           >  2:451410/1‑150 (MQ=255)
                                                          aaGTACGGCTGCGGCCTGGGCGAATGTGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCAcccc                                                                                          <  2:1189274/150‑1 (MQ=255)
                                                               cGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCa                                                                                     >  1:605762/1‑150 (MQ=255)
                                                                   tGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAg                                                                                 <  1:712346/150‑1 (MQ=255)
                                                                        ccTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGtt                                                                            >  1:755876/1‑150 (MQ=255)
                                                                         cTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTc                                                                           <  1:451410/150‑1 (MQ=255)
                                                                         cTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTc                                                                           >  2:733573/1‑150 (MQ=255)
                                                                          tGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGAGGTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCa                                                                          >  1:661334/1‑150 (MQ=255)
                                                                            ggCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATc                                                                        <  1:887050/150‑1 (MQ=255)
                                                                               gAATGCGGGGCGTGTACTGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATcgac                                                                     <  2:575871/150‑1 (MQ=255)
                                                                               gAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATcgac                                                                     <  1:653614/150‑1 (MQ=255)
                                                                                    cGGGGCGTGTACGGTGATCATCGACGGCGTGTCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCa                                                                <  1:1296/150‑1 (MQ=255)
                                                                                             tACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGcccc                                                       >  2:1220539/1‑149 (MQ=255)
                                                                                             tACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCa                                                       <  2:284796/150‑1 (MQ=255)
                                                                                               cGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAgt                                                     <  2:1227291/150‑1 (MQ=255)
                                                                                                ggTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGaa                                                                                                                        <  2:1084577/82‑1 (MQ=255)
                                                                                                ggTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGaa                                                                                                                        >  1:1084577/1‑82 (MQ=255)
                                                                                                  tGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGtgtg                                                  <  1:503708/150‑1 (MQ=255)
                                                                                                  tGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGtgtg                                                  >  1:502846/1‑150 (MQ=255)
                                                                                                    atcatcGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGt                                                <  2:1393841/150‑1 (MQ=255)
                                                                                                     tcatcGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGtt                                               <  1:706759/150‑1 (MQ=255)
                                                                                                        tagACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACt                                            >  1:429734/3‑150 (MQ=255)
                                                                                                         cGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTg                                           >  2:1011964/1‑150 (MQ=255)
                                                                                                           aCGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCa                                         <  2:365610/150‑1 (MQ=255)
                                                                                                            cGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTTCTCCAt                                        <  2:766443/150‑1 (MQ=255)
                                                                                                                gTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAAc                                    <  2:429734/150‑1 (MQ=255)
                                                                                                                gTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGACCCAGTGTGGTTACTGCATGAAc                                    <  2:546069/150‑1 (MQ=255)
                                                                                                                   gccgccCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGc                                 <  1:1165533/150‑1 (MQ=255)
                                                                                                                           cTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCAt                         <  1:875292/150‑1 (MQ=255)
                                                                                                                            tCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCACTATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCACGTGCAACAAGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATg                        >  2:159736/1‑150 (MQ=255)
                                                                                                                                  gTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATcgacga                                                                   >  1:1055541/1‑101 (MQ=255)
                                                                                                                                  gTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATcgacga                                                                   <  2:1055541/101‑1 (MQ=255)
                                                                                                                                  gTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGcc                  >  2:40805/1‑150 (MQ=255)
                                                                                                                                   tCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCa                 >  2:561169/1‑150 (MQ=255)
                                                                                                                                     aTCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGcatcatcat                         >  2:989031/1‑140 (MQ=255)
                                                                                                                                     aTCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCAt                         <  1:989031/140‑1 (MQ=255)
                                                                                                                                         ccTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGccc           <  1:733573/150‑1 (MQ=255)
                                                                                                                                         ccTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGccc           >  1:801979/1‑150 (MQ=255)
                                                                                                                                                cGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCAccc                                                                                           <  1:317644/63‑1 (MQ=255)
                                                                                                                                                cGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCAccc                                                                                           >  2:317644/1‑63 (MQ=255)
                                                                                                                                                 ggCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGAt   >  1:992804/1‑150 (MQ=255)
                                                                                                                                                 ggCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGAt   <  2:795565/150‑1 (MQ=255)
                                                                                                                                                  gcgcCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATc  <  1:1325659/150‑1 (MQ=255)
                                                                                                                                                  gcgcCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATc  >  1:90825/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATC  >  NC_002947/4450024‑4450319

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: