Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 296,697 A→G F107L (TTC→CTC)  yagP ← putative transcriptional regulator LYSR‑type

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913296,6970AG99.3% 981.5 / ‑4.2 286F107L (TTC→CTC) yagPputative transcriptional regulator LYSR‑type
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (140/143);  total (142/144)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.12e-01

Only 100 of 293 total aligned reads displayed.

AAGCCAATTGCTGCCGCCATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAAGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGAGTGGCAATTG  >  NC_000913/296607‑296786
                                                                                          |                                                                                         
 aGCCAATTGCTGCCGCCATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGcc                                                                                        <  1:4216882/93‑1 (MQ=255)
 aGCCAATTGCTGCCGCCATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGcc                                                                                        <  2:8031949/93‑1 (MQ=255)
  gCCAATTGCTGCCGCCATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGccc                                                                                       <  1:5681452/93‑1 (MQ=255)
      aTTGCTGCCGCCATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGaa                                                                                   >  2:7107076/1‑93 (MQ=255)
      aTTGCTGCCGCCATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGaa                                                                                   >  1:5818666/1‑93 (MQ=255)
       ttGCTGCCGCCATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGaaa                                                                                  <  2:5929400/93‑1 (MQ=255)
           tgcCGCCATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATtgt                                                                              <  1:4109788/93‑1 (MQ=255)
             ccgccATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTgg                                                                            >  1:8789881/1‑93 (MQ=255)
             ccgccATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTgg                                                                            >  2:5252658/1‑93 (MQ=255)
             ccgccATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTgg                                                                            <  2:6207284/93‑1 (MQ=255)
              cgccATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGga                                                                           <  2:5630980/93‑1 (MQ=255)
              cgccATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGga                                                                           <  2:8760616/93‑1 (MQ=255)
              cgccATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGga                                                                           <  1:4353939/93‑1 (MQ=255)
               gccATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGgaa                                                                          <  2:7502651/93‑1 (MQ=255)
                ccATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGgaag                                                                         <  2:304389/93‑1 (MQ=255)
                ccATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGgaag                                                                         <  2:8146422/93‑1 (MQ=255)
                 cATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGgaaga                                                                        <  1:6654311/93‑1 (MQ=255)
                 cATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGAGgaaga                                                                        >  1:4561375/1‑93 (MQ=255)
                 cATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGAGGGAAATTCAGATAGAGGCCCGGAAATTGTGgaaga                                                                        >  1:6697332/1‑93 (MQ=255)
                  aTTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGCAATTGTGgaagaa                                                                       >  2:3752597/1‑93 (MQ=255)
                   ttCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGaaa                                                                      <  2:2980045/93‑1 (MQ=255)
                   ttCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGaaa                                                                      <  2:5560521/93‑1 (MQ=255)
                     cTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGt                                                                    <  1:7618638/93‑1 (MQ=255)
                     cTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGt                                                                    <  2:6546468/93‑1 (MQ=255)
                       ttGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGtca                                                                  >  1:5000385/1‑93 (MQ=255)
                        tGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAAGTGTGGAAGAAAGtcat                                                                 >  1:7787866/1‑93 (MQ=255)
                         gACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGtcatc                                                                <  1:6255716/93‑1 (MQ=255)
                         gACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGtcatc                                                                >  1:1630151/1‑93 (MQ=255)
                           cGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCaa                                                              >  1:3936159/1‑93 (MQ=255)
                           cGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCaa                                                              <  1:2283583/93‑1 (MQ=255)
                            gTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAag                                                             >  2:1176283/1‑93 (MQ=255)
                             tAGTCAATCAGGGCGCAGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAagc                                                            <  1:1956676/93‑1 (MQ=255)
                              aGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAagca                                                           >  1:7218457/1‑93 (MQ=255)
                               gTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAagcag                                                          >  1:1064944/1‑93 (MQ=255)
                               gTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAagcag                                                          >  1:3377756/1‑93 (MQ=255)
                                 cAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCg                                                        <  1:2228573/93‑1 (MQ=255)
                                  aaTCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGTAGAAAGTCATCAAGCAGCGa                                                       >  2:8367354/1‑93 (MQ=255)
                                  aaTCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGa                                                       <  1:7944733/93‑1 (MQ=255)
                                   aTCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGAt                                                      >  1:1880580/1‑93 (MQ=255)
                                      aGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACa                                                   <  2:1817998/93‑1 (MQ=255)
                                      aGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACa                                                   >  2:3446513/1‑93 (MQ=255)
                                      aGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACa                                                   <  2:1442759/93‑1 (MQ=255)
                                      aGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACa                                                   >  1:8202665/1‑93 (MQ=255)
                                       gggCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACaa                                                  >  1:978295/1‑93 (MQ=255)
                                         gcgcGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGc                                                <  2:4056973/93‑1 (MQ=255)
                                         gcgcGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGc                                                >  2:4784188/1‑93 (MQ=255)
                                          ggcgGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCt                                               <  2:6796004/92‑1 (MQ=255)
                                           gcgGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCtt                                              <  1:2137080/93‑1 (MQ=255)
                                           gcgGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGAtt                                              >  2:790627/1‑93 (MQ=255)
                                            cgGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTa                                             >  2:8819085/1‑93 (MQ=255)
                                            cgGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTa                                             <  2:2363497/93‑1 (MQ=255)
                                              gAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTAcc                                           <  1:3515110/93‑1 (MQ=255)
                                              gAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTAcc                                           >  1:1376347/1‑93 (MQ=255)
                                              gAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTAcc                                           <  2:7457901/93‑1 (MQ=255)
                                               aGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCg                                          <  1:1743365/93‑1 (MQ=255)
                                                gCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGc                                         >  2:5777688/1‑93 (MQ=255)
                                                gCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGc                                         <  1:8548051/93‑1 (MQ=255)
                                                  tttGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCtt                                       >  2:603398/1‑93 (MQ=255)
                                                  tttGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCtt                                       <  2:784019/93‑1 (MQ=255)
                                                    tGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTc                                     >  2:8472834/1‑93 (MQ=255)
                                                     ggTGCAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCa                                    <  1:2405364/93‑1 (MQ=255)
                                                       ggCAATATTGCGACGCTCTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAat                                  <  2:432807/92‑1 (MQ=255)
                                                         cAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAatat                                <  1:8348561/93‑1 (MQ=255)
                                                         cAATATTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAatat                                >  1:8613999/1‑93 (MQ=255)
                                                           atatTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAatatat                              <  2:3353839/93‑1 (MQ=255)
                                                           atatTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAatatat                              <  1:543195/93‑1 (MQ=255)
                                                            tatTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATg                             >  2:2082428/1‑93 (MQ=255)
                                                             atTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATgg                            <  1:5234702/93‑1 (MQ=255)
                                                             atTGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATgg                            <  2:1551794/93‑1 (MQ=255)
                                                              ttGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGc                           >  2:6690109/1‑93 (MQ=255)
                                                              ttGCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGc                           <  2:2961667/93‑1 (MQ=255)
                                                                gCGACGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCt                         <  2:281918/93‑1 (MQ=255)
                                                                   aCGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGaa                      <  2:6961448/93‑1 (MQ=255)
                                                                    cGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGaaa                     <  2:3985707/93‑1 (MQ=255)
                                                                    cGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGaaa                     >  2:5723931/1‑93 (MQ=255)
                                                                    cGCTGTGGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGaaa                     <  1:5310984/93‑1 (MQ=255)
                                                                     gCTGTGGGAAATACAGATAGAGGCCCGGAAATCGTGGAAGAAATTCCTAAAGCAGGGATACAAGCTTACAGCTTTCACTATATGGCCGGaaat                    >  1:7563967/1‑92 (MQ=255)
                                                                        gtgGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGttt                 >  1:4634389/1‑93 (MQ=255)
                                                                        gtgGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGttt                 >  1:612571/1‑93 (MQ=255)
                                                                        gtgGGAAATACAGATAGAAGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGttt                 >  2:7241814/1‑93 (MQ=255)
                                                                         tgGGAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTc                >  2:8620308/1‑93 (MQ=255)
                                                                          gggAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTcc               <  1:6025771/93‑1 (MQ=255)
                                                                           ggAAATACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCt              <  2:7946595/93‑1 (MQ=255)
                                                                             aaaTACAGAGAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGa            <  2:2403098/93‑1 (MQ=255)
                                                                              aaTACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGAg           >  2:4551397/1‑93 (MQ=255)
                                                                              aaTACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGAg           >  1:5935437/1‑93 (MQ=255)
                                                                               aTACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATTCAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGAgg          >  1:1109809/1‑92 (MQ=255)
                                                                               aTACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGAGt          >  2:3210350/1‑93 (MQ=255)
                                                                               aTACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGAGt          <  1:8431578/93‑1 (MQ=255)
                                                                                tACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGGGTg         >  1:2890256/1‑93 (MQ=255)
                                                                                tACAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGAGTg         >  1:4029073/1‑93 (MQ=255)
                                                                                 aCAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGAGTgg        >  2:307834/1‑93 (MQ=255)
                                                                                  cAGATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGAGTGGc       <  1:3799871/93‑1 (MQ=255)
                                                                                    gATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGAGTGGCaa     >  2:5018383/1‑93 (MQ=255)
                                                                                    gATAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGAGTGGCaa     >  1:2040724/1‑93 (MQ=255)
                                                                                      tAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGAGTGGCAAtt   >  2:4182060/1‑93 (MQ=255)
                                                                                      tAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGAGTGGCAAtt   <  2:2682919/93‑1 (MQ=255)
                                                                                      tAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGAGTGGCAAtt   <  1:3159383/93‑1 (MQ=255)
                                                                                      tAGAGGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGAGTGGCAAtt   >  2:905071/1‑93 (MQ=255)
                                                                                       agagGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGAGTGGCAATTg  <  2:7633052/93‑1 (MQ=255)
                                                                                          |                                                                                         
AAGCCAATTGCTGCCGCCATTCTTTGACGTAGTCAATCAGGGCGCGGAGCTTTGGTGCAATATTGCGACGCTGTGGGAAATACAGATAGAAGCCCGGAAATTGTGGAAGAAAGTCATCAAGCAGCGATACAAGCTTACCGCTTTCAATATATGGCCTGAAAGTTTCCTGAGTGGCAATTG  >  NC_000913/296607‑296786

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