Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 575,944 C→T 5.3% Y327Y (TAC→TAT fdoG → formate dehydrogenase‑O major subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947575,9440CT 5.3% 229.6 / 2.1 75Y327Y (TAC→TATfdoGformate dehydrogenase‑O major subunit
Reads supporting (aligned to +/- strand):  ref base C (39/32);  new base T (2/2);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

ACCGACATCGCCTTCATGGGCGGGCTGATCAACTACCTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTCCAGCTGATGAAGCAGCATTACAGCCGCTATACCCCGG  >  NC_002947/575801‑576089
                                                                                                                                               |                                                                                                                                                 
accgACATCGCCTTCATGGGCGGGCTGATCAACTAGCTGCTCAGCAACGACAATATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTGATCGTCATAGAGAACTATGGCTTCGAGGATGGCCTTTTCAGCGGCTacgtcgt                                                                                                                                              <  6:153215/149‑2 (MQ=255)
 ccgACATCGCCTTCATGGGCGGGCTGATCAACTACCTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAg                                                                                                                                             <  4:320952/149‑1 (MQ=255)
       tCGCCTTCATGGGCGGGCTGATCAACTACCTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAg                                                                                                                                       >  1:212485/1‑149 (MQ=255)
           cTTCATGGGCGGGCTGATCAACTACCTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCgtg                                                                                                                                   >  7:334647/1‑149 (MQ=255)
              cATGGGCGGGCTGATCAACTACCTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCgtgtg                                                                                                                                 >  4:219368/1‑148 (MQ=255)
                   gcgggcTGATCAACTACCTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGc                                                                                                                            <  3:39473/148‑1 (MQ=255)
                      ggcTGATCAACTACCTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCaa                                                                                                                                        >  3:315241/1‑133 (MQ=255)
                      ggcTGATCAACTACCTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCaa                                                                                                                                        <  4:315241/133‑1 (MQ=255)
                      ggcTGATCAACTACCTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGAc                                                                                                                        <  8:270474/149‑1 (MQ=255)
                         tGATCAACTACCTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTAAGAGGCTACGACGATGCCAAGCGTGTGTACGCCGACAAg                                                                                                                     >  8:166984/1‑149 (MQ=255)
                           aTCAACTACCTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGCCGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTc                                                                                                                   >  2:232413/1‑149 (MQ=255)
                              aaCTACCTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCgg                                                                                                                >  7:95972/1‑149 (MQ=255)
                                cTACCTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCt                                                                                                              <  5:262655/149‑1 (MQ=255)
                                 tACCTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTg                                                                                                             >  4:264433/1‑149 (MQ=255)
                                  aCCTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTgg                                                                                                            >  5:137478/1‑149 (MQ=255)
                                  aCCTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTgg                                                                                                            <  3:219368/149‑1 (MQ=255)
                                   ccTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTggg                                                                                                           <  1:198389/149‑1 (MQ=255)
                                   ccTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTggg                                                                                                           <  3:264433/149‑1 (MQ=255)
                                      gCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGt                                                                                                         >  8:298285/1‑148 (MQ=255)
                                       cTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACCCTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGtt                                                                                                        <  7:166984/148‑1 (MQ=255)
                                       cTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTa                                                                                                       >  8:247481/1‑149 (MQ=255)
                                          aGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGa                                                                                                    <  6:137478/149‑1 (MQ=255)
                                              aCGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTcc                                                                                                                  >  6:8077/1‑131 (MQ=255)
                                              aCGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTcc                                                                                                                  <  5:8077/131‑1 (MQ=255)
                                               cGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTgg                                                                                               <  2:212485/149‑1 (MQ=255)
                                                gACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTATGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTgg                                                                                               >  3:321949/1‑148 (MQ=255)
                                                gACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTATGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTgg                                                                                               <  4:321949/148‑1 (MQ=255)
                                                 aCAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGt                                                                                             >  4:314124/1‑149 (MQ=255)
                                                  cAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGt                                                                                             >  2:307278/1‑148 (MQ=255)
                                                   aaGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGt                                                                                             >  1:248650/1‑147 (MQ=255)
                                                      aTCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGa                                                                                        >  6:186936/1‑149 (MQ=255)
                                                         cAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCgt                                                                                                                                    >  8:101571/1‑102 (MQ=255)
                                                         cAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCgt                                                                                                                                    <  7:101571/102‑1 (MQ=255)
                                                          aGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGc                                                                                    >  1:307912/1‑149 (MQ=255)
                                                           gCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCt                                                                                   <  7:298285/149‑1 (MQ=255)
                                                           gCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCt                                                                                   >  8:281486/1‑149 (MQ=255)
                                                           gCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCATCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGc                                                                                    >  1:180949/1‑148 (MQ=255)
                                                            cACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCgtgtgt                                                                                                                                <  6:138404/103‑1 (MQ=255)
                                                             aCGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTAt                                                                                 >  6:77989/1‑149 (MQ=255)
                                                               gAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGc                                                                               >  5:274616/1‑149 (MQ=255)
                                                                    cGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAgg                                                                          <  7:247481/149‑1 (MQ=255)
                                                                     gTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGt                                                                         >  2:280198/1‑149 (MQ=255)
                                                                     gTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGt                                                                         <  8:334647/149‑1 (MQ=255)
                                                                       gcgcAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTc                                                                        >  5:53850/1‑148 (MQ=255)
                                                                        cgcAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGa                                                                      <  8:95972/149‑1 (MQ=255)
                                                                         gcAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGAc                                                                     <  3:314124/149‑1 (MQ=255)
                                                                           aaCTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGAccc                                                                   <  1:232413/149‑1 (MQ=255)
                                                                                cacTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGacac                                                              >  7:96138/1‑149 (MQ=255)
                                                                                cacTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGaca                                                               <  5:186936/148‑1 (MQ=255)
                                                                                     aCGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTg                                                                                                <  5:64946/110‑1 (MQ=255)
                                                                                     aCGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTg                                                                                                >  6:64946/1‑110 (MQ=255)
                                                                                                aTCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTg                                              >  1:334425/1‑149 (MQ=255)
                                                                                                 tCGTCAAAGATAACTTTGGCTTCGAGGATGGCCTGTTCAGCGGCTACGCCGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGc                                             <  2:334425/149‑1 (MQ=255)
                                                                                                    tCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCgtg                                          >  6:335346/1‑149 (MQ=255)
                                                                                                     cAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCgt                                           >  8:232646/1‑147 (MQ=255)
                                                                                                     cAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCgt                                           <  7:232646/147‑1 (MQ=255)
                                                                                                       aaGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTc                                       <  6:274616/149‑1 (MQ=255)
                                                                                                         gagaACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTc                                       >  7:300201/1‑147 (MQ=255)
                                                                                                               tATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGAc                                                                                                                        <  6:120091/60‑1 (MQ=255)
                                                                                                               tATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGAc                                                                                                                        >  5:120091/1‑60 (MQ=255)
                                                                                                                 tGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTCCAGCtgatg                              >  2:84582/1‑148 (MQ=255)
                                                                                                                 tGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTCCAGCtgatg                              >  2:268987/1‑148 (MQ=255)
                                                                                                                   gCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTCCAGCTGATGAag                           <  1:307278/149‑1 (MQ=255)
                                                                                                                        gAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTCCAGCTGATGAagcagc                       <  1:280198/148‑1 (MQ=255)
                                                                                                                          ggATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTCCAGCTGATGAAGCAGCAtt                    <  7:281486/149‑1 (MQ=255)
                                                                                                                          ggATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTCCAGCTGATGAAGCAGCAt                     >  5:288181/1‑148 (MQ=255)
                                                                                                                               gCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTCCAGCTGATGAAGCAGCATTACAg                <  5:77989/148‑1 (MQ=255)
                                                                                                                                 cTGTTCAGCGGCTATGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTCCAGCTGATGAAGCAGCATTACAGCCg             >  3:162302/1‑149 (MQ=255)
                                                                                                                                      cAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTCCAGCTGATGAAGCAGCATTACAGCCGCtata        >  6:352877/1‑149 (MQ=255)
                                                                                                                                       aGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTCCAGCTGATGAAGCAGCATTACAGcc              >  5:189403/1‑142 (MQ=255)
                                                                                                                                       aGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTCCAGCTGATGAAGCAGCATTACAGcc              <  6:189403/142‑1 (MQ=255)
                                                                                                                                        gCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTCCAGCTGATGAAGCAGCATTACAGCCGCTATAcc      >  2:149493/1‑149 (MQ=255)
                                                                                                                                         cGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTCCAGCTGATGAAGCAGCATTACAGCCGCTATAccc     >  4:210274/1‑149 (MQ=255)
                                                                                                                                           gCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTCCAGCTGATGAAGCAGCATTACAGCCGCTATACCCCg   <  3:210274/149‑1 (MQ=255)
                                                                                                                                            cTATGACGAGTCCATGCGTGTGTAAGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGTCCCGACACTGCAGCACCCTCGCTGCGTGTTCCAGCTGATGAAGCAGCATTACAGCCGCTATACCCCg   <  4:162302/148‑1 (MQ=255)
                                                                                                                                             tACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTCCAGCTGATGAAGCAGCATTACAGCCGCTATACCCCgg  >  2:166496/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
ACCGACATCGCCTTCATGGGCGGGCTGATCAACTACCTGCTCAGCAACGACAAGATCCAGCACGAGTACGTGCGCAACTACACTGACGTGTCGTTCATCGTCAAAGAGAACTATGGCTTCGAGGATGGCCTGTTCAGCGGCTACGACGAGGCCAAGCGTGTGTACGCCGACAAGTCCGGCTGGGGTTACGAGCTGGGTGAGGACGGCTATGCCAAGGTCGACCCGACACTGCAGCACCCTCGCTGCGTGTTCCAGCTGATGAAGCAGCATTACAGCCGCTATACCCCGG  >  NC_002947/575801‑576089

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