Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,465,938 A→T V132V (GTA→GTT)  gspO → bifunctional prepilin leader peptidase/ methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,465,9380AT100.0% 274.5 / NA 84V132V (GTA→GTT) gspObifunctional prepilin leader peptidase/ methylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (42/42);  total (42/42)

GCTCAGTGTTATTGATTTTCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTAGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACGTTGCCAATG  >  NC_000913/3465809‑3466070
                                                                                                                                 |                                                                                                                                    
gCTCAGTGTTATTGATTTTCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTAt                                                                                                                             >  1:1367405/1‑139 (MQ=255)
 cTCAGTGTTATTGATTTTCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATg                                                                                                                            <  2:734802/139‑1 (MQ=255)
      tgtTATTGATTTTCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTg                                                                                                                       >  1:666794/1‑139 (MQ=255)
           ttGATTTTCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCtgtg                                                                                                                  <  2:1505273/139‑1 (MQ=255)
                tttCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGgtgt                                                                                                              <  1:701665/138‑1 (MQ=255)
                tttCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGgtgt                                                                                                              >  2:701665/1‑138 (MQ=255)
                  tCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATggg                                                                                                                          >  2:1375069/1‑124 (MQ=255)
                  tCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATggg                                                                                                                          <  1:1375069/124‑1 (MQ=255)
                  tCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGtt                                                                                                           >  2:761889/1‑139 (MQ=255)
                        tCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTgg                                                                                                     <  2:1214060/139‑1 (MQ=255)
                        tCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTgg                                                                                                     <  1:494300/139‑1 (MQ=255)
                          aGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTgggg                                                                                                   >  2:150029/1‑139 (MQ=255)
                               cTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCt                                                                                              >  2:209438/1‑139 (MQ=255)
                                  cccGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGt                                                                                           <  2:236674/139‑1 (MQ=255)
                                  cccGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGt                                                                                           <  2:426472/139‑1 (MQ=255)
                                          aCTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGt                                                                                   <  1:1183752/139‑1 (MQ=255)
                                                 ttACCGCTTCTCTGGCTTGGCTTGGGATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTGACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAag                                                                            <  1:62874/139‑1 (MQ=255)
                                                 ttACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAag                                                                            <  1:209438/139‑1 (MQ=255)
                                                        tgctCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGAtt                                                                     >  1:1541584/1‑139 (MQ=255)
                                                               gcttggcttgGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGt                                                                                       >  1:1239500/1‑114 (MQ=255)
                                                               gcttggcttgGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGt                                                                                       <  2:1239500/114‑1 (MQ=255)
                                                                    gcttggTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCg                                                                                                                                   <  1:1232426/65‑1 (MQ=255)
                                                                    gcttggTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCg                                                                                                                                   >  2:1232426/1‑65 (MQ=255)
                                                                         gTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGttt                                                                                                          >  1:1248586/1‑85 (MQ=255)
                                                                         gTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGttt                                                                                                          <  2:1248586/85‑1 (MQ=255)
                                                                                   cgcAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGc                                          >  1:162182/1‑139 (MQ=255)
                                                                                   cgcAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGc                                          <  2:1541584/139‑1 (MQ=255)
                                                                                    gcAGTATGGATTAATGGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGcc                                         >  2:429472/1‑139 (MQ=255)
                                                                                             aTTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTgtggt                                                                                                                >  2:695071/1‑59 (MQ=255)
                                                                                             aTTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTgtggt                                                                                                                <  1:695071/59‑1 (MQ=255)
                                                                                               tAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCt                                                                                                                     <  1:1408545/52‑1 (MQ=255)
                                                                                               tAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCt                                                                                                                     >  2:1408545/1‑52 (MQ=255)
                                                                                                  ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATggg                                                                                                                          <  2:994349/44‑1 (MQ=255)
                                                                                                  ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATggg                                                                                                                          >  1:994349/1‑44 (MQ=255)
                                                                                                  ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCtgt                                                                                                                   >  2:993338/1‑51 (MQ=255)
                                                                                                  ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCtgt                                                                                                                   <  1:993338/51‑1 (MQ=255)
                                                                                                  ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGt                                                                                                <  1:228529/70‑1 (MQ=255)
                                                                                                  ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGt                                                                                                >  2:228529/1‑70 (MQ=255)
                                                                                                  ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGttt                                                                                     >  2:739270/1‑81 (MQ=255)
                                                                                                  ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGttt                                                                                     <  1:739270/81‑1 (MQ=255)
                                                                                                  ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACaaa                                                                             <  2:573126/89‑1 (MQ=255)
                                                                                                  ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACaaa                                                                             >  1:573126/1‑89 (MQ=255)
                                                                                                  ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGc                                          <  2:529980/124‑1 (MQ=255)
                                                                                                  ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGc                                          >  1:529980/1‑124 (MQ=255)
                                                                                                  ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCa                                   <  2:369732/131‑1 (MQ=255)
                                                                                                  ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCa                                   >  1:369732/1‑131 (MQ=255)
                                                                                                  ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATg                           >  1:769845/1‑139 (MQ=255)
                                                                                                  ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATg                           >  1:1203554/1‑139 (MQ=255)
                                                                                                  ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATg                           >  2:358511/1‑139 (MQ=255)
                                                                                                  ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATg                           >  2:340151/1‑139 (MQ=255)
                                                                                                      tttACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGgtgt                       >  2:1386003/1‑139 (MQ=255)
                                                                                                       ttACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGa                                                                       <  1:1252785/90‑1 (MQ=255)
                                                                                                       ttACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGa                                                                       >  2:1252785/1‑90 (MQ=255)
                                                                                                               tgCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTgtggt                                                                                                                >  2:956122/1‑41 (MQ=255)
                                                                                                               tgCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTgtggt                                                                                                                <  1:956122/41‑1 (MQ=255)
                                                                                                               tgCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGt                                                                                           >  1:652227/1‑62 (MQ=255)
                                                                                                               tgCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGt                                                                                           <  2:652227/62‑1 (MQ=255)
                                                                                                               tgCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCaa              >  2:905697/1‑139 (MQ=255)
                                                                                                                gCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCt                                                                <  1:295231/88‑1 (MQ=255)
                                                                                                                gCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCt                                                                >  2:295231/1‑88 (MQ=255)
                                                                                                                  ggTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACg           <  2:1203554/139‑1 (MQ=255)
                                                                                                                    tttACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTca                                                                                 <  2:817829/67‑1 (MQ=255)
                                                                                                                    tttACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTca                                                                                 >  1:817829/1‑67 (MQ=255)
                                                                                                                    tttACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACGtt         <  1:150029/139‑1 (MQ=255)
                                                                                                                       aCGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGc                                          <  2:1328226/103‑1 (MQ=255)
                                                                                                                       aCGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGc                                          >  1:1328226/1‑103 (MQ=255)
                                                                                                                        cGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACGTTGCCa     >  1:666036/1‑139 (MQ=255)
                                                                                                                         ggCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACt                                                                                                       <  2:350775/40‑1 (MQ=255)
                                                                                                                         ggCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACt                                                                                                       >  1:350775/1‑40 (MQ=255)
                                                                                                                         ggCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGc                                          >  2:411641/1‑101 (MQ=255)
                                                                                                                         ggCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGc                                          <  1:411641/101‑1 (MQ=255)
                                                                                                                         ggCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACGTTGCCaa    <  1:1269592/139‑1 (MQ=255)
                                                                                                                         ggCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACGTTGCCaa    >  2:1269592/1‑139 (MQ=255)
                                                                                                                         ggCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACGTTGCCaa    <  1:761889/139‑1 (MQ=255)
                                                                                                                           cgcgGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACGTTGCCAATg  <  2:1466412/139‑1 (MQ=255)
                                                                                                                           cgcgGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACGTTGCCAATg  <  2:162182/139‑1 (MQ=255)
                                                                                                                             cgGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGt                                                                                                >  2:486859/1‑43 (MQ=255)
                                                                                                                             cgGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGt                                                                                                <  1:486859/43‑1 (MQ=255)
                                                                                                                             cgGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTca                                                                                 >  1:1506367/1‑58 (MQ=255)
                                                                                                                             cgGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTca                                                                                 <  2:1506367/58‑1 (MQ=255)
                                                                                                                             cgGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACGtt         >  2:476875/1‑130 (MQ=255)
                                                                                                                             cgGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACGtt         <  1:476875/130‑1 (MQ=255)
                                                                                                                              ggTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGa                                                                       <  2:592086/67‑1 (MQ=255)
                                                                                                                              ggTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGa                                                                       >  1:592086/1‑67 (MQ=255)
                                                                                                                                 |                                                                                                                                    
GCTCAGTGTTATTGATTTTCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTAGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACGTTGCCAATG  >  NC_000913/3465809‑3466070

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: ‑