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%	1025.8 / ‑5.5	275	I70I (ATC→ATT)	yiaB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (105/169); total (106/169)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-01

Only 100 of 277 total aligned reads displayed.

GGGATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAGATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCAAG  >  NC_000913/3725090‑3725269

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

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

   aTAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATg                                                                                      <  2:260809/93‑1 (MQ=255)

   aTAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATg                                                                                      <  1:2252037/93‑1 (MQ=255)

   aTAGATCACCATTTCATAAGTATTGATAGGCGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATg                                                                                      <  2:5635720/93‑1 (MQ=255)

     aGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATggg                                                                                    <  1:877095/93‑1 (MQ=255)

     aGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATggg                                                                                    <  1:444663/93‑1 (MQ=255)

     aGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATggg                                                                                    <  2:5378158/93‑1 (MQ=255)

     aGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATggg                                                                                    >  1:2521680/1‑93 (MQ=255)

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

        tCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTAGCAAATATGGGTaa                                                                                 >  2:3037464/1‑93 (MQ=255)

         cACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaa                                                                                <  2:2482802/93‑1 (MQ=255)

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

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

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

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

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

             aTTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaaaaa                                                                            <  2:1528703/93‑1 (MQ=255)

             aTTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaaaaa                                                                            <  2:4450882/93‑1 (MQ=255)

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

                tCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAatca                                                                         <  1:4705377/93‑1 (MQ=255)

                  aTAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAatcatc                                                                       <  2:4849659/93‑1 (MQ=255)

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

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

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

                    aaGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCa                                                                     >  1:1131182/1‑93 (MQ=255)

                     aGTATTGATAGGAGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAg                                                                    >  2:3592802/1‑93 (MQ=255)

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

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

                            aTAGGTGCGGTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCg                                                             >  1:2568035/1‑93 (MQ=255)

                             tAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGa                                                            >  2:6487817/1‑93 (MQ=255)

                              aGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAg                                                           >  2:94773/1‑93 (MQ=255)

                              aGGTGAGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAg                                                           >  1:4219775/1‑93 (MQ=255)

                               ggggCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGt                                                          >  2:1688328/4‑93 (MQ=255)

                               ggTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGt                                                          >  2:6529325/1‑93 (MQ=255)

                               ggTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGt                                                          <  2:4414794/93‑1 (MQ=255)

                               ggTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGt                                                          >  2:287857/1‑93 (MQ=255)

                                gTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGtt                                                         <  1:1710244/93‑1 (MQ=255)

                                gTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGtt                                                         <  2:3784717/93‑1 (MQ=255)

                                   cGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGttttt                                                      <  2:146353/93‑1 (MQ=255)

                                   cGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGttttt                                                      <  1:6447554/93‑1 (MQ=255)

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

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

                                    gTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTc                                                     >  2:1028109/1‑93 (MQ=255)

                                      ttAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTg                                                   <  1:6525147/93‑1 (MQ=255)

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

                                        aaaaCACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTgcg                                                 >  1:3070708/1‑93 (MQ=255)

                                          aaCACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGcgcg                                               <  2:4183933/93‑1 (MQ=255)

                                          aaCACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGcgcg                                               <  2:2853606/93‑1 (MQ=255)

                                           acacCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGa                                              >  1:1806772/1‑93 (MQ=255)

                                             acCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAg                                            >  2:202229/1‑93 (MQ=255)

                                              ccGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGa                                           >  2:323829/1‑93 (MQ=255)

                                                gATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGata                                         >  2:4476087/1‑93 (MQ=255)

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

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

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

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

                                                 aTAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGatat                                        >  1:703524/1‑93 (MQ=255)

                                                  tAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATg                                       <  1:2607223/93‑1 (MQ=255)

                                                   aaaCAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGc                                      <  1:557289/93‑1 (MQ=255)

                                                   aaaCAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGc                                      <  2:3639644/93‑1 (MQ=255)

                                                     aCAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCAt                                    <  2:4556672/93‑1 (MQ=255)

                                                       aaGAGACCGATAGTGATTAACGCAACTAACTGGCAGATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAg                                  >  2:4069434/1‑93 (MQ=255)

                                                        agagACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGt                                 >  2:819988/1‑93 (MQ=255)

                                                        agagACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGt                                 >  1:4976425/1‑93 (MQ=255)

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

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

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

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

                                                                aTAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCaaaa                         >  2:6184722/1‑93 (MQ=255)

                                                                 tAGTGATTAACGCAACTAACTGGCAAATATGGGGAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAg                        >  2:2161525/1‑93 (MQ=255)

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

                                                                  aGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGt                       >  1:5916446/1‑93 (MQ=255)

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

                                                                   gTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGtt                      <  1:2981235/93‑1 (MQ=255)

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

                                                                   gTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGtt                      <  2:5743895/93‑1 (MQ=255)

                                                                    tGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTg                     >  1:5479216/1‑93 (MQ=255)

                                                                    tGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTg                     >  1:5692050/1‑93 (MQ=255)

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

                                                                      aTTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCt                   >  1:95820/1‑93 (MQ=255)

                                                                        tAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGt                 <  1:4789422/93‑1 (MQ=255)

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

                                                                           cGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCAt              <  1:1129097/93‑1 (MQ=255)

                                                                           cGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCAt              <  1:621744/93‑1 (MQ=255)

                                                                           cGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCAt              <  1:5110498/93‑1 (MQ=255)

                                                                           cGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCAt              <  1:3459319/93‑1 (MQ=255)

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

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

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

                                                                              aactaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTaa           >  2:3995951/1‑93 (MQ=255)

                                                                               actaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAac          >  1:3361071/1‑93 (MQ=255)

                                                                                ctaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAaca         <  1:1798921/93‑1 (MQ=255)

                                                                                ctaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAaca         <  1:5608957/93‑1 (MQ=255)

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

                                                                                   actGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGc      <  1:1466082/93‑1 (MQ=255)

                                                                                   actGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGc      <  2:6572853/93‑1 (MQ=255)

                                                                                    ctGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGcc     >  1:5388561/1‑93 (MQ=255)

                                                                                     tGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACGCGCCa    <  1:2063341/93‑1 (MQ=255)

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

GGGATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAGATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCAAG  >  NC_000913/3725090‑3725269

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