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% 275.1 / NA 85A51V (GCC→GTC) xisEe14 prophage, putative excisionase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (42/43);  total (42/43)

CTACTATAAAGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCA  >  NC_000913/1200745‑1201005
                                                                                                                                        |                                                                                                                            
ctactaTAAAGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACg                                                                                                                            <  2:552999/139‑1 (MQ=255)
         aGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGc                                                                                                                   >  1:174962/1‑139 (MQ=255)
           cGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCAc                                                                                                                 >  2:868315/1‑139 (MQ=255)
              atatAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCg                                                                                                              >  2:744184/1‑139 (MQ=255)
              atatAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCg                                                                                                              >  1:844678/1‑139 (MQ=255)
              atatAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCg                                                                                                              >  2:306910/1‑139 (MQ=255)
                       aGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCg                                                                                                     >  2:1360826/1‑139 (MQ=255)
                       aGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCg                                                                                                     >  1:306479/1‑139 (MQ=255)
                         ggACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGc                                                                                                   <  1:275530/139‑1 (MQ=255)
                             tGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTg                                                                                               >  1:848760/1‑139 (MQ=255)
                             tGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTg                                                                                               >  1:1339014/1‑139 (MQ=255)
                              gAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGt                                                                                              >  2:614916/1‑139 (MQ=255)
                               aaaCATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGtt                                                                                             >  1:1571860/1‑139 (MQ=255)
                                    ttattCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATg                                                                                        >  2:893592/1‑139 (MQ=255)
                                     tattCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAAcc                                                                                                      <  2:1116694/124‑1 (MQ=255)
                                     tattCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAAcc                                                                                                      >  1:1116694/1‑124 (MQ=255)
                                      attCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGc                                                                                      <  1:614916/139‑1 (MQ=255)
                                       ttCTTACGCGGTCTGACTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGcc                                                                                     >  2:316167/1‑139 (MQ=255)
                                        tCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCg                                                                                    <  1:1360826/139‑1 (MQ=255)
                                          ttACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGt                                                                                  <  2:174962/139‑1 (MQ=255)
                                            aCGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTgg                                                                                >  2:416939/1‑139 (MQ=255)
                                                 gTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACaaa                                                                           >  1:815243/1‑139 (MQ=255)
                                                 gTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACaaa                                                                           >  2:288363/1‑139 (MQ=255)
                                                   cTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAaaat                                                                         <  1:306910/139‑1 (MQ=255)
                                                    tGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAaaata                                                                        >  1:191052/1‑139 (MQ=255)
                                                        tGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                    <  2:1339014/139‑1 (MQ=255)
                                                        tGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATGTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                    <  1:1194920/139‑1 (MQ=255)
                                                            aTTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                    >  2:79896/1‑135 (MQ=255)
                                                            aTTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                    <  1:79896/135‑1 (MQ=255)
                                                             ttACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCt                                                               <  1:416939/139‑1 (MQ=255)
                                                                cTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCa                                                                                                                  <  2:227297/85‑1 (MQ=255)
                                                                cTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCa                                                                                                                  >  1:227297/1‑85 (MQ=255)
                                                                cTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCt                                                                                                 >  2:1565747/1‑102 (MQ=255)
                                                                cTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCt                                                                                                 <  1:1565747/102‑1 (MQ=255)
                                                                  cAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCa                                                          <  2:93185/139‑1 (MQ=255)
                                                                    aaaTCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATa                                                        <  2:844678/139‑1 (MQ=255)
                                                                          gTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCa                                                                                                                  <  2:575062/75‑1 (MQ=255)
                                                                          gTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCa                                                                                                                  >  1:575062/1‑75 (MQ=255)
                                                                          gTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgacg                                                  <  1:868315/139‑1 (MQ=255)
                                                                          gTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgacg                                                  <  2:1571860/139‑1 (MQ=255)
                                                                           ttGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGc                                                 <  1:1549368/139‑1 (MQ=255)
                                                                              cAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTGACTCGACGTCCTCACGCACTCGCCATATCCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAAt                                              <  1:1227284/139‑1 (MQ=255)
                                                                               aaaaGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAAtg                                             >  1:672961/1‑139 (MQ=255)
                                                                                   gTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGa                                         <  1:744184/139‑1 (MQ=255)
                                                                                        gAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGacacac                                    >  2:431478/1‑139 (MQ=255)
                                                                                         aGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACt                                   <  2:191052/139‑1 (MQ=255)
                                                                                           tCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATaa                                                                                                        >  1:1360112/1‑68 (MQ=255)
                                                                                           tCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATaa                                                                                                        <  2:1360112/68‑1 (MQ=255)
                                                                                           tCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTg                                 >  2:391317/1‑139 (MQ=255)
                                                                                           tCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTg                                 <  2:848760/139‑1 (MQ=255)
                                                                                            cATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCg                                                                                    <  1:898717/87‑1 (MQ=255)
                                                                                            cATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCg                                                                                    >  2:898717/1‑87 (MQ=255)
                                                                                            cATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTgg                                <  1:316167/139‑1 (MQ=255)
                                                                                              tttttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACg                                                                                                                            <  2:1555228/45‑1 (MQ=255)
                                                                                              tttttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACg                                                                                                                            >  1:1555228/1‑45 (MQ=255)
                                                                                                     tGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCg                                                                                                     >  2:1438203/1‑61 (MQ=255)
                                                                                                     tGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCg                                                                                                     <  1:1438203/61‑1 (MQ=255)
                                                                                                     tGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgac                                                   <  2:1155284/111‑1 (MQ=255)
                                                                                                     tGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACt                       >  1:1357928/1‑139 (MQ=255)
                                                                                                     tGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCAGGTGGACAAAATAAATTCTGCTTAGCATAacgac                                                   >  1:1155284/1‑111 (MQ=255)
                                                                                                      gATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACtt                      <  2:439325/139‑1 (MQ=255)
                                                                                                         tACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTg                   <  1:431478/139‑1 (MQ=255)
                                                                                                                 ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTg                   >  2:928827/1‑131 (MQ=255)
                                                                                                                 ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTg                   <  1:928827/131‑1 (MQ=255)
                                                                                                                  tACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTg                   <  1:584357/130‑1 (MQ=255)
                                                                                                                  tACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTg                   >  2:584357/1‑130 (MQ=255)
                                                                                                                        cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGa                                         >  2:926080/1‑102 (MQ=255)
                                                                                                                        cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGa                                         <  1:926080/102‑1 (MQ=255)
                                                                                                                         tCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCAt                                                                                         >  2:548209/1‑53 (MQ=255)
                                                                                                                         tCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCAt                                                                                         <  1:548209/53‑1 (MQ=255)
                                                                                                                         tCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgacg                                                  <  2:1164836/92‑1 (MQ=255)
                                                                                                                         tCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgacg                                                  >  1:1164836/1‑92 (MQ=255)
                                                                                                                          ccccAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                    <  2:727738/73‑1 (MQ=255)
                                                                                                                          ccccAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                    >  1:727738/1‑73 (MQ=255)
                                                                                                                          ccccAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCa  <  1:288363/139‑1 (MQ=255)
                                                                                                                           cccAACTAACTCGACGTCCTCACGTACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                    >  1:1094817/1‑72 (MQ=255)
                                                                                                                           cccAACTAACTCGACGTCCTCACGTACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                    <  2:1094817/72‑1 (MQ=255)
                                                                                                                             cAAGTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCt                                                                                                 <  2:1078417/41‑1 (MQ=255)
                                                                                                                             cAACTAACTCGACGTTCTCACGCACTCGCCATAACCGGCCt                                                                                                 >  1:1078417/1‑41 (MQ=255)
                                                                                                                                   actCGACGTCCTCACGCACTCGCCATAACCGGCCtt                                                                                                <  2:1171628/36‑1 (MQ=255)
                                                                                                                                   actCGACGTCCTCACGCACTCGCCATAACCGGCCtt                                                                                                >  1:1171628/1‑36 (MQ=255)
                                                                                                                                   actCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGc                                                 <  1:868289/83‑1 (MQ=255)
                                                                                                                                   actCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGc                                                 >  2:868289/1‑83 (MQ=255)
                                                                                                                                   actCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGaca                                       >  1:1570126/1‑93 (MQ=255)
                                                                                                                                   actCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGaca                                       <  2:1570126/93‑1 (MQ=255)
                                                                                                                                        |                                                                                                                            
CTACTATAAAGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCA  >  NC_000913/1200745‑1201005

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