BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,782,476	A→G	V56A (GTT→GCT)	yibN ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,782,476	0	A	G	99.7%	1143.9 / ‑5.3	333	V56A (GTT→GCT)	yibN	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (143/189); total (144/189)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 359 total aligned reads displayed.

ACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACA  >  NC_000913/3782386‑3782565

aCATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGcc                                                                                         <  2:6589554/93‑1 (MQ=255)

   ttgttgGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCcaca                                                                                      <  2:8481791/93‑1 (MQ=255)

   ttgttgGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCcaca                                                                                      <  2:6007975/93‑1 (MQ=255)

     gttgGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCcacagccacagc                                                                                    >  2:8109192/1‑93 (MQ=255)

      ttgGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCg                                                                                   <  1:5014394/93‑1 (MQ=255)

       tgGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGt                                                                                  >  1:574992/1‑93 (MQ=255)

        ggCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGtc                                                                                 <  2:807990/93‑1 (MQ=255)

        ggCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGtc                                                                                 <  1:2416850/93‑1 (MQ=255)

           tttGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGtcttc                                                                              <  2:7945554/93‑1 (MQ=255)

            ttGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGtcttct                                                                             >  1:4688157/1‑93 (MQ=255)

              gATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCttt                                                                           <  2:4648743/93‑1 (MQ=255)

              gATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCttt                                                                           <  1:8475102/93‑1 (MQ=255)

               aTTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTttg                                                                          <  1:8504460/93‑1 (MQ=255)

               aTTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTttg                                                                          <  2:3186194/93‑1 (MQ=255)

                tttCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTttgt                                                                         >  2:2780356/1‑93 (MQ=255)

                tttCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTttgt                                                                         <  2:2248485/93‑1 (MQ=255)

                  tCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTttgttg                                                                       <  2:1016745/93‑1 (MQ=255)

                   cGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGa                                                                      >  1:2998150/1‑93 (MQ=255)

                   cGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGa                                                                      <  2:7486982/93‑1 (MQ=255)

                    gCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGAt                                                                     <  2:2457968/93‑1 (MQ=255)

                     cTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATc                                                                    <  2:264445/93‑1 (MQ=255)

                     cTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATc                                                                    >  2:7740711/1‑93 (MQ=255)

                      tCGGCAACAGATTAATAGAACCTGCGGTATGGCCTTTACGGAAGTCGTCACGCTGACGTAAATCCACAGCCACAGCCGCTTCTTTGTTGagcg                                                                   >  1:6043345/1‑90 (MQ=255)

                        ggCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCGTCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGa                                                                 <  2:2972383/93‑1 (MQ=255)

                        ggCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGa                                                                 <  2:744040/93‑1 (MQ=255)

                        ggCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGa                                                                 <  2:6086587/93‑1 (MQ=255)

                        ggCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGa                                                                 >  2:193596/1‑93 (MQ=255)

                        ggCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCANNGCGTCTTCTTTGTTGATANGa                                                                 >  1:1471308/1‑93 (MQ=255)

                         gCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGAc                                                                <  1:1827012/93‑1 (MQ=255)

                         gCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCAAAGCGTCTTCTTTGTTGATCAGAc                                                                >  2:4289547/1‑93 (MQ=255)

                            aCAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGAcgct                                                             >  1:874188/1‑92 (MQ=255)

                            aCAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGAcgcg                                                             >  1:211009/1‑93 (MQ=255)

                             cAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGt                                                            <  1:3430487/93‑1 (MQ=255)

                             cAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGt                                                            >  2:6735463/1‑93 (MQ=255)

                              aGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTa                                                           <  1:6863716/93‑1 (MQ=255)

                              aGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTa                                                           <  2:2832841/93‑1 (MQ=255)

                              aGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTa                                                           <  1:8146492/93‑1 (MQ=255)

                              aGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTa                                                           <  2:2289417/93‑1 (MQ=255)

                               gATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAg                                                          >  1:1004876/1‑93 (MQ=255)

                                aTTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGc                                                         <  2:5258870/93‑1 (MQ=255)

                                aTTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATAACGCTGACCTAAATCCACGGCCACAGCGTCTCCTTTGTTGAGCAGACGGGtcgt                                                         >  2:6407438/1‑90 (MQ=255)

                                 ttAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCt                                                        <  2:7715675/93‑1 (MQ=255)

                                 ttAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCt                                                        >  1:2007258/1‑93 (MQ=255)

                                   aaTAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTc                                                      <  2:5890431/93‑1 (MQ=255)

                                    aTAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTca                                                     >  2:8420287/1‑93 (MQ=255)

                                       gAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTcacca                                                  <  1:230006/93‑1 (MQ=255)

                                        aaCCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTcaccac                                                 <  1:663965/93‑1 (MQ=255)

                                         aCCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACg                                                >  1:7179218/1‑93 (MQ=255)

                                          ccTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGa                                               >  2:1199048/1‑93 (MQ=255)

                                          ccTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGa                                               >  2:3744304/1‑93 (MQ=255)

                                           cTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAg                                              >  1:88507/1‑93 (MQ=255)

                                           cTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAg                                              >  2:414722/1‑93 (MQ=255)

                                           cTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAg                                              >  1:2971176/1‑93 (MQ=255)

                                             gCGATATGGCCTTTACGGAAGTCATCACCCTGACGTAAATCCACAGCCACAGCGTCTTTTTTGTTGTTCAGACGCGAAGCTTCACCACGAGTa                                            >  1:1891340/1‑93 (MQ=255)

                                               gATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAAt                                          >  2:1379531/1‑93 (MQ=255)

                                                atatGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATc                                         <  2:1726191/93‑1 (MQ=255)

                                                    ggCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCt                                     >  2:8633250/1‑93 (MQ=255)

                                                       cTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCa                                  >  1:1730673/1‑93 (MQ=255)

                                                        tttACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCAc                                 <  2:3538760/93‑1 (MQ=255)

                                                        tttACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCAc                                 >  2:2115616/1‑93 (MQ=255)

                                                        tttACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCAc                                 <  1:6283678/93‑1 (MQ=255)

                                                        tttACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCAc                                 <  1:3125641/93‑1 (MQ=255)

                                                          tACGGAAGTCATCACGCTGACCTAAACCCACAGCCACAGCGTCTTCTTTTTTGATCAGACGCGCAGCTTCACCACGAGTAATCACCTTCTCtt                               <  1:8248989/93‑1 (MQ=255)

                                                           aCGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACttt                              <  1:3655759/93‑1 (MQ=255)

                                                            cGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTc                             <  1:1157558/93‑1 (MQ=255)

                                                             ggAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCg                            <  2:1564325/93‑1 (MQ=255)

                                                              gAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGa                           <  2:8601259/93‑1 (MQ=255)

                                                               aaGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAg                          <  1:6179849/93‑1 (MQ=255)

                                                                 gTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGt                        >  2:51219/1‑93 (MQ=255)

                                                                   catcaCGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCa                      <  1:4125416/93‑1 (MQ=255)

                                                                   catcaCGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCa                      <  1:4240963/93‑1 (MQ=255)

                                                                   catcaCGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCa                      >  1:7852494/1‑93 (MQ=255)

                                                                    atcaCGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAg                     <  2:5156893/93‑1 (MQ=255)

                                                                     tcaCGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAgg                    >  1:4821321/1‑93 (MQ=255)

                                                                      caCGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGc                   <  2:3937856/93‑1 (MQ=255)

                                                                      caCGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGc                   >  2:4949649/1‑93 (MQ=255)

                                                                      caCGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGc                   <  1:7798852/93‑1 (MQ=255)

                                                                       aCGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGct                  >  1:2672126/1‑93 (MQ=255)

                                                                       aCGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGct                  >  2:7310750/1‑93 (MQ=255)

                                                                        cGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGctc                 >  2:8765676/1‑93 (MQ=255)

                                                                        cGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGctc                 >  1:743488/1‑93 (MQ=255)

                                                                         gCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGctct                <  1:5537768/93‑1 (MQ=255)

                                                                         gCTGACGTAAACCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGctct                <  2:5459031/93‑1 (MQ=255)

                                                                          cTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCtt               <  1:7601757/93‑1 (MQ=255)

                                                                           tGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTa              <  2:5627690/93‑1 (MQ=255)

                                                                           tGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTCTCGAGGTCAGGCTCTTa              >  2:7811971/1‑93 (MQ=255)

                                                                            gACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTaa             >  2:1286527/1‑93 (MQ=255)

                                                                              cGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAAc           <  1:3220002/93‑1 (MQ=255)

                                                                               gTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACg          <  1:5055760/93‑1 (MQ=255)

                                                                               gTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACTCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACg          <  1:8094669/93‑1 (MQ=255)

                                                                                tAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGt         >  2:8238784/1‑93 (MQ=255)

                                                                                 aaaTCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTa        <  1:4039952/93‑1 (MQ=255)

                                                                                 aaaTCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTa        >  2:4382939/1‑93 (MQ=255)

                                                                                  aaTCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAg       >  1:207438/1‑93 (MQ=255)

                                                                                   aTCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGt      >  1:4329790/1‑93 (MQ=255)

                                                                                    tCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTc     >  2:7128404/1‑93 (MQ=255)

                                                                                    tCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTc     >  2:477392/1‑93 (MQ=255)

                                                                                      cacagccacagcGTCTTCTTTGTTGATCAGACTCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTcac   >  1:3759846/1‑93 (MQ=255)

                                                                                       acaGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTcaca  >  1:3479015/1‑93 (MQ=255)

                                                                                       acaGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTcaca  >  2:3294629/1‑93 (MQ=255)

ACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACA  >  NC_000913/3782386‑3782565

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