BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	677,120	G→A	A161A (GCG→GCA)	ybeS →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	677,120	0	G	A	99.1%	810.3 / ‑4.2	216	A161A (GCG→GCA)	ybeS	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/2); new base A (121/93); total (121/95)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 224 total aligned reads displayed.

GTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCG  >  NC_000913/677030‑677209

tctattaGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCaaa                                                                                         <  1:2039359/91‑1 (MQ=255)

gtTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCaaa                                                                                         >  1:101963/1‑93 (MQ=255)

 ttattaGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAg                                                                                        >  2:4226390/1‑93 (MQ=255)

  tattaGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAgg                                                                                       >  2:4014611/1‑93 (MQ=255)

   attaGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAggg                                                                                      <  1:1113554/93‑1 (MQ=255)

   attaGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAggg                                                                                      <  1:3288636/93‑1 (MQ=255)

   attaGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAggg                                                                                      <  2:5716308/93‑1 (MQ=255)

   attaGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAggg                                                                                      >  1:162319/1‑93 (MQ=255)

    ttaGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAgggg                                                                                     <  2:6066393/93‑1 (MQ=255)

     taGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGa                                                                                    <  1:657974/93‑1 (MQ=255)

      aGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGaa                                                                                   <  1:1479350/93‑1 (MQ=255)

       gCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAAt                                                                                  <  1:377222/93‑1 (MQ=255)

        cAGAAAGATTGCGGTGGTTACAGGAGGACAATACCGGGGAAATAGACGAAGAAGAACTGGCATCCTTTTTATATGCCATTGCAAAGGGGAATg                                                                                 >  2:22570/1‑93 (MQ=255)

        cAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATg                                                                                 <  1:556455/93‑1 (MQ=255)

          gAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGtt                                                                               >  2:48552/1‑93 (MQ=255)

          gAAAGATTGCGGTGGTTACAGGAGGAAAATACAGGGGAAATAGACGAATAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGTGGAATGtt                                                                               >  1:799341/1‑93 (MQ=255)

           aaaGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGttt                                                                              <  1:484455/93‑1 (MQ=255)

            aaGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGtttt                                                                             <  2:121068/93‑1 (MQ=255)

             aGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGttttt                                                                            <  2:1401803/93‑1 (MQ=255)

             aGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGttttt                                                                            >  1:323374/1‑93 (MQ=255)

             aGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGttttt                                                                            >  1:4561851/1‑93 (MQ=255)

             aGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGttttt                                                                            >  1:1169652/1‑93 (MQ=255)

              gATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTa                                                                           <  1:1949233/93‑1 (MQ=255)

               aTTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTaa                                                                          <  2:302123/93‑1 (MQ=255)

               aTTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTaa                                                                          >  1:695406/1‑93 (MQ=255)

               aTTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTaa                                                                          <  2:918799/93‑1 (MQ=255)

                ttGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAAc                                                                         >  2:2280519/1‑93 (MQ=255)

                    ggtggtTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTcc                                                                     <  1:5711297/93‑1 (MQ=255)

                     gtggtTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTTTATGCCATTGCAAAGGGGAAGGTTTTTAATTTCCa                                                                    >  1:3528352/1‑93 (MQ=255)

                        gtTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTACAGAc                                                                 >  1:3259962/1‑93 (MQ=255)

                         ttACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGAcc                                                                <  1:4403531/93‑1 (MQ=255)

                          tACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCa                                                               >  2:3365491/1‑93 (MQ=255)

                          tACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCa                                                               >  1:5887312/1‑93 (MQ=255)

                           aCAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCAt                                                              >  2:4840850/1‑93 (MQ=255)

                           aCAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTCCCAGACCAt                                                              >  1:3196223/1‑93 (MQ=255)

                           aCAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACATCCAGACCAt                                                              <  2:5482643/93‑1 (MQ=255)

                            cAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCAtt                                                             <  2:3671576/93‑1 (MQ=255)

                            cAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCAtt                                                             <  1:2318068/93‑1 (MQ=255)

                            cAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCAtt                                                             <  2:4678887/93‑1 (MQ=255)

                            aaggaggAAAATACCGGGGAAAGAGACGAAGAAGAACTGGACTCCTTTTGATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCAtt                                                             <  1:1883917/92‑1 (MQ=255)

                             aggaggAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTc                                                            >  2:1723879/1‑93 (MQ=255)

                               gaggaAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTg                                                          <  1:6007527/93‑1 (MQ=255)

                                aggaAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGc                                                         <  1:5860049/93‑1 (MQ=255)

                                 ggaAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCa                                                        <  1:4746255/93‑1 (MQ=255)

                                 ggaAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCa                                                        >  1:4690680/1‑93 (MQ=255)

                                  gaAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCAt                                                       <  2:5065079/93‑1 (MQ=255)

                                   aaaaTACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATc                                                      <  2:2453774/93‑1 (MQ=255)

                                   aaaaTACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATc                                                      <  2:1748604/93‑1 (MQ=255)

                                        aCCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGccc                                                 >  2:5249132/1‑93 (MQ=255)

                                        aCCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGccc                                                 >  2:5405498/1‑93 (MQ=255)

                                         ccGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCg                                                <  2:2946068/93‑1 (MQ=255)

                                         ccGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCg                                                <  2:5895667/93‑1 (MQ=255)

                                          cGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGt                                               >  1:2225069/1‑93 (MQ=255)

                                           ggggAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGtt                                              >  1:2184084/1‑93 (MQ=255)

                                           ggggAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGtt                                              <  1:5190021/93‑1 (MQ=255)

                                               aaaTAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCg                                          >  2:2722812/1‑93 (MQ=255)

                                               aaaTAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCg                                          >  2:4185066/1‑93 (MQ=255)

                                                aaTAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGt                                         >  1:2821836/1‑93 (MQ=255)

                                                aaTAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGt                                         >  1:3480544/1‑93 (MQ=255)

                                                      cGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCaaaa                                   >  1:3115291/1‑93 (MQ=255)

                                                      cGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCaaaa                                   <  1:3016160/93‑1 (MQ=255)

                                                      cGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCaaaa                                   >  2:2298369/1‑93 (MQ=255)

                                                       gaagaagaaCTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAAt                                  <  1:4424956/93‑1 (MQ=255)

                                                        aagaagaaCTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATg                                 >  2:1315095/1‑93 (MQ=255)

                                                          gaagaaCTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGac                               >  1:3691611/1‑93 (MQ=255)

                                                             gaaCTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACAcc                            >  2:222716/1‑93 (MQ=255)

                                                               aCTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCAt                          <  1:5144111/93‑1 (MQ=255)

                                                                cTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATc                         <  1:1975827/93‑1 (MQ=255)

                                                                 tGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCg                        >  1:1448163/1‑93 (MQ=255)

                                                                 tGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCg                        <  2:3921556/93‑1 (MQ=255)

                                                                 tGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCg                        <  1:1340167/93‑1 (MQ=255)

                                                                   gAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGAt                      >  2:4616025/1‑93 (MQ=255)

                                                                   gAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGAt                      >  1:6097570/1‑93 (MQ=255)

                                                                   gAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGAt                      >  1:4265387/1‑93 (MQ=255)

                                                                     aTCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGTCACCATCGAttt                    >  1:5477142/1‑93 (MQ=255)

                                                                     aTCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGAttt                    >  1:3860647/1‑93 (MQ=255)

                                                                     aTCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGAttt                    >  1:2958181/1‑93 (MQ=255)

                                                                      tCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGAtttt                   >  2:4388068/1‑93 (MQ=255)

                                                                      tCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGAtttt                   <  2:1662447/93‑1 (MQ=255)

                                                                        cTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTa                 >  2:6002593/1‑93 (MQ=255)

                                                                         tttttATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTAc                >  2:2780519/1‑93 (MQ=255)

                                                                         tgtttATATGCCATTGCAAAGGGGAATGTGTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTAc                >  1:6033400/3‑93 (MQ=255)

                                                                          ttttATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTAcc               <  1:4234925/93‑1 (MQ=255)

                                                                          ttttATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTAcc               <  2:5546488/93‑1 (MQ=255)

                                                                            ttATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCaa             <  1:386966/93‑1 (MQ=255)

                                                                              atatGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAAt           >  1:6128530/1‑93 (MQ=255)

                                                                              atatGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAAt           <  2:395181/93‑1 (MQ=255)

                                                                              atatGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAAt           >  1:3415621/1‑93 (MQ=255)

                                                                               tatGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATg          >  1:254364/1‑93 (MQ=255)

                                                                               tatGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATg          >  1:4718355/1‑93 (MQ=255)

                                                                                atGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGt         >  1:4921525/1‑93 (MQ=255)

                                                                                atGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAAGGt         <  2:2095321/93‑1 (MQ=255)

                                                                                atGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCCGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGt         >  1:1373209/1‑93 (MQ=255)

                                                                                 tGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGtt        <  2:1552811/93‑1 (MQ=255)

                                                                                   ccATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCg      >  1:3002009/1‑93 (MQ=255)

                                                                                   ccATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCg      <  1:5337662/93‑1 (MQ=255)

                                                                                    cATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGc     >  2:2461673/1‑93 (MQ=255)

                                                                                    cATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGc     <  2:2382270/93‑1 (MQ=255)

                                                                                     aTTGCAAAGGGGAATGTTTTTAAATTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCt    <  1:3497336/93‑1 (MQ=255)

                                                                                      ttGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTc   <  1:3984235/93‑1 (MQ=255)

GTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCG  >  NC_000913/677030‑677209

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