BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,520,082	A→G	H527H (CAT→CAC)	yjhG ←	putative dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,520,082	0	A	G	99.7%	1129.8 / ‑5.4	356	H527H (CAT→CAC)	yjhG	putative dehydratase
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (161/194); total (162/194)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

Only 100 of 362 total aligned reads displayed.

GGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACATGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGAC  >  NC_000913/4519992‑4520171

ggggCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTg                                                                                         >  1:390774/1‑93 (MQ=255)

 gggCCTCCGCCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGc                                                                                        <  2:5544962/93‑1 (MQ=255)

    ccTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTa                                                                                     >  2:1255021/1‑93 (MQ=255)

    ccTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTa                                                                                     >  2:3444927/1‑93 (MQ=255)

     cTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTAc                                                                                    >  1:5006022/1‑93 (MQ=255)

     cTCCGGACAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTAc                                                                                    <  2:2990953/93‑1 (MQ=255)

      tCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTAcc                                                                                   >  2:109016/1‑93 (MQ=255)

      tCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTAcc                                                                                   >  2:2238511/1‑93 (MQ=255)

       ccGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCg                                                                                  >  1:2477637/1‑93 (MQ=255)

        cGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGt                                                                                 >  1:1392780/1‑93 (MQ=255)

         ggCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTa                                                                                >  1:6512161/1‑93 (MQ=255)

           ccAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATg                                                                              >  1:3737862/1‑93 (MQ=255)

            cAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGa                                                                             <  2:3113831/93‑1 (MQ=255)

              gcgcTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACa                                                                           <  2:3831763/93‑1 (MQ=255)

                gcTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGa                                                                         <  1:762284/93‑1 (MQ=255)

                gcTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGa                                                                         <  1:4731617/93‑1 (MQ=255)

                gcTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGa                                                                         <  1:7234481/93‑1 (MQ=255)

                gcTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGa                                                                         <  2:6631337/93‑1 (MQ=255)

                gcTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGa                                                                         <  1:2769130/93‑1 (MQ=255)

                 cTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGAt                                                                        >  1:4594170/1‑93 (MQ=255)

                  ttCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATg                                                                       >  2:262451/1‑93 (MQ=255)

                    cTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCt                                                                     <  2:7931821/93‑1 (MQ=255)

                     tGGCCCCAAATGGCCGATGCAAGCGCCAGTAGAAACGCCCGAGAAACGGGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCtt                                                                    <  2:5342256/93‑1 (MQ=255)

                      ggCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTc                                                                   <  1:1813412/93‑1 (MQ=255)

                       gCCCCACATGTCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCa                                                                  <  1:6334121/93‑1 (MQ=255)

                       gCCCCACATGGCCGATGCACGCGCGAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCa                                                                  <  2:7145025/93‑1 (MQ=255)

                       gCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCa                                                                  >  2:18094/1‑93 (MQ=255)

                       gCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCa                                                                  <  2:6812043/93‑1 (MQ=255)

                       gCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCa                                                                  <  1:5280064/93‑1 (MQ=255)

                        ccccACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAg                                                                 <  2:1437855/93‑1 (MQ=255)

                         cccACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAgg                                                                <  2:2669812/93‑1 (MQ=255)

                         cccACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAgg                                                                >  2:6694131/1‑93 (MQ=255)

                         ccaacaTGGCCGATGCGCGCGCCAGTAGAACCGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAgg                                                                <  1:1977933/90‑1 (MQ=255)

                          ccACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAggg                                                               >  1:2291617/1‑93 (MQ=255)

                           cacaTGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGc                                                              >  1:1666752/1‑93 (MQ=255)

                             caTGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGCCAGATGCTTCAGGGCAc                                                            >  2:6114971/1‑93 (MQ=255)

                             caTGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCAc                                                            <  1:1360042/93‑1 (MQ=255)

                             caTGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCAc                                                            >  1:1925826/1‑93 (MQ=255)

                                ggCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTgg                                                         <  1:3149606/93‑1 (MQ=255)

                                ggCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTgg                                                         <  1:6743196/93‑1 (MQ=255)

                                    gATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAAc                                                     >  2:3795299/1‑93 (MQ=255)

                                    gATGCACGCGCCAGTAGAAACGACCGAGAAACGTGCATCGGTGATTAACGAAACATGCTTATCGTATGACAGATGCTTCAGGGCACTGGTAAc                                                     >  2:4686002/1‑93 (MQ=255)

                                     aTGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAAcc                                                    <  1:3492095/93‑1 (MQ=255)

                                     aTGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAAcc                                                    <  2:1563732/93‑1 (MQ=255)

                                         aCGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGt                                                >  2:723415/1‑93 (MQ=255)

                                          cgcgCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCa                                                                  <  2:5900038/74‑1 (MQ=255)

                                          cgcgCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTa                                               <  2:3675580/93‑1 (MQ=255)

                                          cgcgCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTa                                               <  2:44936/93‑1 (MQ=255)

                                           gcgcCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAGCGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAg                                              <  1:2915477/93‑1 (MQ=255)

                                              ccAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGtt                                           >  1:7290863/1‑93 (MQ=255)

                                              ccAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGtt                                           >  2:4367353/1‑93 (MQ=255)

                                                aGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGCCAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttc                                         >  1:961773/1‑93 (MQ=255)

                                                aGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttc                                         >  1:1043339/1‑93 (MQ=255)

                                                  tAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACATGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttctt                                       >  1:7540983/1‑93 (MQ=255)

                                                   aGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttcttc                                      <  1:4131429/93‑1 (MQ=255)

                                                   aGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttcttc                                      <  2:3497338/93‑1 (MQ=255)

                                                    gAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTcc                                     <  2:4195221/93‑1 (MQ=255)

                                                    gAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTcc                                     >  2:5148495/1‑93 (MQ=255)

                                                    gAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTcc                                     >  2:2732108/1‑93 (MQ=255)

                                                      aaCGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCAt                                   >  2:3281757/1‑93 (MQ=255)

                                                            cGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGt                             >  1:6184961/1‑93 (MQ=255)

                                                            cGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGt                             <  1:2996575/93‑1 (MQ=255)

                                                            cGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGt                             >  1:3836309/1‑93 (MQ=255)

                                                            cGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGt                             <  2:6887107/93‑1 (MQ=255)

                                                              agaAAGGTGAAT‑GGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTAcc                          <  2:7055583/93‑1 (MQ=255)

                                                              agaAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTAc                           <  1:7846631/93‑1 (MQ=255)

                                                              agaAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTAc                           <  2:5472831/93‑1 (MQ=255)

                                                               gaAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTAcc                          >  1:4326055/1‑93 (MQ=255)

                                                               gaAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTAcc                          <  1:5390277/93‑1 (MQ=255)

                                                               gaAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTAcc                          >  1:5119817/1‑93 (MQ=255)

                                                                aaaCGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCt                         <  2:2817802/93‑1 (MQ=255)

                                                                 aaCGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTg                        >  1:6611242/1‑93 (MQ=255)

                                                                  aCGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGa                       <  1:2604480/93‑1 (MQ=255)

                                                                  aCGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGa                       >  2:7583330/1‑93 (MQ=255)

                                                                  aCGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGa                       >  1:5622553/1‑93 (MQ=255)

                                                                   cGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGaa                      <  1:3634253/93‑1 (MQ=255)

                                                                    gTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAg                     <  1:3672960/93‑1 (MQ=255)

                                                                    gTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAg                     <  2:5005419/93‑1 (MQ=255)

                                                                     tGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAgg                    >  1:2075222/1‑93 (MQ=255)

                                                                        aTCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTcc                 <  1:1251493/93‑1 (MQ=255)

                                                                        aTCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTcc                 >  1:6930250/1‑93 (MQ=255)

                                                                        aTCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTcc                 >  1:6223269/1‑93 (MQ=255)

                                                                         tCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCa                >  2:4810844/1‑93 (MQ=255)

                                                                         tCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCa                >  2:7481173/1‑93 (MQ=255)

                                                                         tCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCa                >  2:2517153/1‑93 (MQ=255)

                                                                          ctgTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCaa               <  2:5791456/91‑1 (MQ=255)

                                                                          cGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCaa               <  1:3435227/93‑1 (MQ=255)

                                                                          cGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCaa               <  2:2376252/93‑1 (MQ=255)

                                                                            gTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAAcc             >  1:908129/1‑92 (MQ=255)

                                                                                ttAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaa         >  1:5444471/1‑93 (MQ=255)

                                                                                ttAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaa         >  2:6478319/1‑93 (MQ=255)

                                                                                ttAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaa         <  1:7366476/93‑1 (MQ=255)

                                                                                ttAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaa         <  1:5755872/93‑1 (MQ=255)

                                                                                   aCGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaataa      <  1:2693976/93‑1 (MQ=255)

                                                                                   aCGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaataa      >  2:4078403/1‑93 (MQ=255)

                                                                                     gAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCCTCCCTGTACCTGAAGGTACAACGCCAATAATg    >  2:2082074/1‑93 (MQ=255)

                                                                                     gAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATACTg    >  2:4602230/1‑93 (MQ=255)

                                                                                      aaaCGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGa   >  1:4068682/1‑93 (MQ=255)

                                                                                       aaCGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGAc  <  1:6876305/93‑1 (MQ=255)

                                                                                       aaCGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGAc  <  2:3906691/93‑1 (MQ=255)

GGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACATGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGAC  >  NC_000913/4519992‑4520171

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