Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 293 A→G intergenic (+38/‑44) thrL → / → thrA thr operon leader peptide/aspartokinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132930AG99.1% 336.6 / ‑4.8 110intergenic (+38/‑44)thrL/thrAthr operon leader peptide/aspartokinase I
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (38/71);  total (39/71)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

Only 100 of 119 total aligned reads displayed.

GCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAA  >  NC_000913/203‑381
                                                                                          |                                                                                         
gCACCACTATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGt                                                                                         <  1:2616408/93‑1 (MQ=255)
gCACCACCATCACCCCCACCATCACCATTCCCACAGGTAACGGTGCGGGCTGCCGCGTACAGAAAACACAGAAAAAAGCCCGCACCTGACGGt                                                                                         <  2:3389439/93‑1 (MQ=255)
  accaccATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGc                                                                                       <  1:4171383/93‑1 (MQ=255)
    caccaTCACCCCCACCACCACCATTCCCCAGGGTAACGGTGCGGGCTGACGCGTACAGGAAACAAAGAAAAAAGCCCGCACCTGACGGTGCgg                                                                                     <  1:1672786/93‑1 (MQ=255)
     accaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCggg                                                                                    >  2:2710060/1‑93 (MQ=255)
        aTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCtt                                                                                 >  1:392734/1‑93 (MQ=255)
        aTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCtt                                                                                 <  2:2715749/93‑1 (MQ=255)
        aTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCtt                                                                                 >  1:2370323/1‑93 (MQ=255)
         ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCttt                                                                                >  2:484743/1‑93 (MQ=255)
         ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCttt                                                                                >  2:2624136/1‑93 (MQ=255)
          tACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCtttt                                                                               <  2:602526/93‑1 (MQ=255)
          tACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCtttt                                                                               >  1:746801/1‑93 (MQ=255)
           tccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCttttt                                                                              <  2:4265858/92‑1 (MQ=255)
             caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCttttttt                                                                            >  1:4338833/1‑93 (MQ=255)
              accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCtttttttt                                                                           <  2:2164230/93‑1 (MQ=255)
              accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCCGACGGGGCGGGCtttttttt                                                                           >  2:927627/1‑93 (MQ=255)
               ccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCCCCTGACGGTGCGGGCttttttttc                                                                          >  1:1542469/1‑92 (MQ=255)
               ccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCttttttttt                                                                          >  2:352181/1‑93 (MQ=255)
               ccaccaCCCCCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCttttttttt                                                                          <  2:918100/93‑1 (MQ=255)
                 accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCg                                                                        >  2:654855/1‑93 (MQ=255)
                  ccatccccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGa                                                                       <  2:924858/93‑1 (MQ=255)
                  ccaccaccatTCCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGa                                                                       <  2:958071/93‑1 (MQ=255)
                   catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGGTTTTTTTTTTGAc                                                                     >  2:1695208/1‑93 (MQ=255)
                   catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGAc                                                                     >  1:2761002/1‑93 (MQ=255)
                     ccacaatAACCACAGGTAAGGGTGCGGGCTGACGCGAACAGAAAACAAAGAAAAAAGCCCGCCCCTGACGGTGCGGGCTTTTTTTTTCGACCa                                                                   <  2:1637736/92‑1 (MQ=255)
                       cccatCACCACAGGTAACGGGGGGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCGGCACCTGACGGTGCGGGCTTTTTTTTTCGACCaaa                                                                 <  1:3187849/92‑1 (MQ=255)
                        ccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAg                                                                >  1:2885097/1‑93 (MQ=255)
                        ccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAg                                                                <  1:3414050/93‑1 (MQ=255)
                         catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAgg                                                               >  1:3786849/1‑93 (MQ=255)
                          atTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTTGACCCAAAgg                                                               >  1:2724918/1‑93 (MQ=255)
                          atTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGt                                                              <  1:4192081/93‑1 (MQ=255)
                           ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTa                                                             <  2:565631/93‑1 (MQ=255)
                           ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTa                                                             <  1:1094201/93‑1 (MQ=255)
                            tcccACAGGTACCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTaa                                                            <  1:2536682/91‑1 (MQ=255)
                            tACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTaa                                                            <  1:3143339/93‑1 (MQ=255)
                             aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAAc                                                           <  1:2374567/93‑1 (MQ=255)
                               cacaGGTAACGGTGCGGGCTGACGCTTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGa                                                         <  1:1820772/93‑1 (MQ=255)
                                acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGCCCAAAGGTAACGAg                                                        >  1:1088731/1‑93 (MQ=255)
                                acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAg                                                        <  2:3646765/93‑1 (MQ=255)
                                  aGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTTGACCAAAAGGAACGAGGt                                                      >  2:4123231/1‑93 (MQ=255)
                                  aGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAggg                                                      >  1:1866315/1‑92 (MQ=255)
                                    gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaa                                                    >  1:322044/1‑93 (MQ=255)
                                     taccGGTGCGGTCTGACCCGTACAGAAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaac                                                   <  2:4352587/90‑1 (MQ=255)
                                     tAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaac                                                   >  2:4635255/1‑93 (MQ=255)
                                       aCGGTGGGGGCTGACGCGTACAGAAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                 <  1:3516239/93‑1 (MQ=255)
                                        cGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaac                                                >  1:699793/1‑93 (MQ=255)
                                         ggTGCGGGCTGACGCGTAAAGGAAACACAGAAAAAACCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAAcc                                               <  2:4154548/93‑1 (MQ=255)
                                          gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                              <  2:2880840/93‑1 (MQ=255)
                                           tGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                             <  1:4283724/93‑1 (MQ=255)
                                            gCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATg                                            <  1:1984902/93‑1 (MQ=255)
                                             cGGGCTGACGCATAAAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                           <  1:3272062/93‑1 (MQ=255)
                                               gtcTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGa                                         <  1:3624344/91‑1 (MQ=255)
                                               ggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCCTGCGa                                         >  1:1184564/1‑93 (MQ=255)
                                               ggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGa                                         <  1:2625449/93‑1 (MQ=255)
                                                gCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTTGACCAAAGGTAACGAGGTAACAACCATGCGAg                                        >  1:2696075/1‑93 (MQ=255)
                                                gCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                                        <  1:4003006/93‑1 (MQ=255)
                                                 cTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgt                                       <  2:1510907/93‑1 (MQ=255)
                                                      gcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGa                                  <  1:2448833/93‑1 (MQ=255)
                                                      gcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGAGGGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGGGGTGa                                  >  2:1241290/1‑93 (MQ=255)
                                                      gcgTACAGAAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGa                                  <  1:938048/93‑1 (MQ=255)
                                                       cgTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                 <  2:1561405/93‑1 (MQ=255)
                                                         tACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGt                               <  2:1617524/93‑1 (MQ=255)
                                                           cAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                             <  2:1988745/93‑1 (MQ=255)
                                                           cAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                             >  2:2895934/1‑93 (MQ=255)
                                                           cAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                             <  2:4588225/93‑1 (MQ=255)
                                                           aaGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                             <  1:3352384/92‑1 (MQ=255)
                                                            aagAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                            >  2:3733572/3‑93 (MQ=255)
                                                            aGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                            >  2:4312847/1‑93 (MQ=255)
                                                             gtaaaCAAAGAAAAAAGCCGGCACCTGACGGTGCGTGCTTTTTTTTCGACCAAAGATAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                          <  1:3077877/91‑1 (MQ=255)
                                                             ggAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                           <  1:2986580/93‑1 (MQ=255)
                                                                 acacAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGt                       <  1:652304/93‑1 (MQ=255)
                                                                 acacAGAAAAAAGCCCGCACCTGACGGGGCGGGCTTTTTTTTTTGACCAAAGGGAACGAGGTAACAACCATGCGAGTGTTGAAGGTCGGCGGt                       >  2:2718208/1‑93 (MQ=255)
                                                                  cacaGAAAAAGGCCGGAACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                      <  1:4300753/93‑1 (MQ=255)
                                                                  cacaAAAAAAACCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                      <  2:4588609/93‑1 (MQ=255)
                                                                     aGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAt                   <  1:1398476/93‑1 (MQ=255)
                                                                      aaaaaaaGCCCGCCCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                  <  1:1155111/92‑1 (MQ=255)
                                                                      aaaaaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                  <  1:270484/92‑1 (MQ=255)
                                                                         aaaaGCCCGCCCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt               <  2:1893552/93‑1 (MQ=255)
                                                                          agaccccGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg              <  2:1769374/89‑1 (MQ=255)
                                                                          aaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg              <  2:516626/93‑1 (MQ=255)
                                                                          aaaGCCCACACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg              <  1:637424/93‑1 (MQ=255)
                                                                            aGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc            <  2:3140725/93‑1 (MQ=255)
                                                                               ccGCCCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaaa         <  2:3307484/93‑1 (MQ=255)
                                                                                cGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAt        <  1:3246075/93‑1 (MQ=255)
                                                                                cGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAt        <  1:3132171/93‑1 (MQ=255)
                                                                                cGCACCTGACGGTGCGGGCTTTTATTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAt        <  1:1075140/93‑1 (MQ=255)
                                                                                 gCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATg       >  1:74381/1‑93 (MQ=255)
                                                                                 gCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATg       <  1:4200638/93‑1 (MQ=255)
                                                                                  cACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc      >  2:1693104/1‑93 (MQ=255)
                                                                                   aCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa     >  1:2916909/1‑93 (MQ=255)
                                                                                   aCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa     >  1:992855/1‑93 (MQ=255)
                                                                                   aCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGGGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa     >  1:3046151/1‑93 (MQ=255)
                                                                                    ccTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg    >  2:2754448/1‑93 (MQ=255)
                                                                                    ccTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg    <  2:2186397/93‑1 (MQ=255)
                                                                                    ccTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg    <  1:4173478/93‑1 (MQ=255)
                                                                                     cTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGa   >  1:1845674/1‑93 (MQ=255)
                                                                                     cTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGAGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGa   >  2:3256523/1‑93 (MQ=255)
                                                                                      tGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa  <  1:2622718/93‑1 (MQ=255)
                                                                                      tGACGGTGCGGGCTTTTTTTTTCGACCAAAGGGAACGAGGTAACAACCATGCGAGTGCTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa  >  1:2247421/1‑93 (MQ=255)
                                                                                      tGACGGTGCGCGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa  <  1:676716/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAA  >  NC_000913/203‑381

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