BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,945,896	T→C	V63A (GTA→GCA)	yifE →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,945,896	0	T	C	99.3%	957.7 / ‑3.7	296	V63A (GTA→GCA)	yifE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (155/139); total (156/140)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 313 total aligned reads displayed.

AACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACA  >  NC_000913/3945806‑3945985

aaCTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCAt                                                                                         >  1:1620734/1‑93 (MQ=255)

aaCTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCAt                                                                                         <  2:1944322/93‑1 (MQ=255)

 aCTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATg                                                                                        >  1:452121/1‑93 (MQ=255)

  cTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGc                                                                                       >  1:50188/1‑93 (MQ=255)

   tGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGcc                                                                                      >  2:5710173/1‑93 (MQ=255)

    gCTTGAGCGTCTTGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCg                                                                                     <  1:2051837/93‑1 (MQ=255)

    gCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCg                                                                                     >  2:3893344/1‑93 (MQ=255)

    gCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCg                                                                                     >  1:3785406/1‑93 (MQ=255)

    gCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGACg                                                                                     >  2:5518554/1‑93 (MQ=255)

     cTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGt                                                                                    >  2:33409/1‑93 (MQ=255)

     cTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGt                                                                                    >  2:270984/1‑93 (MQ=255)

     cTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGt                                                                                    <  1:4734729/93‑1 (MQ=255)

      ttGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTg                                                                                   >  2:4391781/1‑93 (MQ=255)

       tGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTgg                                                                                  >  2:4176400/1‑93 (MQ=255)

       tGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTgg                                                                                  <  1:3955626/93‑1 (MQ=255)

          gCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGa                                                                               >  2:2084379/1‑93 (MQ=255)

           cGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGaa                                                                              <  2:2467954/93‑1 (MQ=255)

             tCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACg                                                                            >  1:1284812/1‑93 (MQ=255)

              cATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGt                                                                           <  2:4494484/93‑1 (MQ=255)

                   ttATGCCTTCAATGAGTTGGATCTTGGTAAACGCGAGCCGGTTACCGAGGAAGAGAAAATCTTCGTAGCAGCATGCCGTGGCGAACGTGAGcc                                                                      >  2:3601768/1‑93 (MQ=255)

                   ttATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGcc                                                                      <  1:1810769/93‑1 (MQ=255)

                     aTGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAg                                                                    >  2:5570134/1‑93 (MQ=255)

                     aTGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAg                                                                    >  1:5980548/1‑93 (MQ=255)

                      tGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGt                                                                   >  1:2899947/1‑93 (MQ=255)

                      tGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGt                                                                   <  1:638527/93‑1 (MQ=255)

                       gCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTg                                                                  >  1:2653989/1‑93 (MQ=255)

                       gCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTg                                                                  >  1:5757116/1‑93 (MQ=255)

                         cTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGAc                                                                >  2:6614699/1‑93 (MQ=255)

                             aaTGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGaa                                                            >  1:609357/1‑93 (MQ=255)

                                gAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCa                                                         <  2:6067255/93‑1 (MQ=255)

                                 aGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAg                                                        <  1:3365948/93‑1 (MQ=255)

                                   ttGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGaa                                                      >  2:4876274/1‑93 (MQ=255)

                                   ttGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGaa                                                      <  1:517829/93‑1 (MQ=255)

                                     ggATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAcg                                                    >  1:864183/1‑93 (MQ=255)

                                     ggATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAcg                                                    >  1:4590858/1‑93 (MQ=255)

                                      gATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGGCAGTGACAGAAGCAGAAcgc                                                   >  1:4827856/1‑93 (MQ=255)

                                       aTCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAcgcg                                                  >  2:5679931/1‑93 (MQ=255)

                                       aTCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAcgcg                                                  >  2:3015761/1‑93 (MQ=255)

                                        tCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCgt                                                 >  1:1940547/1‑93 (MQ=255)

                                           tGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGtgtg                                              <  1:1060752/93‑1 (MQ=255)

                                               aaaCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTcc                                          <  2:5076979/93‑1 (MQ=255)

                                               aaaCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTcc                                          <  2:457966/93‑1 (MQ=255)

                                               aaaCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTcc                                          <  1:3168057/93‑1 (MQ=255)

                                                aaCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCa                                         <  1:5063090/93‑1 (MQ=255)

                                                aaCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCa                                         >  2:378126/1‑93 (MQ=255)

                                                aaCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCa                                         <  1:1481411/93‑1 (MQ=255)

                                                aaCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGGGTGGTCCa                                         >  1:74899/1‑93 (MQ=255)

                                                aaCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTCGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCa                                         >  1:3129156/1‑93 (MQ=255)

                                                 aCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCaa                                        >  2:5453814/1‑93 (MQ=255)

                                                 aCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCaa                                        <  1:2191823/93‑1 (MQ=255)

                                                   gcgAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGt                                      <  1:940728/93‑1 (MQ=255)

                                                   gcgAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGt                                      <  1:1540134/93‑1 (MQ=255)

                                                    cgAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGta                                     <  1:1771386/93‑1 (MQ=255)

                                                    cgAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGta                                     <  1:2580996/93‑1 (MQ=255)

                                                    cgAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGta                                     <  1:4238658/93‑1 (MQ=255)

                                                     gAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGtat                                    >  1:1888089/1‑93 (MQ=255)

                                                      aGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGtata                                   <  1:404439/93‑1 (MQ=255)

                                                      aGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGtata                                   >  1:3406611/1‑93 (MQ=255)

                                                      aGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGtata                                   >  2:2843437/1‑93 (MQ=255)

                                                      aGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGtata                                   >  1:6583284/1‑93 (MQ=255)

                                                         cGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGa                                >  2:786558/1‑93 (MQ=255)

                                                           gTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAcg                              >  2:5867519/1‑93 (MQ=255)

                                                            ttACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAcgc                             <  2:1350208/93‑1 (MQ=255)

                                                             tACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAcgcg                            <  1:6145721/93‑1 (MQ=255)

                                                             tACCGAGAAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAcgcg                            <  2:4115158/93‑1 (MQ=255)

                                                               ccGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTa                          >  1:2731392/1‑93 (MQ=255)

                                                                cGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTAt                         <  2:1739167/93‑1 (MQ=255)

                                                                  aGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTa                       >  2:348598/1‑93 (MQ=255)

                                                                  aGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATCACGCGTATTa                       <  1:4076064/93‑1 (MQ=255)

                                                                  aGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGAGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTa                       <  2:2606151/93‑1 (MQ=255)

                                                                   ggAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTaa                      >  2:659227/1‑93 (MQ=255)

                                                                   ggAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTaa                      >  2:1991763/1‑93 (MQ=255)

                                                                   ggAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTaa                      <  1:566531/93‑1 (MQ=255)

                                                                   ggAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTaa                      >  1:5293269/1‑93 (MQ=255)

                                                                    gAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAg                     <  1:2355069/93‑1 (MQ=255)

                                                                     aaGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGc                    <  1:4345497/93‑1 (MQ=255)

                                                                      agagAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCg                   <  1:3000338/93‑1 (MQ=255)

                                                                      agagAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGAGTGTGGTCCAAGTATATGACGCGTATTAAGCg                   >  2:3698695/1‑93 (MQ=255)

                                                                          aaaCTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCa               >  2:1633513/1‑93 (MQ=255)

                                                                          aaaCTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCa               >  2:225131/1‑93 (MQ=255)

                                                                           aaCTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCaa              <  2:5305817/93‑1 (MQ=255)

                                                                            aCTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCaaa             >  2:1879818/1‑93 (MQ=255)

                                                                            aCTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCaaa             >  1:5120695/1‑93 (MQ=255)

                                                                             ctctTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCaaaa            <  2:842585/93‑1 (MQ=255)

                                                                             ctctTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCaaaa            >  2:6470450/1‑93 (MQ=255)

                                                                             ctctTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCaaaa            <  2:6720681/93‑1 (MQ=255)

                                                                              tctTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCaaaag           >  2:3268496/1‑92 (MQ=255)

                                                                              tctTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAAc           >  1:5470620/1‑93 (MQ=255)

                                                                              tctTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCATCCAAAAc           >  1:4141535/1‑93 (MQ=255)

                                                                                 tCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGtt        <  2:6408115/93‑1 (MQ=255)

                                                                                 tCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGtt        <  1:2509018/93‑1 (MQ=255)

                                                                                 tCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGtt        >  1:3483777/1‑93 (MQ=255)

                                                                                  cGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGttt       <  1:3258661/93‑1 (MQ=255)

                                                                                  cGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGttt       >  2:1506315/1‑93 (MQ=255)

                                                                                  cGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGttt       <  2:1126254/93‑1 (MQ=255)

                                                                                  cGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGttt       <  1:4877995/93‑1 (MQ=255)

                                                                                    tAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGCGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTc     <  2:4688233/93‑1 (MQ=255)

                                                                                      gcagcaTGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTcac   <  2:4443499/93‑1 (MQ=255)

                                                                                      gcagcaTGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTcac   <  2:2799877/93‑1 (MQ=255)

                                                                                      gcagcaTGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTcac   <  2:6147914/93‑1 (MQ=255)

AACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACA  >  NC_000913/3945806‑3945985

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