Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,231,466 G→C 6.1% intergenic (+227/‑199) ygjK → / → fadH alpha‑glucosidase/2,4‑dienoyl‑CoA reductase, NADH and FMN‑linked

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,231,4660GC 6.1% 344.5 / 11.7 99intergenic (+227/‑199)ygjK/fadHalpha‑glucosidase/2,4‑dienoyl‑CoA reductase, NADH and FMN‑linked
Reads supporting (aligned to +/- strand):  ref base G (41/52);  new base C (2/4);  total (43/56)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.14e-01

AAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGCTAA  >  NC_000913/3231319‑3231613
                                                                                                                                                   |                                                                                                                                                   
aaGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTgg                                                                                                                                                    >  3:514277/1‑149 (MQ=255)
 aGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGa                                                                                                                                                   >  4:706081/1‑149 (MQ=255)
    gTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGc                                                                                                                                                >  1:193661/1‑149 (MQ=255)
    gTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCCCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGc                                                                                                                                                <  3:768351/149‑1 (MQ=255)
         aTCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTAt                                                                                                                                           <  3:762951/149‑1 (MQ=255)
          tCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATc                                                                                                                                          <  3:918019/149‑1 (MQ=255)
          gccggagaaaagggaaaaggggagtgcggggtttgtgatgggattcaaaaagaatacgcccggtgataacgccggatgcgacgcttgtcgcgttttatccaggttacggaattgcacatgtagGGCGGATGCGGCGTGGACGCCTTATc                                                                                                                                          >  4:443005‑M1/124‑149 (MQ=255)
              gCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGcc                                                                                                                                      >  3:1027480/1‑149 (MQ=255)
              gCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGcc                                                                                                                                      <  1:129473/149‑1 (MQ=255)
                 aCTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTa                                                                                                                                   <  4:653874/149‑1 (MQ=255)
                    gCACAGCGCGACTGGCGGATGCGGCGTGAACGCCTTATCAGCCCGACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACAt                                                                                                                                <  4:588760/149‑1 (MQ=255)
                      aCCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACAtgt                                                                                                                              <  4:440/149‑1 (MQ=255)
                          cgcCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCaa                                                                                                                          <  1:938677/149‑1 (MQ=255)
                          cgcCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCaa                                                                                                                          <  1:933734/149‑1 (MQ=255)
                          cgcCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCaa                                                                                                                          <  3:403055/149‑1 (MQ=255)
                           gcCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAAt                                                                                                                         >  1:994159/1‑149 (MQ=255)
                              aCTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTcc                                                                                                                      <  4:1071504/149‑1 (MQ=255)
                              aCTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTcc                                                                                                                      <  3:856391/149‑1 (MQ=255)
                                  gCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTaa                                                                                                                  <  4:657216/149‑1 (MQ=255)
                                    ggATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAAcc                                                                                                                >  4:670164/1‑149 (MQ=255)
                                     gATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCt                                                                                                               <  4:585697/149‑1 (MQ=255)
                                     gATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCt                                                                                                               <  1:953801/149‑1 (MQ=255)
                                     gATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCt                                                                                                               <  1:292821/149‑1 (MQ=255)
                                       tGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTgg                                                                                                             <  1:794980/149‑1 (MQ=255)
                                        gcggcgTGAACGCCTTATCCGCCCTTCATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTCGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGa                                                                                                            <  3:706081/149‑1 (MQ=255)
                                        gcggcgTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGa                                                                                                            <  3:251453/149‑1 (MQ=255)
                                             gttcaCGCCGCATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATaaaa                                                                                                       <  3:195723/145‑1 (MQ=255)
                                             gTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCAACTGGCGGATGCGGCGTGGACCCCTTATCCGCCCTACATTTGCAATTCCGTAACCTGGATaaaa                                                                                                       >  2:875812/1‑149 (MQ=255)
                                                 aCGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAAcgcg                                                                                                   >  2:589057/1‑149 (MQ=255)
                                                  cGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGa                                                                                                  >  4:508099/1‑149 (MQ=255)
                                                  cGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGACCTACATGTGCAATTCCGGAACCTGAATTAATCGCGa                                                                                                  <  2:502346/149‑1 (MQ=18)
                                                   gCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGAc                                                                                                 >  1:125545/1‑149 (MQ=255)
                                                     cTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACaa                                                                                               <  4:322069/149‑1 (MQ=255)
                                                        aTCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCg                                                                                            <  4:829583/149‑1 (MQ=255)
                                                             cccTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCa                                                                                       >  4:238233/1‑149 (MQ=255)
                                                                   aTGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCt                                                                                                               >  3:343748/1‑119 (MQ=255)
                                                                   aTGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCt                                                                                                               <  4:343748/119‑1 (MQ=255)
                                                                   aTGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGc                                                                                 <  2:482444/149‑1 (MQ=255)
                                                                   aTGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGc                                                                                 <  4:986723/149‑1 (MQ=255)
                                                                   aTGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGc                                                                                 >  1:380564/1‑149 (MQ=255)
                                                                     gtgtgtTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGt                                                                               <  2:151222/149‑1 (MQ=255)
                                                                     gtgtgtTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGt                                                                               <  4:970821/149‑1 (MQ=255)
                                                                     gtgtgtTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGt                                                                               >  1:150467/1‑149 (MQ=255)
                                                                        tgtTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGa                                                                                                  <  1:863279/127‑1 (MQ=255)
                                                                        tgtTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGa                                                                                                  >  2:863279/1‑127 (MQ=255)
                                                                        tgtTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTAt                                                                            >  3:33643/1‑149 (MQ=255)
                                                                          ttCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCa                                                                          >  3:233876/1‑149 (MQ=255)
                                                                             ccGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCg                                                                       >  4:772239/1‑149 (MQ=255)
                                                                             ccGTAGATCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCg                                                                       >  3:1002510/1‑149 (MQ=255)
                                                                              cGTAGGTCGGATAAGACGCGACAAGCGTCGCATCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCggg                                                                     >  2:42361/1‑149 (MQ=255)
                                                                                tAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCGCCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGc                                                                    <  2:771090/149‑1 (MQ=255)
                                                                                tAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGc                                                                    <  2:588693/149‑1 (MQ=255)
                                                                                tAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGc                                                                    <  1:630688/149‑1 (MQ=255)
                                                                                      ggATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGtt                                                                              <  4:267095/133‑1 (MQ=255)
                                                                                      ggATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGtt                                                                              >  3:267095/1‑133 (MQ=255)
                                                                                       gATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCt                                                             >  2:770735/1‑149 (MQ=255)
                                                                                         tAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCttt                                                           <  3:670164/149‑1 (MQ=255)
                                                                                          aaGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCtttt                                                          <  2:426875/149‑1 (MQ=255)
                                                                                           aGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCttttt                                                         <  4:947509/149‑1 (MQ=255)
                                                                                           aGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCttttt                                                         <  2:403799/149‑1 (MQ=255)
                                                                                           aGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCttttt                                                         <  1:770735/149‑1 (MQ=255)
                                                                                           aGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCttttt                                                         <  1:1058727/149‑1 (MQ=255)
                                                                                           aGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCttttt                                                         <  2:327575/149‑1 (MQ=255)
                                                                                             aCGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATc                                                                           <  3:481873/129‑1 (MQ=255)
                                                                                             aCGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATc                                                                           >  4:481873/1‑129 (MQ=255)
                                                                                                 gACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATcc                                                   >  2:930212/1‑149 (MQ=255)
                                                                                                    aaGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCAt                                                >  2:71856/1‑149 (MQ=255)
                                                                                                    aaGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCAt                                                >  3:908852/1‑149 (MQ=255)
                                                                                                     aGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATc                                               <  1:795380/149‑1 (MQ=255)
                                                                                                       cGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATcac                                             >  1:656782/1‑149 (MQ=255)
                                                                                                             aTCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAAccc                                       >  3:444501/1‑149 (MQ=255)
                                                                                                               ccGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCg                                     <  4:1002510/149‑1 (MQ=255)
                                                                                                                 ggCATCTGCACCGCGCCACTGGCGGATGCGGCGGGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCa                                   >  2:330235/1‑149 (MQ=255)
                                                                                                                  gCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTaa                                                                                                                  >  2:425919/1‑69 (MQ=255)
                                                                                                                  gCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTaa                                                                                                                  <  1:425919/69‑1 (MQ=255)
                                                                                                                  gCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCAc                                  >  3:608173/1‑149 (MQ=255)
                                                                                                                  gCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCAc                                  <  2:125545/149‑1 (MQ=255)
                                                                                                                    aTCTGCACCGTGTCACAGGCGGATGCGGCGTGGACGCCTTATCCGCTCTACTTGTGCAATTCCGTAAACTGGATAAAACGAGACAAGCGTCGCATCCGGCGTTATAACCGGGCGTATTCTTTTTGAATCCCACCACAAACCCCGCACTc                                <  2:419824/149‑1 (MQ=255)
                                                                                                                       tGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTcccc                             >  4:484412/1‑149 (MQ=255)
                                                                                                                        gCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCccct                            <  3:772239/149‑1 (MQ=255)
                                                                                                                          aCCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttt                          >  3:914402/1‑149 (MQ=255)
                                                                                                                            cgcgCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttttc                        <  2:656782/149‑1 (MQ=255)
                                                                                                                             gcgcCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttttcc                       <  1:490346/149‑1 (MQ=255)
                                                                                                                             gcgcCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttttcc                       <  1:149664/149‑1 (MQ=255)
                                                                                                                                ccACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATcc                                                   >  1:181716/1‑118 (MQ=255)
                                                                                                                                ccACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATcc                                                   <  2:181716/118‑1 (MQ=255)
                                                                                                                                ccACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCccctcttccctt                    <  2:965787/149‑1 (MQ=255)
                                                                                                                                 cACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttttcccttt                   >  4:959088/1‑149 (MQ=255)
                                                                                                                                 cACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttttcccttt                   >  3:27565/1‑149 (MQ=255)
                                                                                                                                  aCTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttttccctttt                  >  1:483649/1‑149 (MQ=255)
                                                                                                                                   agGGCGGATAAGGCGTCCACGCCGCATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGc                                                                                 >  4:1021329‑M1/3‑85 (MQ=255)
                                                                                                                                   agGGCGGATAAGGCGTCCACGCCGCATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGc                                                                                 <  3:1021329‑M1/83‑1 (MQ=255)
                                                                                                                                    gGGCGGATAAGGCGTCCACGCCGCATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCtttt                                                          <  4:202929‑M1/106‑1 (MQ=255)
                                                                                                                                    gGGCGGATAAGGCGTCCACGCCGCATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCtttt                                                          >  3:202929‑M1/2‑107 (MQ=255)
                                                                                                                                    gGGCGGATAAGGCGTCCACGCCGCATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTtct                <  2:789527‑M1/148‑1 (MQ=255)
                                                                                                                                     GGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTtctc               >  1:1073476‑M1/1‑149 (MQ=255)
                                                                                                                                      GCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTcc              >  3:238233‑M1/1‑149 (MQ=255)
                                                                                                                                         GATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGc           <  4:999025‑M1/149‑1 (MQ=255)
                                                                                                                                                 gTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGCTa   >  4:968180/1‑149 (MQ=255)
                                                                                                                                                  tGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGCTaa  >  2:548648/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGCTAA  >  NC_000913/3231319‑3231613

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