Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,082,381 C→A 6.2% D127Y (GAC→TAC)  fdoH ← formate dehydrogenase‑O, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,082,3810CA 6.2% 290.8 / 4.8 81D127Y (GAC→TAC) fdoHformate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base C (36/40);  new base A (3/2);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

TTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAA  >  NC_000913/4082244‑4082516
                                                                                                                                         |                                                                                                                                       
ttCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTcga                                                                                                                                       <  2:155941/140‑1 (MQ=255)
ttCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTcga                                                                                                                                       <  1:1128081/140‑1 (MQ=255)
    tGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGAt                                                                                                                                   <  1:1061258/140‑1 (MQ=255)
     gCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATa                                                                                                                                  >  1:845484/1‑140 (MQ=255)
     gCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATa                                                                                                                                  >  2:392718/1‑140 (MQ=255)
     gCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTAGACGATa                                                                                                                                  >  1:701052/1‑140 (MQ=255)
     gCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTAGACGATa                                                                                                                                  >  1:701089/1‑140 (MQ=255)
     gCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGGTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCACTAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATa                                                                                                                                  >  1:1106609/1‑140 (MQ=255)
      cccAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATAc                                                                                                                                 <  1:747845/140‑1 (MQ=255)
      cccAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATAc                                                                                                                                 <  1:935554/140‑1 (MQ=255)
         aaCCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGt                                                                                                                              >  2:251733/1‑140 (MQ=255)
          accaccACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGtt                                                                                                                             >  1:899357/1‑140 (MQ=255)
           ccaccaCGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTg                                                                                                                            <  2:46874/140‑1 (MQ=255)
           ccaccaCGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTg                                                                                                                            <  2:46872/140‑1 (MQ=255)
           ccaccaCGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTg                                                                                                                            <  2:44890/140‑1 (MQ=255)
              ccacGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCa                                                                                                                         <  2:515457/140‑1 (MQ=255)
              ccacGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCa                                                                                                                         <  1:713956/140‑1 (MQ=255)
                   cgGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTg                                                                                                                    <  1:715303/140‑1 (MQ=255)
                   cgGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTg                                                                                                                    <  1:715287/140‑1 (MQ=255)
                      tCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGaat                                                                                                                 >  2:402079/1‑140 (MQ=255)
                        aaCGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGaataa                                                                                                               >  2:656753/1‑140 (MQ=255)
                        aaCGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGaataa                                                                                                               >  2:656758/1‑140 (MQ=255)
                        aaCGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGaataa                                                                                                               >  2:656686/1‑140 (MQ=255)
                        aaCGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGaataa                                                                                                               >  2:531083/1‑140 (MQ=255)
                         aCGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGaataat                                                                                                              >  2:1032850/1‑140 (MQ=255)
                                gCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCt                                                                                                       >  1:518537/1‑140 (MQ=255)
                                gCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCt                                                                                                       >  2:732353/1‑140 (MQ=255)
                                gCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCt                                                                                                       >  2:556691/1‑140 (MQ=255)
                                    aCATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTccg                                                                                                   >  1:478314/1‑140 (MQ=255)
                                           tAGACGCGGTTGTCTTCCGGTTTTGGGCGCGGAATGTCTAACGTACAGCCCGCAATGCTATAACCGCAGCCAATGCACTGCTCGGCATGGAAGACGTCGATACCGTTGGCATCCTGAATCATTGCACCTTCCGCCGGGCa                                                                                            >  2:333054/1‑140 (MQ=255)
                                           tAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCa                                                                                            >  2:988849/1‑140 (MQ=255)
                                           tAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCa                                                                                            >  2:988791/1‑140 (MQ=255)
                                           tAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCa                                                                                            <  2:547196/140‑1 (MQ=255)
                                           tAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCa                                                                                            >  1:982753/1‑140 (MQ=255)
                                           tAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCa                                                                                            >  1:417600/1‑140 (MQ=255)
                                             gACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGACGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACg                                                                                          >  2:400082/1‑140 (MQ=255)
                                                 cgGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCttt                                                                                      >  1:400609/1‑140 (MQ=255)
                                                  ggTTATCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTAGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTc                                                                                     >  2:206784/1‑140 (MQ=255)
                                                     tGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAgg                                                                                  <  1:471912/140‑1 (MQ=255)
                                                     tGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAgg                                                                                  <  1:795173/140‑1 (MQ=255)
                                                     tGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAgg                                                                                  <  1:795120/140‑1 (MQ=255)
                                                     tGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAgg                                                                                  <  1:373839/140‑1 (MQ=255)
                                                     tGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAgg                                                                                  <  2:169621/140‑1 (MQ=255)
                                                          tCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGcc                                                                             >  1:1190658/1‑140 (MQ=255)
                                                          tCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGcc                                                                             >  1:1190738/1‑140 (MQ=255)
                                                          tCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGcc                                                                             <  1:1231313/140‑1 (MQ=255)
                                                          tCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGcc                                                                             <  1:380158/140‑1 (MQ=255)
                                                              ggTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGa                                                                         >  2:896804/1‑140 (MQ=255)
                                                               gTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGAt                                                                        <  2:526181/140‑1 (MQ=255)
                                                               gTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGAt                                                                        <  1:1118588/140‑1 (MQ=255)
                                                                 tGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGAGGATACCGTTGGCATAGTGAATAATTGCACCTTCCGCCGGGCACGCGTTCAGGCAGCCGGGATCg                                                                      >  2:1135535/1‑140 (MQ=255)
                                                                    ggCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATTCTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGaa                                                                   <  1:162393/140‑1 (MQ=255)
                                                                    ggCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGaa                                                                   <  1:162307/140‑1 (MQ=255)
                                                                      cgcgGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACa                                                                 <  1:1226413/140‑1 (MQ=255)
                                                                        cgGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGt                                                               >  1:578931/1‑140 (MQ=255)
                                                                        cgGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGt                                                               <  1:634993/140‑1 (MQ=255)
                                                                        cgGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGt                                                               >  1:90925/1‑140 (MQ=255)
                                                                        cgGAATGTCGAACGGACAGCCCGCAAAGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGACGACGATACCGATGGCATACTGAATAAATGCACCTTCCGCCCGGCACGCTTTCAGGCAGCCTGGATCGTAACAGt                                                               <  2:1210419/140‑1 (MQ=255)
                                                                          gAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGc                                                             <  1:714152/140‑1 (MQ=255)
                                                                               tCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTAGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACa                                                        <  2:701089/140‑1 (MQ=255)
                                                                               tCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTAGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACa                                                        <  2:701052/140‑1 (MQ=255)
                                                                                         aGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTg                                              >  2:449284/1‑140 (MQ=255)
                                                                                               caatgcaatAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATc                                        <  2:899357/140‑1 (MQ=255)
                                                                                                 atgcaatAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAg                                      <  2:960685/140‑1 (MQ=255)
                                                                                                 atgcaatAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAg                                      <  1:1127558/140‑1 (MQ=255)
                                                                                                     aatAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCAt                                  >  1:442463/1‑140 (MQ=255)
                                                                                                     aatAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCAt                                  >  1:442894/1‑140 (MQ=255)
                                                                                                       tAACCGCAGCCAATGCACTGCTCGGACTGGAAGGCGACGATACCGTTGGCCTACTGAATAATTGCACCGTCCGCCGGGCCCGCTTTCAGGCAGCCGGGATCGGAACAGGGCAGACAGCCGGCTTGGCGGAGCAGcccgtc                                >  2:596871/1‑136 (MQ=255)
                                                                                                         aCCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCa                              <  2:835720/140‑1 (MQ=255)
                                                                                                                         tGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTc              <  1:190351/140‑1 (MQ=255)
                                                                                                                         tGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTc              <  1:977137/140‑1 (MQ=255)
                                                                                                                          gCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTcc             >  2:559740/1‑140 (MQ=255)
                                                                                                                           cTCGGACTGGAAGTCGACGATACCGTTGGCATACTGACTAATAGCACCTTCCGCCGGGCACGCTTTCAGGCAGGCTGGATCGGAACACTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCa            <  2:500197/140‑1 (MQ=255)
                                                                                                                           cTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCa            <  1:1013280/140‑1 (MQ=255)
                                                                                                                           cTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCa            >  1:388404/1‑140 (MQ=255)
                                                                                                                              ggACTGGAAGTCGACGATACCGTTGGCATACGGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGTACTGTGCATACAGCCGTCTGTGCGGATCTGCCAGTCGAGTGTGTCGGTCTGCTCCACtt         >  2:1222169/1‑140 (MQ=255)
                                                                                                                                   ggAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCgag    <  2:1179691/140‑1 (MQ=255)
                                                                                                                                   ggAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCgag    <  2:893113/140‑1 (MQ=255)
                                                                                                                                     aaGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATGCAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGaa  <  2:1190658/140‑1 (MQ=255)
                                                                                                                                     aaGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGaa  <  2:1190738/140‑1 (MQ=255)
                                                                                                                                     aaGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGaa  >  1:120432/1‑140 (MQ=255)
                                                                                                                                     aaGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGaa  <  2:975176/140‑1 (MQ=255)
                                                                                                                                     aaGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTACATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGaa  <  2:61817/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
TTCTTGCCCAACCACCACGCGGTCAACGCACAGCGTACATTTGTAGACGCGGTTGTCTTCCGGGTTGAGGCGCGGAATGTCGAACGGACAGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAA  >  NC_000913/4082244‑4082516

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