BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,725,180	G→A	I70I (ATC→ATT)	yiaB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,725,180	0	G	A	99.6%	930.4 / ‑5.3	264	I70I (ATC→ATT)	yiaB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (81/182); total (81/183)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

Only 100 of 266 total aligned reads displayed.

GGGATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAGATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCAAG  >  NC_000913/3725090‑3725269

gggATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAat                                                                                         <  2:2219851/93‑1 (MQ=255)

 ggATAGATCTCCATTTCATAAGTATTGATAGGTGCGTTTAACACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAata                                                                                        <  2:4223008/93‑1 (MQ=255)

 ggATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAata                                                                                        <  1:973975/93‑1 (MQ=255)

 ggATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAata                                                                                        <  2:774174/93‑1 (MQ=255)

  gATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAatat                                                                                       <  2:2927946/93‑1 (MQ=255)

    tAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATgg                                                                                     <  1:1998101/93‑1 (MQ=255)

    tAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATgg                                                                                     <  2:3004202/93‑1 (MQ=255)

     aGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCATCTAACTGGCAAATATggg                                                                                    <  2:4645561/93‑1 (MQ=255)

      gATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGt                                                                                   >  1:1542789/1‑93 (MQ=255)

       aTCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTa                                                                                  >  1:365100/1‑93 (MQ=255)

       aTCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTa                                                                                  >  2:1321118/1‑93 (MQ=255)

       aTCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTa                                                                                  <  2:1432163/93‑1 (MQ=255)

        tCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaa                                                                                 <  2:9042/93‑1 (MQ=255)

         cACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaa                                                                                >  1:1656325/1‑93 (MQ=255)

          aCCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaa                                                                               >  1:4297293/1‑93 (MQ=255)

           ccATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaaa                                                                              >  1:1033122/1‑93 (MQ=255)

            cATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaaaa                                                                             <  2:3923626/93‑1 (MQ=255)

            cATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaaaa                                                                             >  1:4032034/1‑93 (MQ=255)

             aTTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaaaaa                                                                            >  1:4722617/1‑93 (MQ=255)

               ttCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAatc                                                                          <  1:4071444/93‑1 (MQ=255)

                tCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAatca                                                                         >  2:387976/1‑93 (MQ=255)

                 cATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAatcat                                                                        >  2:2527070/1‑93 (MQ=255)

                  aTAAGTATTGATAGGTGCGTTTACAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAatcatc                                                                       <  2:3223427/93‑1 (MQ=255)

                  aTAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAatcatc                                                                       >  1:2883970/1‑93 (MQ=255)

                   tAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATcc                                                                      >  2:117370/1‑93 (MQ=255)

                   tAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATcc                                                                      >  2:941492/1‑93 (MQ=255)

                   tAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATcc                                                                      <  2:4421718/93‑1 (MQ=255)

                    aaGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCa                                                                     <  1:4579916/93‑1 (MQ=255)

                    aaGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCa                                                                     <  1:105814/93‑1 (MQ=255)

                      gTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGt                                                                   <  1:3467036/93‑1 (MQ=255)

                      gTATTGATAGGTGCGTTGAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGt                                                                   <  2:4118269/93‑1 (MQ=255)

                       tATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGtt                                                                  <  1:2902843/93‑1 (MQ=255)

                       tATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGtt                                                                  <  2:1871421/93‑1 (MQ=255)

                       tATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGtt                                                                  <  1:3275452/93‑1 (MQ=255)

                         ttGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGGAAAAAAATCATCCAGTTGc                                                                >  1:234480/1‑93 (MQ=255)

                         ttGATAGGTGCATTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGc                                                                <  2:4356991/93‑1 (MQ=255)

                           gATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGccc                                                              <  1:3771106/93‑1 (MQ=255)

                           gATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGccc                                                              >  1:2473624/1‑93 (MQ=255)

                             tAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGa                                                            <  2:1225343/93‑1 (MQ=255)

                             tAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGa                                                            <  1:1267747/93‑1 (MQ=255)

                             tAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGa                                                            <  1:770043/93‑1 (MQ=255)

                              aGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAg                                                           <  2:2319383/93‑1 (MQ=255)

                                gTGCTTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGtt                                                         <  1:2031372/93‑1 (MQ=255)

                                    gTTTAAAACACCGATAAGCAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTc                                                     >  2:2474664/1‑93 (MQ=255)

                                    gTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTc                                                     <  1:888879/93‑1 (MQ=255)

                                     tttAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCt                                                    <  1:599743/93‑1 (MQ=255)

                                     tttAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCt                                                    >  1:4108382/1‑93 (MQ=255)

                                       tAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTgc                                                  <  2:3079671/93‑1 (MQ=255)

                                        aaaaCACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTgcg                                                 <  1:658866/93‑1 (MQ=255)

                                         aaaCACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTgcgc                                                <  1:72958/93‑1 (MQ=255)

                                         aaaCACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTgcgc                                                >  1:86777/1‑93 (MQ=255)

                                         aaaCACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTgcgc                                                >  2:480866/1‑93 (MQ=255)

                                         aaaCACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTgcgc                                                <  1:3595910/93‑1 (MQ=255)

                                          aaCACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGcgcg                                               >  1:2177860/1‑93 (MQ=255)

                                            cacCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGaa                                             <  2:845886/93‑1 (MQ=255)

                                            cacCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGaa                                             <  1:81287/93‑1 (MQ=255)

                                              ccGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGa                                           <  1:3448480/93‑1 (MQ=255)

                                              ccGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGa                                           <  2:3251187/93‑1 (MQ=255)

                                               cGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGat                                          <  2:3567278/93‑1 (MQ=255)

                                                gATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGata                                         <  2:4501318/93‑1 (MQ=255)

                                                gATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGata                                         <  1:2075974/93‑1 (MQ=255)

                                                gATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGata                                         <  2:1656480/93‑1 (MQ=255)

                                                 aTAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGatat                                        <  1:946228/93‑1 (MQ=255)

                                                   aaaCAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGc                                      >  2:4320855/1‑93 (MQ=255)

                                                    aaCAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCa                                     >  2:3227731/1‑93 (MQ=255)

                                                       aaGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAg                                  >  1:2091886/1‑93 (MQ=255)

                                                        agagACCGATATTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGt                                 <  2:1337684/93‑1 (MQ=255)

                                                        agagACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGt                                 <  2:3813433/93‑1 (MQ=255)

                                                         gagaCCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGtt                                <  1:4432915/93‑1 (MQ=255)

                                                           gaCCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGc                              >  2:1544652/1‑93 (MQ=255)

                                                            aCCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGcc                             <  2:2091289/93‑1 (MQ=255)

                                                             ccGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCa                            >  1:3927624/1‑93 (MQ=255)

                                                              cGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCaa                           <  2:4625859/93‑1 (MQ=255)

                                                               caTAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCaaa                          <  1:26119/92‑1 (MQ=255)

                                                                aTAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCaaaa                         <  2:2677355/93‑1 (MQ=255)

                                                                aTAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCaaaa                         <  2:1803160/93‑1 (MQ=255)

                                                                 tAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAg                        <  1:4396443/93‑1 (MQ=255)

                                                                  aGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGt                       <  1:2831830/93‑1 (MQ=255)

                                                                   gTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGtt                      >  1:3948693/1‑93 (MQ=255)

                                                                     gATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGc                    <  2:1124343/93‑1 (MQ=255)

                                                                      aTTAACGCAACTAACTGGCAAATGTGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCt                   <  2:2580065/93‑1 (MQ=255)

                                                                      aTTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCt                   <  2:3321422/93‑1 (MQ=255)

                                                                      aTTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCt                   <  1:3070488/93‑1 (MQ=255)

                                                                      aTTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCt                   <  2:2541245/93‑1 (MQ=255)

                                                                       ttAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTg                  >  1:2119863/1‑93 (MQ=255)

                                                                         aaCGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTc                <  2:3545027/93‑1 (MQ=255)

                                                                          aCGCAACTACCTGGCAACTATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCGGTCa               <  1:2712825/93‑1 (MQ=255)

                                                                          aCGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCa               <  1:2405103/93‑1 (MQ=255)

                                                                          aCGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCa               <  1:1736480/93‑1 (MQ=255)

                                                                          aCGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCa               <  1:1427552/93‑1 (MQ=255)

                                                                          aCGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTGGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCa               <  2:2789949/93‑1 (MQ=255)

                                                                           cGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCAt              >  1:31770/1‑93 (MQ=255)

                                                                            gCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCAtt             <  1:388515/93‑1 (MQ=255)

                                                                              aactaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTaa           <  1:3138747/93‑1 (MQ=255)

                                                                              aactaactGGCAAATATGGGTAAAAAAATCATACAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTaa           >  2:1944415/1‑93 (MQ=255)

                                                                                  aactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACg       >  2:60061/1‑93 (MQ=255)

                                                                                  aactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACg       >  2:615763/1‑93 (MQ=255)

                                                                                      ggCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCaa   <  1:3723153/93‑1 (MQ=255)

                                                                                      ggCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCaa   <  2:2965028/93‑1 (MQ=255)

                                                                                       gCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCAAg  >  1:2710938/1‑93 (MQ=255)

GGGATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAGATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCAAG  >  NC_000913/3725090‑3725269

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