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.0% 295.2 / ‑4.7 95intergenic (+38/‑44)thrL/thrAthr operon leader peptide/aspartokinase I
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (35/59);  total (36/59)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 101 total aligned reads displayed.

GCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAC  >  NC_000913/203‑382
                                                                                          |                                                                                           
gCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGt                                                                                           <  2:2706770/93‑1 (MQ=255)
gCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGt                                                                                           <  1:4484505/93‑1 (MQ=255)
 caccacCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTg                                                                                          <  2:2702887/93‑1 (MQ=255)
  accaccATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGc                                                                                         <  2:233718/93‑1 (MQ=255)
    caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCgg                                                                                       <  2:1292361/93‑1 (MQ=255)
    caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCgg                                                                                       <  2:4445743/93‑1 (MQ=255)
    caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCgg                                                                                       <  1:3899849/93‑1 (MQ=255)
     accaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCggg                                                                                      >  2:4514423/1‑93 (MQ=255)
       caatACCACCACCCTCCCCATTCCCACAGGTAACGGTGCGGGCTGGCGCAAACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCt                                                                                   <  2:3685431/90‑1 (MQ=255)
           accacccccATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCttttt                                                                              <  2:1042602/93‑1 (MQ=255)
            ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCtttttt                                                                             <  1:2210132/93‑1 (MQ=255)
            ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGCGCGGGGtttttt                                                                             >  1:2601192/1‑93 (MQ=255)
             caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCttttttt                                                                            <  1:616189/93‑1 (MQ=255)
              accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCtttttttt                                                                           >  1:3284868/1‑93 (MQ=255)
               ccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGTAACACAGAACAAAGCCAGCACCTGACGGTGCGGGCttttttttt                                                                          >  1:3973842/1‑93 (MQ=255)
                    atcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCCGACGGGGCGGGGGTTTTTTTTCGAcc                                                                     >  2:1959088/1‑93 (MQ=255)
                     tcaacatTACCACAAGCAACGGTGCGGGCTGACGCGTACAGGAAACACAGGAAAAAGCCCGCACCTGACGGCGCGGGCTTTTATTTTTCGAcc                                                                     >  2:1406765/1‑93 (MQ=255)
                      caccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCaa                                                                   <  1:1104219/93‑1 (MQ=255)
                       accatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTTGACCcaa                                                                  >  1:1114146/1‑93 (MQ=255)
                           ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTa                                                              >  2:600524/1‑93 (MQ=255)
                            tACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTaa                                                             <  1:119903/93‑1 (MQ=255)
                              ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACg                                                           >  2:819431/1‑93 (MQ=255)
                                acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTTTGACCAAAGGGAACGa                                                          >  2:3157977/1‑93 (MQ=255)
                                acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAg                                                         <  1:784584/93‑1 (MQ=255)
                                acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAg                                                         <  1:94747/93‑1 (MQ=255)
                                 caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAgg                                                        <  1:4098523/93‑1 (MQ=255)
                                 caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTAGCCCAAAGGTAACGCgg                                                        >  2:3714186/1‑93 (MQ=255)
                                  aGGTATCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCCCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGt                                                       <  1:2728472/93‑1 (MQ=255)
                                  aGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCCCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGt                                                       <  2:1644338/93‑1 (MQ=255)
                                     tAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaac                                                    <  2:1941503/93‑1 (MQ=255)
                                     tAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTCGACCAAAGGTAACGAGGTaaca                                                   >  2:1673116/1‑93 (MQ=255)
                                      aaCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaaca                                                   <  2:416791/93‑1 (MQ=255)
                                       aCGGTGCGGGCTGACGGGTAAAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                  <  1:2800287/93‑1 (MQ=255)
                                       aCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                  >  2:2899415/1‑93 (MQ=255)
                                        cGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaac                                                 <  2:465711/93‑1 (MQ=255)
                                         ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAAcc                                                <  1:1321492/93‑1 (MQ=255)
                                         ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAAcc                                                <  2:1571306/93‑1 (MQ=255)
                                           tGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                              <  1:3872318/93‑1 (MQ=255)
                                            ggggggTGACGCGTACAGAAAAAAAAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATg                                             <  2:54888/91‑1 (MQ=255)
                                             cGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                            <  1:1839276/93‑1 (MQ=255)
                                                 cTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgt                                        <  2:1247398/93‑1 (MQ=255)
                                                  tGACGCGTAAAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgtg                                       >  1:4674875/1‑93 (MQ=255)
                                                   gACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgtgt                                      >  2:1233650/1‑93 (MQ=255)
                                                    aCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                     <  2:1315502/93‑1 (MQ=255)
                                                      gcgAACAGGAAACACAGAAAAACCCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGa                                   <  2:2565641/93‑1 (MQ=255)
                                                         tACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACACCCATGCGAGTGTTGAAGt                                <  1:4611837/93‑1 (MQ=255)
                                                            aGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAATGAGGTAACAACCATGCGAGGGTTGAAGTTcg                             >  1:3999970/1‑93 (MQ=255)
                                                            aGGAAACACAGAAAAAAGCCCCCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGGAACGAGGCAACAAACATGCGAGTGTTGAAGTTcg                             >  1:4346659/1‑93 (MQ=255)
                                                              gAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                           <  1:1080681/93‑1 (MQ=255)
                                                                aaCACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcgg                         <  2:4329057/93‑1 (MQ=255)
                                                                aaCACAGAAAAAAGCCCGCACCCGACCGTGCGGGCTTTTTTTTTTGACCTACTGGAATGAAGTCACAAACCTGCGAGCGGTGAGGTTcggctg                         >  2:2801476/1‑91 (MQ=255)
                                                                 acacAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGt                        >  1:2143376/1‑93 (MQ=255)
                                                                 acacAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGt                        <  2:2144353/93‑1 (MQ=255)
                                                                   acaGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCACAGGTACCGAGGAAACAACCATGGGGGGGTTAAACTTCGCGGTAcc                     >  1:4095816/1‑92 (MQ=255)
                                                                    caGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACa                     <  1:4688090/93‑1 (MQ=255)
                                                                       aaaaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCa                  <  1:3732863/93‑1 (MQ=255)
                                                                        aaaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAg                 <  2:136833/93‑1 (MQ=255)
                                                                        aaaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAg                 <  1:1813470/93‑1 (MQ=255)
                                                                        aaaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAg                 <  2:4421020/93‑1 (MQ=255)
                                                                         aaaaGCCCGCCCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                <  1:1521095/93‑1 (MQ=255)
                                                                         aaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                <  2:1573462/93‑1 (MQ=255)
                                                                         aaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                <  1:872350/93‑1 (MQ=255)
                                                                          aaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg               <  2:3396806/93‑1 (MQ=255)
                                                                          aaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGGACATCAGTg               >  2:3101745/1‑93 (MQ=255)
                                                                             gCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGGAACGAGGTAACAACCATGCGAGGGTTGAAGTTCGGCGGTACATCAGTGGCa            >  2:3870575/1‑93 (MQ=255)
                                                                              cccGCACCTGACAGTGCGGGCTTTTTTTTTTGACCAAAGGTAACGAGGGAACAACCATGCGCGTGTTGAAGTTCGGCGGTACATCAGTGGCaa           >  2:1993906/1‑93 (MQ=255)
                                                                               ccGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaaa          >  1:4526568/1‑93 (MQ=255)
                                                                               ccGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaaa          >  2:4043180/1‑93 (MQ=255)
                                                                                cGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAt         >  2:2195017/1‑93 (MQ=255)
                                                                                cGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAt         <  1:2316172/93‑1 (MQ=255)
                                                                                 gCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATg        <  2:1973159/93‑1 (MQ=255)
                                                                                 gCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATg        >  1:1272710/1‑93 (MQ=255)
                                                                                 gCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATg        <  1:3125712/93‑1 (MQ=255)
                                                                                 gCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATg        >  2:2952034/1‑93 (MQ=255)
                                                                                  cACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATga       >  1:389270/1‑92 (MQ=255)
                                                                                  cACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc       <  1:4367497/93‑1 (MQ=255)
                                                                                  cACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc       >  1:2803228/1‑93 (MQ=255)
                                                                                  cACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCAAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc       >  2:3465748/1‑93 (MQ=255)
                                                                                   cccTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa      <  1:1512487/92‑1 (MQ=255)
                                                                                   aCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa      >  1:1172277/1‑93 (MQ=255)
                                                                                   aCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa      >  1:3627014/1‑93 (MQ=255)
                                                                                   aCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa      <  2:3073406/93‑1 (MQ=255)
                                                                                    ccTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCaa     >  2:693896/1‑92 (MQ=255)
                                                                                    ccTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg     <  2:2288707/93‑1 (MQ=255)
                                                                                    ccTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg     <  2:89260/93‑1 (MQ=255)
                                                                                    ccTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg     >  2:4301216/1‑93 (MQ=255)
                                                                                    ccTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg     >  1:3488902/1‑93 (MQ=255)
                                                                                    ccTGACGGCGAGGGATTTTTTTTTCGGCCAAAAGCAACGAGGCAACAACCATGCGAGTGTTCAAGTGCTGCAGGACATCAGTTGCAAATGCAg     >  2:1715483/1‑93 (MQ=255)
                                                                                     cTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGa    >  1:3657994/1‑93 (MQ=255)
                                                                                     cTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGa    <  2:2314591/93‑1 (MQ=255)
                                                                                     cTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGa    <  2:3182494/93‑1 (MQ=255)
                                                                                     cTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGa    >  2:4513876/1‑93 (MQ=255)
                                                                                     cTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGa    <  2:1135372/93‑1 (MQ=255)
                                                                                     cTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGa    >  1:3104543/1‑93 (MQ=255)
                                                                                     cTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGa    >  1:79740/1‑93 (MQ=255)
                                                                                      tGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa   <  2:3787858/93‑1 (MQ=255)
                                                                                      tGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa   <  1:4096266/93‑1 (MQ=255)
                                                                                      tGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa   <  1:4467841/93‑1 (MQ=255)
                                                                                      tGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACACCAGTGGCAAATGCAGaa   <  1:1553260/93‑1 (MQ=255)
                                                                                       gACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc  <  2:3171430/93‑1 (MQ=255)
                                                                                          |                                                                                           
GCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAC  >  NC_000913/203‑382

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