Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,353,121 C→T *310* (TAG→TAA yhcC ← putative Fe‑S oxidoreductase, Radical SAM superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,353,1210CT100.0% 279.0 / NA 85*310* (TAG→TAAyhcCputative Fe‑S oxidoreductase, Radical SAM superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (43/42);  total (43/42)

GTCAACATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACCTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGAGCAAGCAGCGTCGGGCGACG  >  NC_000913/3352987‑3353255
                                                                                                                                      |                                                                                                                                      
gacAACATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAACTCCCCTCAACGACTTACt                                                                                                                                    <  1:853541/137‑1 (MQ=255)
    aCATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGt                                                                                                                                <  2:241294/139‑1 (MQ=255)
         ataCTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                           <  1:100930/139‑1 (MQ=255)
         ataCTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                           <  2:8760/139‑1 (MQ=255)
          taCTGCGCCAGCAGACGAATTTGCTTCATTAGGTAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAgg                                                                                                                          >  2:825926/1‑139 (MQ=255)
          taCTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACt                                                                                                                                    <  2:267524/129‑1 (MQ=255)
          taCTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACt                                                                                                                                    >  1:267524/1‑129 (MQ=255)
          taCTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAgg                                                                                                                          <  1:1280738/139‑1 (MQ=255)
               cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                           >  1:1047926/1‑133 (MQ=255)
               cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                           <  2:1047926/133‑1 (MQ=255)
               cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGcc                                                                                                                     >  1:1024648/1‑139 (MQ=255)
               cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGcc                                                                                                                     >  2:528031/1‑139 (MQ=255)
                gcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAACAATACCCAATTCTGACAGTTACGTTATAAGTTGTGTGAGAAATGCGGAATTCCCCTCAACTACTTACTCCGTTGGAGGtctgtc                                                                                                                    >  2:853541/1‑134 (MQ=255)
                  cagcagACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACg                                                                                                                  <  2:491065/139‑1 (MQ=255)
                    gcagACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGa                                                                                                                <  2:999480/139‑1 (MQ=255)
                    gcagACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGa                                                                                                                >  2:603999/1‑139 (MQ=255)
                                gCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGcc                                                                                                    <  2:271422/139‑1 (MQ=255)
                                gCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGcc                                                                                                    <  2:1272012/139‑1 (MQ=255)
                                gCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGcc                                                                                                    <  1:528031/139‑1 (MQ=255)
                                          ggAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTa                                                                                          <  1:603999/139‑1 (MQ=255)
                                             aTTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGc                                                                                       <  2:186788/139‑1 (MQ=255)
                                              ttCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                           >  2:121404/1‑102 (MQ=255)
                                              ttCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                           <  1:121404/102‑1 (MQ=255)
                                              ttCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGcc                                                                                      >  2:721225/1‑139 (MQ=255)
                                                      cGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcat                                                                              >  1:421516/1‑139 (MQ=255)
                                                      cGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcat                                                                              >  1:883810/1‑139 (MQ=255)
                                                       gACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcatt                                                                             <  1:121409/139‑1 (MQ=255)
                                                        aCAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGtt                                                                                 >  2:469634/1‑134 (MQ=255)
                                                        aCAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGtt                                                                                 <  1:469634/134‑1 (MQ=255)
                                                         cAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcattca                                                                           >  2:951818/1‑139 (MQ=255)
                                                           aCCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGa                                                                         >  2:54014/1‑139 (MQ=255)
                                                            ccTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGc                                                                                       <  2:1274244/124‑1 (MQ=255)
                                                            ccTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGc                                                                                       >  1:1274244/1‑124 (MQ=255)
                                                            ccTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcat                                                                              <  2:127065/133‑1 (MQ=255)
                                                            ccTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcat                                                                              >  1:127065/1‑133 (MQ=255)
                                                            ccTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGAt                                                                        >  2:253237/1‑139 (MQ=255)
                                                             cTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATa                                                                       <  1:721225/139‑1 (MQ=255)
                                                                     ataataCCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATg                                                                                   >  1:917813/1‑119 (MQ=255)
                                                                     ataataCCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTAGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATg                                                                                   <  2:917813/119‑1 (MQ=255)
                                                                      taataCCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                           >  1:84078/1‑78 (MQ=255)
                                                                      taataCCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                           <  2:84078/78‑1 (MQ=255)
                                                                          aCCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCg                                                          >  1:897899/1‑139 (MQ=255)
                                                                           cccAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGAt                                                                  >  2:491506/1‑130 (MQ=255)
                                                                           cccAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGAt                                                                  <  1:491506/130‑1 (MQ=255)
                                                                           cccAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGa                                                         >  2:874043/1‑139 (MQ=255)
                                                                           cccAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGa                                                         <  1:874043/139‑1 (MQ=255)
                                                                               aTTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAgcg                                                                                                      >  2:1151262/1‑90 (MQ=255)
                                                                               aTTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAgcg                                                                                                      <  1:1151262/90‑1 (MQ=255)
                                                                                ttCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGa                                                         >  2:441847/1‑134 (MQ=255)
                                                                                ttCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGa                                                         <  1:441847/134‑1 (MQ=255)
                                                                                     aCAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATACCGATCCAGCTCGACCATTCCCGt                                               >  2:999175/1‑139 (MQ=255)
                                                                                         ttACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGAt                                                                                                 <  1:1216624/85‑1 (MQ=255)
                                                                                         ttACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGAt                                                                                                 >  2:1216624/1‑85 (MQ=255)
                                                                                          tACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTACGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATc                                          >  2:690366/1‑139 (MQ=255)
                                                                                                      tgtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGa                                                                                                                <  1:197274/57‑1 (MQ=255)
                                                                                                      tgtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGa                                                                                                                >  2:197274/1‑57 (MQ=255)
                                                                                                       gtgtgAGAAATGCGGAATTCCCCTCAACGACTTAc                                                                                                                                     <  2:430242/35‑1 (MQ=255)
                                                                                                       gtgtgAGAAATGCGGAATTCCCCTCAACGACTTAc                                                                                                                                     >  1:430242/1‑35 (MQ=255)
                                                                                                       gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGAc                                                                                                               >  2:1112815/1‑57 (MQ=255)
                                                                                                       gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGAc                                                                                                               <  1:1112815/57‑1 (MQ=255)
                                                                                                        tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGcc                                                                                                    >  2:473012/1‑67 (MQ=255)
                                                                                                        tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGcc                                                                                                    <  1:473012/67‑1 (MQ=255)
                                                                                                        tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGAt                                       >  1:162739/1‑128 (MQ=255)
                                                                                                        tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGAt                                       <  2:162739/128‑1 (MQ=255)
                                                                                                        tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGcaccac                            >  1:285878/1‑139 (MQ=255)
                                                                                                        tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGcaccac                            >  1:940253/1‑139 (MQ=255)
                                                                                                         gtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACcaca                           <  2:883810/139‑1 (MQ=255)
                                                                                                          tgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAg                          <  2:897899/139‑1 (MQ=255)
                                                                                                           gagaAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATc                                                                                                <  2:939590/68‑1 (MQ=255)
                                                                                                           gagaAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATc                                                                                                >  1:939590/1‑68 (MQ=255)
                                                                                                            agaAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGt                                                                                                                         >  2:595368/1‑42 (MQ=255)
                                                                                                            agaAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGt                                                                                                                         <  1:595368/42‑1 (MQ=255)
                                                                                                              aaaTGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGa                      <  2:285878/139‑1 (MQ=255)
                                                                                                               aaTGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCg                                                                                                   >  1:342015/1‑61 (MQ=255)
                                                                                                               aaTGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCg                                                                                                   <  2:342015/61‑1 (MQ=255)
                                                                                                                aTGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGAt                                                                                                 >  2:73342/1‑62 (MQ=255)
                                                                                                                aTGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGAt                                                                                                 <  1:73342/62‑1 (MQ=255)
                                                                                                                  gCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGagcaa                  >  2:43599/1‑139 (MQ=255)
                                                                                                                  gCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCAACACAGCGGagcaa                  >  1:944982/1‑139 (MQ=255)
                                                                                                                   cGGAATTCCCCTCAACGACTTACTCCGTTGGAGGGAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGagcaag                 <  1:1115831/139‑1 (MQ=255)
                                                                                                                           cccTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGAt                                                                                                 >  2:923155/1‑51 (MQ=255)
                                                                                                                           cccTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGAt                                                                                                 <  1:923155/51‑1 (MQ=255)
                                                                                                                                  cGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGTGCCGATCCCTGTACGCCATg                                                                                   <  2:1098388/58‑1 (MQ=255)
                                                                                                                                  cGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGTGCCGATCCCTGTACGCCATg                                                                                   >  1:1098388/1‑58 (MQ=255)
                                                                                                                                  cGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGAGCAAGCAGCGTCGGGCGACg  >  2:94013/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                      
GTCAACATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACCTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGAGCAAGCAGCGTCGGGCGACG  >  NC_000913/3352987‑3353255

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