Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,200,881 G→A A51V (GCC→GTC)  xisE ← e14 prophage, putative excisionase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,200,8810GA100.0% 355.0 / NA 107A51V (GCC→GTC) xisEe14 prophage, putative excisionase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (55/52);  total (55/52)

Only 100 of 107 total aligned reads displayed.

AGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCA  >  NC_000913/1200754‑1201012
                                                                                                                               |                                                                                                                                   
aGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGc                                                                                                                          >  2:1508300/1‑139 (MQ=255)
aGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGc                                                                                                                          >  2:772822/1‑139 (MQ=255)
  cGGATATAAGTTAGGGACTGCAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTTCTCACGctc                                                                                                                        >  2:1019373/1‑137 (MQ=255)
  cGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCAc                                                                                                                        >  2:664016/1‑139 (MQ=255)
    gATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTc                                                                                                                      >  2:1176996/1‑139 (MQ=255)
     atatAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCg                                                                                                                     >  2:278299/1‑139 (MQ=255)
               gggACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCgg                                                                                                           >  2:1349198/1‑139 (MQ=255)
                  aCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCt                                                                                                        <  1:772822/139‑1 (MQ=255)
                    tGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTg                                                                                                      >  1:572955/1‑139 (MQ=255)
                     gAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGt                                                                                                     >  1:346340/1‑139 (MQ=255)
                      aaaCATTATTCTTACGCGGTATGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGTCATAACCGGCCTTGtt                                                                                                    >  1:1176239/1‑139 (MQ=255)
                             attCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGc                                                                                             <  1:1508300/139‑1 (MQ=255)
                                 ttACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGt                                                                                         <  1:977380/139‑1 (MQ=255)
                                        gTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                     <  2:363390/136‑1 (MQ=255)
                                        gTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                     >  1:363390/1‑136 (MQ=255)
                                         tCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCg                                                                                           <  1:500991/129‑1 (MQ=255)
                                         tCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACaaaa                                                                                 >  2:910608/1‑139 (MQ=255)
                                           tGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACaaataa                                                                               >  2:46120/1‑136 (MQ=255)
                                             gCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAaaataaa                                                                             <  1:664016/139‑1 (MQ=255)
                                               tGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                           <  2:572955/139‑1 (MQ=255)
                                               tGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                           <  1:210872/139‑1 (MQ=255)
                                                  cATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTg                                                                        <  1:1730064/139‑1 (MQ=255)
                                                  cATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTg                                                                        <  2:346340/139‑1 (MQ=255)
                                                   aTTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGc                                                                                                                    >  2:422669/1‑94 (MQ=255)
                                                   aTTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGc                                                                                                                    <  1:422669/94‑1 (MQ=255)
                                                   aTTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGcca                                                                                                        <  2:1193046/106‑2 (MQ=255)
                                                   aTTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGcca                                                                                                        >  1:1193046/1‑105 (MQ=255)
                                                   aTTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGc                                                                       >  2:266780/1‑139 (MQ=255)
                                                   aTTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGc                                                                       >  1:1508477/1‑139 (MQ=255)
                                                                  ttGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTg                                                                        <  1:603736/123‑1 (MQ=255)
                                                                  ttGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTg                                                                        >  2:603736/1‑123 (MQ=255)
                                                                  ttGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGc                                                        >  1:638576/1‑139 (MQ=255)
                                                                  ttGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGc                                                        >  1:744256/1‑139 (MQ=255)
                                                                  ttGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGc                                                        >  1:966361/1‑139 (MQ=255)
                                                                  ttGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGc                                                        <  2:638576/139‑1 (MQ=255)
                                                                    gCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCaa                                                      >  2:629233/1‑139 (MQ=255)
                                                                     cAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAAt                                                     <  1:1019373/139‑1 (MQ=255)
                                                                         aGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAATAAATTCTGCTTAGCATAACGACGCAATGTGGa                                                <  1:46120/139‑1 (MQ=255)
                                                                          gTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGa                                                <  1:629233/139‑1 (MQ=255)
                                                                            aaTGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGaca                                              <  1:1349198/139‑1 (MQ=255)
                                                                            aaTGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGaca                                              <  2:966361/139‑1 (MQ=255)
                                                                            aaTGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGaca                                              <  2:744256/139‑1 (MQ=255)
                                                                            aaTGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGaca                                              <  2:1508477/139‑1 (MQ=255)
                                                                             aTGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACGCGACGTCCTCACGCACTCGCCATAACCGACCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGacac                                             >  1:1162411/1‑139 (MQ=255)
                                                                                aGTCATTTTTCTTGATTACAGGTGTTACCAACTCCACAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACt                                          <  1:1046831/139‑1 (MQ=255)
                                                                                     tttttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATaa                                                                                                               >  1:551346/1‑65 (MQ=255)
                                                                                     tttttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATaa                                                                                                               <  2:551346/65‑1 (MQ=255)
                                                                                     tttttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTgg                                       >  2:1603998/1‑137 (MQ=255)
                                                                                     tttttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTgg                                       <  1:1603998/137‑1 (MQ=255)
                                                                                     tttttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAAAGTGGa                                                >  2:1597761/1‑128 (MQ=255)
                                                                                     tttttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAAAGTGGa                                                <  1:1597761/128‑1 (MQ=255)
                                                                                        ttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGATGCAATGTGGACACACTTGGAGGAt                                  <  2:1176239/139‑1 (MQ=255)
                                                                                                tACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGc                                                                                                                    <  2:1729542/49‑1 (MQ=255)
                                                                                                tACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGc                                                                                                                    >  1:1729542/1‑49 (MQ=255)
                                                                                                       gtTACCAACTCCCCAACTAACTCGACGTCCTCACg                                                                                                                           >  2:1719258/1‑35 (MQ=255)
                                                                                                       gtTACCAACTCCCCAACTAACTCGACGTCCTCACg                                                                                                                           <  1:1719258/35‑1 (MQ=255)
                                                                                                        ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAAc                                                                                                              >  2:1689578/1‑47 (MQ=255)
                                                                                                        ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAAc                                                                                                              <  1:1689578/47‑1 (MQ=255)
                                                                                                        ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCg                                                                                           >  1:1252295/1‑66 (MQ=255)
                                                                                                        ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCg                                                                                           <  2:1252295/66‑1 (MQ=255)
                                                                                                        ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGa                                                                                      <  2:973177/71‑1 (MQ=255)
                                                                                                        ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGa                                                                                      >  2:329380/1‑71 (MQ=255)
                                                                                                        ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGa                                                                                      >  1:973177/1‑71 (MQ=255)
                                                                                                        ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACaaa                                                                                  >  1:1330282/1‑75 (MQ=255)
                                                                                                        ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACaaa                                                                                  <  2:1330282/75‑1 (MQ=255)
                                                                                                        ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                           >  2:637527/1‑82 (MQ=255)
                                                                                                          aCCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATaa                                                                                                               <  2:41975/44‑1 (MQ=255)
                                                                                                          aCCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATaa                                                                                                               >  1:41975/1‑44 (MQ=255)
                                                                                                          aCCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                     >  2:487999/1‑70 (MQ=255)
                                                                                                          aCCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                     <  1:487999/70‑1 (MQ=255)
                                                                                                          aCCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACaaa                                                                                  <  2:364250/73‑1 (MQ=255)
                                                                                                          aCCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACaaa                                                                                  >  1:364250/1‑73 (MQ=255)
                                                                                                          aCCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgacg                                                         >  2:1561025/1‑98 (MQ=255)
                                                                                                          aCCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTcag                <  1:1176996/139‑1 (MQ=255)
                                                                                                          aCCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTcag                <  1:910608/139‑1 (MQ=255)
                                                                                                              aCTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTg                                                                        <  2:214713/79‑1 (MQ=255)
                                                                                                              aCTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTg                                                                        >  1:214713/1‑79 (MQ=255)
                                                                                                               cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGcc                                                                                                         >  1:581427/1‑45 (MQ=255)
                                                                                                               cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGcc                                                                                                         <  2:581427/45‑1 (MQ=255)
                                                                                                               cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCg                                                                                           >  2:1592611/1‑59 (MQ=255)
                                                                                                               cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCg                                                                                           <  1:1592611/59‑1 (MQ=255)
                                                                                                               cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacga                                                           <  2:301733/91‑1 (MQ=255)
                                                                                                               cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacga                                                           >  1:301733/1‑91 (MQ=255)
                                                                                                               cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGcc           >  1:1259087/1‑139 (MQ=255)
                                                                                                                tCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGccc          >  1:1697268/1‑139 (MQ=255)
                                                                                                                 ccccAACTAACTCGACGTACTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCa         <  2:1162411/139‑1 (MQ=255)
                                                                                                                  cccAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAt        >  2:378375/1‑139 (MQ=255)
                                                                                                                  cccAACTAACTCGACGTCCTCACGCACTCGCCATAACCGACCTTGTTTCATGGCCGGTGGACaaa                                                                                  >  1:1257169/1‑65 (MQ=255)
                                                                                                                     aactaactCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAttct     >  1:156665/1‑139 (MQ=255)
                                                                                                                        taactCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCa  <  1:278299/139‑1 (MQ=255)
                                                                                                                          actCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCg                                                                                           >  1:1385500/1‑48 (MQ=255)
                                                                                                                          actCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCg                                                                                           <  2:1385500/48‑1 (MQ=255)
                                                                                                                          actCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGc                                                                       >  1:1389266/1‑68 (MQ=255)
                                                                                                                          actCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGc                                                                       <  2:1389266/68‑1 (MQ=255)
                                                                                                                          actCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCt                                                                      <  1:1630508/69‑1 (MQ=255)
                                                                                                                          actCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCt                                                                      >  2:1630508/1‑69 (MQ=255)
                                                                                                                          actCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGa                                                >  2:66479/1‑91 (MQ=255)
                                                                                                                          actCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGa                                                <  1:66479/91‑1 (MQ=255)
                                                                                                                           ctCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgacg                                                         <  2:488962/81‑1 (MQ=255)
                                                                                                                           ctCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgacg                                                         >  1:488962/1‑81 (MQ=255)
                                                                                                                               |                                                                                                                                   
AGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCA  >  NC_000913/1200754‑1201012

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