BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,217	T→C	D627D (GAT→GAC)	thrA →	aspartokinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,217	0	T	C	99.7%	1028.4 / ‑5.2	314	D627D (GAT→GAC)	thrA	aspartokinase I
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (124/189); total (125/189)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

Only 100 of 323 total aligned reads displayed.

AAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGAC  >  NC_000913/2127‑2306

 aaaTCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGaa                                                                                        <  1:4617413/93‑1 (MQ=255)

 aaaTCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGaa                                                                                        <  2:520524/93‑1 (MQ=255)

 aaaTCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGaa                                                                                        <  2:527418/93‑1 (MQ=255)

   aTCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAAtt                                                                                      <  1:5989892/93‑1 (MQ=255)

   aTCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAAtt                                                                                      <  1:2257811/93‑1 (MQ=255)

    tCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATtg                                                                                     >  1:5231551/1‑93 (MQ=255)

      gcggcgTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATtgat                                                                                   <  2:3729526/93‑1 (MQ=255)

      gcggcgTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATtgat                                                                                   <  1:2975466/93‑1 (MQ=255)

       cggcgTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATtgatg                                                                                  <  2:4866024/93‑1 (MQ=255)

         gcgTCAATTCCTCTATGACACCCACGTGTGGTCTGGATTACCGGTTATTGAGAACGTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGaa                                                                                <  1:1122182/93‑1 (MQ=255)

          cgTAAATTCCTCTATGACACCAATGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAg                                                                               <  1:1759871/93‑1 (MQ=255)

           gTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGt                                                                              >  2:3049633/1‑93 (MQ=255)

           gTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGt                                                                              <  1:2381768/93‑1 (MQ=255)

             aaaTTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTtc                                                                            <  2:1215418/93‑1 (MQ=255)

             aaaTTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTtc                                                                            <  2:1383986/93‑1 (MQ=255)

              aaTTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTtct                                                                           <  1:5521151/93‑1 (MQ=255)

               aTTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTtctc                                                                          <  2:3000270/93‑1 (MQ=255)

               aTTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTtctc                                                                          <  1:1002620/93‑1 (MQ=255)

               aTTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTtctc                                                                          >  1:1820163/1‑93 (MQ=255)

                ttCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTcc                                                                         >  2:4622293/1‑93 (MQ=255)

                ttCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTcc                                                                         <  2:3664215/93‑1 (MQ=255)

                   ctctATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGc                                                                      >  2:5576030/1‑93 (MQ=255)

                    tctATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCa                                                                     >  1:2431794/1‑93 (MQ=255)

                     ctATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCAt                                                                    <  1:6492375/93‑1 (MQ=255)

                     ctATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCAt                                                                    >  1:1040311/1‑93 (MQ=255)

                      tATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCAtt                                                                   <  2:5840444/93‑1 (MQ=255)

                      tATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCAtt                                                                   <  2:4446837/93‑1 (MQ=255)

                       aTGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCAttc                                                                  <  2:3822653/93‑1 (MQ=255)

                       aTGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCAttc                                                                  >  2:6653031/1‑93 (MQ=255)

                       aTGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCAttc                                                                  <  2:4210227/93‑1 (MQ=255)

                        tGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCAttct                                                                 <  2:2848924/93‑1 (MQ=255)

                         gACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCAttctt                                                                >  2:2341425/1‑93 (MQ=255)

                           cacCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCAttctttc                                                              <  2:2366692/93‑1 (MQ=255)

                            acCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCAttctttct                                                             <  1:2561774/93‑1 (MQ=255)

                            acCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCAttctttct                                                             <  1:5858361/93‑1 (MQ=255)

                              cAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTgg                                                           <  2:6597199/93‑1 (MQ=255)

                               aaCGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGt                                                          <  2:1544068/93‑1 (MQ=255)

                                aCGTTGGGTCTGGATTAGCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGtt                                                         <  1:1968522/93‑1 (MQ=255)

                                aCGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGtt                                                         >  1:284449/1‑93 (MQ=255)

                                 cGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTc                                                        >  1:6185666/1‑93 (MQ=255)

                                   ttGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGc                                                      >  2:3289172/1‑93 (MQ=255)

                                   ttGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGc                                                      <  2:4924615/93‑1 (MQ=255)

                                   ttGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGc                                                      >  1:1424299/1‑93 (MQ=255)

                                   ttGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGc                                                      >  1:5657946/1‑93 (MQ=255)

                                     ggggCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCtt                                                    <  2:3928752/93‑1 (MQ=255)

                                      gggCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCttt                                                   <  2:2086528/93‑1 (MQ=255)

                                      gggCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCttt                                                   <  2:1460719/93‑1 (MQ=255)

                                        gCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCt                                                 <  2:3551996/93‑1 (MQ=255)

                                        gCTAGATTACCTGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATTAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCt                                                 <  1:4058111/93‑1 (MQ=255)

                                          tGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTta                                               >  1:5558442/1‑93 (MQ=255)

                                           ggATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTtat                                              >  1:986569/1‑93 (MQ=255)

                                            gATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTtata                                             >  1:510601/1‑93 (MQ=255)

                                             aTTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTatat                                            >  2:5480420/1‑93 (MQ=255)

                                              ttACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATc                                           <  2:1831987/93‑1 (MQ=255)

                                               tACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCt                                          >  2:224602/1‑93 (MQ=255)

                                               tACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCt                                          <  1:4349781/93‑1 (MQ=255)

                                                  cGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCg                                       <  1:4072837/93‑1 (MQ=255)

                                                  cGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCg                                       >  2:6315908/1‑93 (MQ=255)

                                                    gTTATTGAGAACTTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGc                                     <  2:4283841/93‑1 (MQ=255)

                                                     ttATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATATTCGGCa                                    >  1:3700690/1‑93 (MQ=255)

                                                     gtATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCa                                    <  1:4954990/92‑1 (MQ=255)

                                                      tATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCaa                                   >  2:6209015/1‑93 (MQ=255)

                                                      tATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCaa                                   <  1:6577738/93‑1 (MQ=255)

                                                       aTTGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAg                                  >  1:548602/1‑93 (MQ=255)

                                                        ttGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGt                                 <  1:2129703/93‑1 (MQ=255)

                                                        ttGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGt                                 <  1:1547864/93‑1 (MQ=255)

                                                         tGAGAACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGtt                                <  1:3555776/93‑1 (MQ=255)

                                                          gagaACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTa                               >  2:3122480/1‑93 (MQ=255)

                                                            gatccTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGa                             <  2:1987711/90‑1 (MQ=255)

                                                            gaACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGa                             >  1:3980699/1‑93 (MQ=255)

                                                            gaACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGa                             >  1:632865/1‑93 (MQ=255)

                                                            gaACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGa                             >  2:4689634/1‑93 (MQ=255)

                                                            gaACCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGa                             >  1:4779349/1‑93 (MQ=255)

                                                              aCCTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACg                           <  1:4495663/93‑1 (MQ=255)

                                                               ccTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGa                          <  2:6450845/93‑1 (MQ=255)

                                                                cTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGaa                         <  2:6007849/93‑1 (MQ=255)

                                                                cTGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGaa                         <  2:2451645/93‑1 (MQ=255)

                                                                 tGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAg                        >  2:5029696/1‑93 (MQ=255)

                                                                 tGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAg                        <  2:4104592/93‑1 (MQ=255)

                                                                 tGCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAg                        <  2:288048/93‑1 (MQ=255)

                                                                  gCAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAgg                       <  2:2763550/93‑1 (MQ=255)

                                                                   cAAAATCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGc                      <  1:3314543/93‑1 (MQ=255)

                                                                    aaaaTCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGCAGGCa                     >  1:3080453/1‑93 (MQ=255)

                                                                     aaaTCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCAt                    >  1:607641/1‑93 (MQ=255)

                                                                     aaaTCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCAt                    <  2:3492465/93‑1 (MQ=255)

                                                                      aaTCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATg                   >  2:816336/1‑93 (MQ=255)

                                                                       aTCTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGa                  <  1:3385134/93‑1 (MQ=255)

                                                                         cTGCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGt                >  1:3153261/1‑93 (MQ=255)

                                                                           gCTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGttt              <  1:3935825/93‑1 (MQ=255)

                                                                            cTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTtc             <  2:4653726/93‑1 (MQ=255)

                                                                            cTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTtc             <  2:3852544/93‑1 (MQ=255)

                                                                            cTCAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTtc             <  1:5048063/93‑1 (MQ=255)

                                                                              cAATGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTtctc           >  2:2595782/1‑93 (MQ=255)

                                                                               aaTGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTcc          >  2:2833423/1‑93 (MQ=255)

                                                                                aTGCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCg         >  2:2420056/1‑93 (MQ=255)

                                                                                  gCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAg       <  1:3506964/93‑1 (MQ=255)

                                                                                  gCAGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAg       <  1:6693229/93‑1 (MQ=255)

                                                                                   cAGGTGACGAATTGATAAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAgg      <  1:2391720/93‑1 (MQ=255)

                                                                                    aGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGc     >  2:3187498/1‑93 (MQ=255)

                                                                                    aGGTGACGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGc     <  2:2300161/93‑1 (MQ=255)

AAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGAC  >  NC_000913/2127‑2306

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