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.4% 513.5 / NA 161intergenic (+38/‑44)thrL/thrAthr operon leader peptide/aspartokinase I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (66/94);  minor base C (1/0);  total (67/94)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 178 total aligned reads displayed.

GCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAC  >  NC_000913/203‑382
                                                                                          |                                                                                         
  accaccATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGc                                                                                       <  2:6408434/93‑1 (MQ=255)
  accaccATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGc                                                                                       <  1:962200/93‑1 (MQ=255)
  accaccATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGc                                                                                       >  1:6721511/1‑93 (MQ=255)
  accaccATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGc                                                                                       <  1:5338172/93‑1 (MQ=255)
    caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCgg                                                                                     <  2:3126348/93‑1 (MQ=255)
    caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCgg                                                                                     >  1:4335614/1‑93 (MQ=255)
    caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCgg                                                                                     <  1:2706954/93‑1 (MQ=255)
       caTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCt                                                                                  <  2:3347824/93‑1 (MQ=255)
        aTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCgggtct                                                                                 >  2:5037213/1‑90 (MQ=255)
        aTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCtt                                                                                 >  2:2780974/1‑93 (MQ=255)
        aTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCtt                                                                                 <  1:2869642/93‑1 (MQ=255)
         ttACCACCACCATCACCAGTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCttt                                                                                <  1:2606459/93‑1 (MQ=255)
           accaccaccATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCttttt                                                                              >  1:890876/1‑93 (MQ=255)
           accaccaccATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCttttt                                                                              >  1:6387009/1‑93 (MQ=255)
            cccccaccaTCCCCATTACCACAGGTACCGGTGCGGGCTGACGCGTACGGGAAACACAGAAAAAAGCCCTCACCTGACGGTGCGGGCtttttt                                                                             <  1:448484/90‑1 (MQ=255)
             caccaccaTCACCATTACCACAGGTAGCGGTGCGGGCTGACGCGTACAGGAAACAAAGAAAAAAGCCCGCACCTGACGGTGCGGGCttttttt                                                                            <  1:5770343/93‑1 (MQ=255)
                caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTc                                                                         <  1:3305214/93‑1 (MQ=255)
                 accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCg                                                                        >  1:4261105/1‑93 (MQ=255)
                  ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGa                                                                       <  1:5275087/93‑1 (MQ=255)
                  ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGa                                                                       >  1:1071035/1‑93 (MQ=255)
                  ccatcaccatACCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGa                                                                       <  1:2481958/93‑1 (MQ=255)
                   catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGGTTTTTTTTTTGAc                                                                      >  2:1245547/1‑93 (MQ=255)
                    atcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGAcc                                                                     >  1:4113598/1‑93 (MQ=255)
                     tcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCCGACGGTGCGGGCTTTTTTTTTCGACCa                                                                    >  1:2501549/1‑93 (MQ=255)
                       cccatTACCACAGGTACCGGTGCGGGCTGACGCGTACAGGAACCACAAAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCaaa                                                                  <  2:4600664/92‑1 (MQ=255)
                        ccatTACCACAGGTAGCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAg                                                                 <  2:3968156/93‑1 (MQ=255)
                        ccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGGCGGTGCGGGCTTTTTTTTTTGATCAAAg                                                                 >  1:2398397/1‑93 (MQ=255)
                          atTCCCACAGGTAGCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAACGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGt                                                               <  2:2330366/93‑1 (MQ=255)
                          atTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGt                                                               >  2:2922710/1‑93 (MQ=255)
                          atTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGC‑‑TTTTTTTAGACCAAAGGTaa                                                             >  2:5224840/1‑93 (MQ=255)
                           ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTa                                                              <  2:5769726/93‑1 (MQ=255)
                           tccccACAGGTAACGGTGCGGGCTGACGCTTACAGAAAACACAGAAAAAACCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTa                                                              <  1:5784428/90‑1 (MQ=255)
                             aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAAc                                                            <  2:2575772/93‑1 (MQ=255)
                                acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAg                                                         <  2:3183224/93‑1 (MQ=255)
                                acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAg                                                         <  2:801058/93‑1 (MQ=255)
                                acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAg                                                         <  1:5573297/93‑1 (MQ=255)
                                  aGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGTACAAAAGTAACGAggg                                                       >  2:1199740/1‑92 (MQ=255)
                                   ggTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTa                                                      <  1:3022168/93‑1 (MQ=255)
                                     tAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaac                                                    <  1:6242607/93‑1 (MQ=255)
                                         ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAAcc                                                <  2:2869119/93‑1 (MQ=255)
                                         ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAAcc                                                <  2:1929695/93‑1 (MQ=255)
                                         ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAAcc                                                >  2:1811143/1‑93 (MQ=255)
                                          gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTTGACCAAAGGTAACGAGGTAAAAACCa                                               >  2:4835324/1‑93 (MQ=255)
                                          gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                               >  1:1406710/1‑93 (MQ=255)
                                           tGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                              <  2:1092877/93‑1 (MQ=255)
                                           ggggggCCGCCGCGTAAAGAAAACAAAAAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                              <  1:4588447/90‑1 (MQ=255)
                                             cGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                            <  1:3383393/93‑1 (MQ=255)
                                              gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTTGACCAAAGGTAACGAGGGAAAAAACAAGCg                                           >  2:20381/1‑93 (MQ=255)
                                               ggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGa                                          <  1:445073/93‑1 (MQ=255)
                                                  tGACGCCTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgtg                                       <  2:3024353/93‑1 (MQ=255)
                                                    acacgTACAGAAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                     <  1:4959849/90‑1 (MQ=255)
                                                    aCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                     <  1:6660452/93‑1 (MQ=255)
                                                    aCGCGTACAGAAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                     <  2:2523548/93‑1 (MQ=255)
                                                      gcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGa                                   >  1:1123629/1‑93 (MQ=255)
                                                       cttACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                  <  2:2831736/91‑1 (MQ=255)
                                                       cgTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                  >  2:4490745/1‑93 (MQ=255)
                                                       cgTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                  <  1:4232536/93‑1 (MQ=255)
                                                         tACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCACGCGAGTGTTGAAGt                                >  1:2022700/1‑93 (MQ=255)
                                                          aCAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGtt                               >  1:6620961/1‑93 (MQ=255)
                                                           cAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                              <  1:4369940/93‑1 (MQ=255)
                                                            aGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                             <  2:4365816/93‑1 (MQ=255)
                                                            aGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGGAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                             >  1:2201490/1‑93 (MQ=255)
                                                             ggAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                            >  2:6662496/1‑93 (MQ=255)
                                                              gAAACACAGAAAAAAGCCCGCTCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                           <  1:4837574/93‑1 (MQ=255)
                                                              gAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                           <  2:3811465/93‑1 (MQ=255)
                                                                 acacAGAAAAACGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGt                        <  1:1067144/93‑1 (MQ=255)
                                                                 acacAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGt                        <  1:6415783/93‑1 (MQ=255)
                                                                   acaGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAc                      <  2:6577476/93‑1 (MQ=255)
                                                                      aaaaaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                   <  2:1314519/92‑1 (MQ=255)
                                                                       aaaaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCa                  >  2:2027970/1‑93 (MQ=255)
                                                                        aaaaaGCCCGCCCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAg                 <  1:3299384/93‑1 (MQ=255)
                                                                        aaaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAg                 <  2:4238070/93‑1 (MQ=255)
                                                                        aaaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAg                 <  1:5282443/93‑1 (MQ=255)
                                                                         aaaaGCCCGCGCCGGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                <  2:5692468/93‑1 (MQ=255)
                                                                         aaaaGCCCGCACCTGACGGTGCGGGGTTTTTTTTTCGACCAAAGGGAACGAGGGAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                >  1:505754/1‑93 (MQ=255)
                                                                         aaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                <  1:5897741/93‑1 (MQ=255)
                                                                         aaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                <  2:2126048/93‑1 (MQ=255)
                                                                          aaaGCCCGCCCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg               <  1:6323508/93‑1 (MQ=255)
                                                                          aaaGCCCGCCCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg               <  2:4166117/93‑1 (MQ=255)
                                                                          aaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg               <  2:6694051/93‑1 (MQ=255)
                                                                          aaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg               <  2:6733165/93‑1 (MQ=255)
                                                                          aaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg               <  1:3768585/93‑1 (MQ=255)
                                                                            aGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc             >  2:896187/1‑93 (MQ=255)
                                                                             ccccGCCCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa            <  2:4675363/92‑1 (MQ=255)
                                                                               ccccACCTGACGGGGCGGGCTTTTTTTTTTGACCAAAGGTAACGAGGGAACAACCATGCGAGGGTTGAAGTTCGGCGGTACATCAGTCGCaaa          >  1:1559138/4‑93 (MQ=255)
                                                                               ccGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaaa          <  1:4653835/93‑1 (MQ=255)
                                                                                cGCCCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAt         <  2:5267979/93‑1 (MQ=255)
                                                                                cGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAt         <  2:3639172/93‑1 (MQ=255)
                                                                                cGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAt         >  2:5567508/1‑93 (MQ=255)
                                                                                 gCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATg        >  1:3226920/1‑93 (MQ=255)
                                                                                  cACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATga       >  2:5517769/1‑92 (MQ=255)
                                                                                   aCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa      <  2:5287560/93‑1 (MQ=255)
                                                                                   aCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa      >  1:983762/1‑93 (MQ=255)
                                                                                   aCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa      >  1:3603972/1‑93 (MQ=255)
                                                                                   aCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa      >  1:4150589/1‑93 (MQ=255)
                                                                                   aCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAATGACAAATGCa      >  2:6228326/1‑93 (MQ=255)
                                                                                    ccTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg     >  1:175620/1‑93 (MQ=255)
                                                                                     cTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGa    >  2:3867060/1‑93 (MQ=255)
                                                                                      tGACGGTGCGGGCTTTTTTTGTCGACCAAAGGTAACGAGGTAACAGCCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa   <  1:3948592/93‑1 (MQ=255)
                                                                                       gACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc  >  1:3850900/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAC  >  NC_000913/203‑382

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