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% 274.2 / NA 84*310* (TAG→TAAyhcCputative Fe‑S oxidoreductase, Radical SAM superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (42/42);  total (42/42)

GGTCAACATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACCTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGAGCAAGCAGCGTCGGGCGAC  >  NC_000913/3352986‑3353254
                                                                                                                                       |                                                                                                                                     
ggTCAACATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTAc                                                                                                                                    >  2:1051032/1‑139 (MQ=255)
      cATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGtt                                                                                                                              <  1:132484/139‑1 (MQ=255)
          ataCTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                          <  1:458619/139‑1 (MQ=255)
          ataCTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                          <  2:1184704/139‑1 (MQ=255)
           taCTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTa                                                                                                                                     >  1:1350785/1‑127 (MQ=255)
           taCTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTa                                                                                                                                     <  2:1350785/127‑1 (MQ=255)
            aCTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGt                                                                                                                        >  1:1364713/1‑139 (MQ=255)
                cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAgg                                                                                                                         >  1:253795/1‑134 (MQ=255)
                cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAgg                                                                                                                         <  2:253795/134‑1 (MQ=255)
                cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGcc                                                                                                                    >  1:349730/1‑139 (MQ=255)
                cgcCAGCAGACGAAGTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGcc                                                                                                                    >  2:1113126/1‑139 (MQ=255)
                  ccAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCAc                                                                                                                  >  1:979844/1‑139 (MQ=255)
                          cGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTcc                                                                                                          <  2:586926/139‑1 (MQ=255)
                                      aTTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATcc                                                                                              >  1:1790506/1‑139 (MQ=255)
                                       ttAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATccc                                                                                             >  1:917062/1‑139 (MQ=255)
                                         aGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTg                                                                                           <  2:917062/139‑1 (MQ=255)
                                          gggAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGt                                                                                          <  2:1364713/139‑1 (MQ=255)
                                          gggAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGt                                                                                          <  2:349730/139‑1 (MQ=255)
                                                 ccTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTTAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCAt                                                                                   <  1:1113126/139‑1 (MQ=255)
                                                       cGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcat                                                                             <  2:1790506/139‑1 (MQ=255)
                                                       cGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcat                                                                             >  2:90064/1‑139 (MQ=255)
                                                        gACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcatt                                                                            <  2:979844/139‑1 (MQ=255)
                                                          cAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcattca                                                                          >  2:640125/1‑139 (MQ=255)
                                                           aaCCTGCCCCAATAATACCCAAGTCTGACAGTTACGTGATAAACTGTGTGAGAAATGCGGAATACCCCTCAACGACTTACTCCGTTGGAGGAAGCCACTGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAg                                                                         >  1:3222/1‑139 (MQ=255)
                                                              cTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATa                                                                      <  1:640125/139‑1 (MQ=255)
                                                                gCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGa                                                                                                               >  2:1087663/1‑96 (MQ=255)
                                                                gCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGa                                                                                                               <  1:1087663/96‑1 (MQ=255)
                                                                       taataCCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAg                                                             <  1:90064/139‑1 (MQ=255)
                                                                           aCCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCa                                                                                    >  2:429372/1‑112 (MQ=255)
                                                                           aCCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCa                                                                                    <  1:429372/112‑1 (MQ=255)
                                                                            cccAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGTGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGa                                                        >  2:1502289/1‑139 (MQ=255)
                                                                            cccAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATa                                                                      >  2:1582731/1‑125 (MQ=255)
                                                                            cccAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATa                                                                      <  1:1582731/125‑1 (MQ=255)
                                                                               aaTTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCa                                                     >  1:1689288/1‑139 (MQ=255)
                                                                                aTTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCg                                                                                                                                >  2:845468/1‑63 (MQ=255)
                                                                                aTTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCg                                                                                                                                <  1:845468/63‑1 (MQ=255)
                                                                                aTTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCa                                                                                                        <  1:408517/87‑1 (MQ=255)
                                                                                aTTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCa                                                                                                        >  2:408517/1‑87 (MQ=255)
                                                                                aTTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGt                                                                                          <  2:274426/101‑1 (MQ=255)
                                                                                aTTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGt                                                                                          >  1:274426/1‑101 (MQ=255)
                                                                                 ttCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATg                                                                                  >  1:1536575/1‑108 (MQ=255)
                                                                                 ttCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATg                                                                                  <  2:1536575/108‑1 (MQ=255)
                                                                                 ttCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCAtt                                                   >  1:872493/1‑139 (MQ=255)
                                                                                      aCAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGa                                                                                                               <  2:653309/74‑1 (MQ=255)
                                                                                      aCAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGa                                                                                                               >  1:653309/1‑74 (MQ=255)
                                                                                      aCAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGAt                                                                                                <  2:1018455/89‑1 (MQ=255)
                                                                                      aCAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGAt                                                                                                >  1:1018455/1‑89 (MQ=255)
                                                                                      aCAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAg                                                                     >  1:1797354/1‑116 (MQ=255)
                                                                                      aCAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAg                                                                     <  2:1797354/116‑1 (MQ=255)
                                                                                       cAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGAt                                                                                                >  1:1658904/1‑88 (MQ=255)
                                                                                       cAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGAt                                                                                                <  2:1658904/88‑1 (MQ=255)
                                                                                         gTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCa                                                                                                                   <  2:1433316/67‑1 (MQ=255)
                                                                                         gTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCa                                                                                                                   >  1:1433316/1‑67 (MQ=255)
                                                                                           tACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGcc                                                                                                   >  1:438401/1‑81 (MQ=255)
                                                                                           tACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGcc                                                                                                   <  2:438401/81‑1 (MQ=255)
                                                                                           tACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGc                                                                                      <  2:1660223/94‑1 (MQ=255)
                                                                                           tACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGc                                                                                      >  1:1660223/1‑94 (MQ=255)
                                                                                           tACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcattca                                                                          <  1:1814954/106‑1 (MQ=255)
                                                                                           tACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcattca                                                                          >  2:1814954/1‑106 (MQ=255)
                                                                                           tACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATc                                         >  1:1233133/1‑139 (MQ=255)
                                                                                              gTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGAt                                      <  2:1339379/139‑1 (MQ=255)
                                                                                                 ataAATTGTGTGAGAAATGTGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGAtttt                                   <  1:1502289/139‑1 (MQ=255)
                                                                                                  taAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                          <  1:1817589/51‑1 (MQ=255)
                                                                                                  taAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                          >  2:1817589/1‑51 (MQ=255)
                                                                                                  taAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATg                                                                                  >  1:1236583/1‑91 (MQ=255)
                                                                                                  taAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATg                                                                                  <  2:1236583/91‑1 (MQ=255)
                                                                                                    aaTTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAgcg                                                                                                     <  1:926635/70‑1 (MQ=255)
                                                                                                    aaTTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAgcg                                                                                                     >  2:926635/1‑70 (MQ=255)
                                                                                                    aaTTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATg                                                                                  <  2:811191/89‑1 (MQ=255)
                                                                                                    aaTTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATg                                                                                  >  1:811191/1‑89 (MQ=255)
                                                                                                      ttGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCCTCCATCGATTTTCGcac                              <  1:1051032/139‑1 (MQ=255)
                                                                                                        gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGAt                                                                                                >  1:859322/1‑71 (MQ=255)
                                                                                                        gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGAt                                                                                                <  2:859322/71‑1 (MQ=255)
                                                                                                        gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGcacca                            >  2:828386/1‑139 (MQ=255)
                                                                                                        gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGca                               >  1:773942/1‑136 (MQ=255)
                                                                                                        gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGca                               <  2:773942/136‑1 (MQ=255)
                                                                                                               aaaTGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGa                     <  2:872493/139‑1 (MQ=255)
                                                                                                                   gCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGagcaa                 >  1:1419120/1‑139 (MQ=255)
                                                                                                                     ggAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGcc                                                                                                   <  1:1590519/55‑1 (MQ=255)
                                                                                                                     ggAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGcc                                                                                                   >  2:1590519/1‑55 (MQ=255)
                                                                                                                              cTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCATCGCCGATCCCTGTACGCCATGtt                                                                                <  2:1325604/65‑1 (MQ=255)
                                                                                                                              cTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGAGCAAGCAGCGTCggg      >  2:790102/1‑139 (MQ=255)
                                                                                                                              cTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAACGCCGATCCCTGTACGCCATGtt                                                                                >  1:1325604/1‑65 (MQ=255)
                                                                                                                                  aCGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGAGCAAGCAGCGTCGGGCGAc  <  1:1730837/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
GGTCAACATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACCTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGAGCAAGCAGCGTCGGGCGAC  >  NC_000913/3352986‑3353254

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