BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,460,324	A→G	D59G (GAC→GGC)	b1398 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,460,324	0	A	G	99.6%	878.9 / NA	257	D59G (GAC→GGC)	b1398	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (122/134); minor base T (0/1); total (122/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 262 total aligned reads displayed.

CGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGC  >  NC_000913/1460234‑1460413

cGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGcc                                                                                         <  2:2530879/93‑1 (MQ=255)

cGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGcc                                                                                         >  2:1538357/1‑93 (MQ=255)

cGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGcc                                                                                         >  2:3757426/1‑93 (MQ=255)

 gCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCt                                                                                        >  1:278096/1‑93 (MQ=255)

  cTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTg                                                                                       >  2:2343519/1‑93 (MQ=255)

    cAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCg                                                                                     <  2:3913758/93‑1 (MQ=255)

       acacGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTaa                                                                                  >  1:1428242/1‑93 (MQ=255)

       acacGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTaa                                                                                  >  1:1975518/1‑93 (MQ=255)

       acacGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTaa                                                                                  >  2:3820151/1‑93 (MQ=255)

         acGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAAt                                                                                <  1:6111961/93‑1 (MQ=255)

           gcgTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAAttt                                                                              >  2:5909270/1‑93 (MQ=255)

             gTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTcc                                                                            >  2:5982455/1‑93 (MQ=255)

             gTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTcc                                                                            <  2:753226/93‑1 (MQ=255)

              tATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTccc                                                                           <  2:3076167/93‑1 (MQ=255)

               aTGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCt                                                                          >  2:4678169/1‑93 (MQ=255)

                    aaTGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTacc                                                                     >  1:3256240/1‑93 (MQ=255)

                     aTGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTacca                                                                    >  1:2951153/1‑93 (MQ=255)

                       gtgCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTaccacc                                                                  <  2:6256060/93‑1 (MQ=255)

                        tgCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTAccacca                                                                 <  1:698193/93‑1 (MQ=255)

                         gCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCaa                                                                >  2:6283682/1‑93 (MQ=255)

                           cGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAAc                                                              >  2:5146840/1‑93 (MQ=255)

                             aTGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAg                                                            <  2:4249790/93‑1 (MQ=255)

                              tGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAgg                                                           <  2:3532578/93‑1 (MQ=255)

                                 aTCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATc                                                        <  2:5669943/93‑1 (MQ=255)

                                 aTCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATc                                                        >  1:1483786/1‑93 (MQ=255)

                                  tCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCt                                                       >  2:1172236/1‑93 (MQ=255)

                                   cGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTg                                                      <  1:3374276/93‑1 (MQ=255)

                                     gcgcAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTgcg                                                    <  1:4016678/93‑1 (MQ=255)

                                      cgcAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTgcgc                                                   <  2:2647172/93‑1 (MQ=255)

                                       gcAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGcgcg                                                  <  2:2718526/93‑1 (MQ=255)

                                         aaaTTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGAc                                                >  2:2236667/1‑93 (MQ=255)

                                         aaaTTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGAc                                                >  2:3343930/1‑93 (MQ=255)

                                          aaTTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACa                                               >  1:3430207/1‑93 (MQ=255)

                                          aaTTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACa                                               >  2:5431943/1‑93 (MQ=255)

                                          aaTTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACa                                               >  1:564141/1‑93 (MQ=255)

                                            ttCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAAc                                             <  2:4888666/93‑1 (MQ=255)

                                               gACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTAt                                          <  1:5256681/93‑1 (MQ=255)

                                               gACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTAt                                          >  2:1353262/1‑93 (MQ=255)

                                               gACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTAt                                          <  2:4034397/93‑1 (MQ=255)

                                               gACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTAt                                          >  2:3383657/1‑93 (MQ=255)

                                                aCGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATc                                         <  2:6340805/93‑1 (MQ=255)

                                                aCGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATc                                         >  2:2053677/1‑93 (MQ=255)

                                                 cGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATcc                                        <  2:3863231/93‑1 (MQ=255)

                                                 cGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATcc                                        <  1:1727154/93‑1 (MQ=255)

                                                  gcagcaGGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATccc                                       <  1:6233946/93‑1 (MQ=255)

                                                   cagcagGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCt                                      <  1:5383404/93‑1 (MQ=255)

                                                    agcagGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCtt                                     <  1:2412684/93‑1 (MQ=255)

                                                    agcagGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCtt                                     >  2:2520012/1‑93 (MQ=255)

                                                     gcagGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCttt                                    >  1:1009780/1‑93 (MQ=255)

                                                     gcagGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCttt                                    <  2:127617/93‑1 (MQ=255)

                                                     gcagGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCttt                                    <  1:3574763/93‑1 (MQ=255)

                                                     gcagGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCttt                                    <  1:5464804/93‑1 (MQ=255)

                                                     gcagGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCttt                                    >  2:1025223/1‑93 (MQ=255)

                                                      cagGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTg                                   <  1:3849968/93‑1 (MQ=255)

                                                       agGCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGa                                  >  2:5788852/1‑93 (MQ=255)

                                                         gCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGaca                                >  2:5005151/1‑93 (MQ=255)

                                                         gCGTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGaca                                >  1:2393805/1‑93 (MQ=255)

                                                           gTACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACAcc                              >  1:3156029/1‑93 (MQ=255)

                                                            tACATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCt                             >  1:5838602/1‑93 (MQ=255)

                                                             aCATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCtt                            <  1:145835/93‑1 (MQ=255)

                                                              cATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTCCCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCttt                           <  1:592496/93‑1 (MQ=255)

                                                              cATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCttt                           <  1:3750540/93‑1 (MQ=255)

                                                              cATCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCttt                           <  1:5594136/93‑1 (MQ=255)

                                                              cATCCTGATGACTTCAGGGAACTTTCAGGCATGCGTAATTTTCCCTGTACCACCAAACTGGATCTGCGCAACAACTATCCCTTTGACACCttt                           <  1:495000/93‑1 (MQ=255)

                                                               aTCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTg                          <  2:2157550/93‑1 (MQ=255)

                                                               aTCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTg                          >  1:3591188/1‑93 (MQ=255)

                                                                tCCTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGc                         <  2:4606700/93‑1 (MQ=255)

                                                                 ccTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGcc                        <  1:1283825/93‑1 (MQ=255)

                                                                 ccTGATGATTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGcc                        >  2:192182/1‑93 (MQ=255)

                                                                   tgatgaTTTCAGGGAACTTTCAGTCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGt                      <  2:6025885/93‑1 (MQ=255)

                                                                   tgatgaTTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGt                      >  1:4839656/1‑93 (MQ=255)

                                                                   tgatgaTTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGt                      <  2:4831537/93‑1 (MQ=255)

                                                                   tgatgaTTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGt                      <  2:6112220/93‑1 (MQ=255)

                                                                      tgatTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGcc                   <  2:2749178/93‑1 (MQ=255)

                                                                      tgatTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGcc                   <  1:3697146/93‑1 (MQ=255)

                                                                      tgatTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGcc                   <  1:788358/93‑1 (MQ=255)

                                                                       gatTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCa                  >  2:953646/1‑93 (MQ=255)

                                                                       gatTTCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCa                  >  2:1703266/1‑93 (MQ=255)

                                                                         tttCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAAt                <  1:5029283/93‑1 (MQ=255)

                                                                         tttCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAAt                <  1:4633711/93‑1 (MQ=255)

                                                                          ttCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATg               <  1:5541297/93‑1 (MQ=255)

                                                                          ttCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATg               >  1:2123259/1‑93 (MQ=255)

                                                                           tCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATgg              >  2:435060/1‑93 (MQ=255)

                                                                           tCAGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATgg              <  1:1877933/93‑1 (MQ=255)

                                                                             aGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGaa            >  1:2960670/1‑93 (MQ=255)

                                                                             aGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGaa            >  2:534731/1‑93 (MQ=255)

                                                                             aGGGAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGaa            >  2:1816746/1‑93 (MQ=255)

                                                                              gggAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAAc           >  1:5774346/1‑93 (MQ=255)

                                                                              gggAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAAc           >  1:2224253/1‑93 (MQ=255)

                                                                               ggAACTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACa          >  1:4894981/1‑93 (MQ=255)

                                                                                 accTTCCAGGCCTGCGTAAATTTCCCTGTACCAGCAAACAGGATCTGCGCGACAACGATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAg        <  1:399659/91‑1 (MQ=255)

                                                                                 aaCTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCTACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAg        <  2:3449969/93‑1 (MQ=255)

                                                                                 aaCTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAg        <  2:3276215/93‑1 (MQ=255)

                                                                                 aaCTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAg        >  2:5827087/1‑93 (MQ=255)

                                                                                   cTTTCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAgtg      >  2:5665976/1‑93 (MQ=255)

                                                                                     ttCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAgtggt    <  2:1234440/93‑1 (MQ=255)

                                                                                     ttCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGAAATCTATCCCTTTGACACCTTTGCCGTGCCAATGGAAGAAgtggt    <  1:5106561/93‑1 (MQ=255)

                                                                                      tCAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAgtggtg   >  2:805612/1‑93 (MQ=255)

                                                                                       cAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTgc  >  1:4413500/1‑93 (MQ=255)

                                                                                       cAGGCCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTgc  <  2:2245433/93‑1 (MQ=255)

CGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGC  >  NC_000913/1460234‑1460413

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