BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,923,857	T→C	Q76R (CAG→CGG)	mioC ←	MioC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,923,857	0	T	C	99.8%	1703.2 / ‑5.3	522	Q76R (CAG→CGG)	mioC	MioC
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (219/302); total (220/302)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 532 total aligned reads displayed.

TCGAGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCTGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAG  >  NC_000913/3923767‑3923946

tCGAGTTTATCGATAGCCCCACACAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGt                                                                                         <  2:5494718/93‑1 (MQ=255)

tCGAGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGt                                                                                         >  2:8716109/1‑93 (MQ=255)

tCGAGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGt                                                                                         <  1:2470385/93‑1 (MQ=255)

   aGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTcc                                                                                      <  2:5910723/93‑1 (MQ=255)

   aGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTcc                                                                                      >  1:9029632/1‑93 (MQ=255)

     tttATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTg                                                                                    >  2:2777689/1‑93 (MQ=255)

       tATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCa                                                                                  >  2:1268765/1‑93 (MQ=255)

       tATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCa                                                                                  <  1:8304604/93‑1 (MQ=255)

       tATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCa                                                                                  >  2:2490848/1‑93 (MQ=255)

        agcGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCaa                                                                                 <  1:8496083/91‑1 (MQ=255)

        aTCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCaa                                                                                 >  1:4207403/1‑93 (MQ=255)

           gATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGc                                                                              <  1:9247793/93‑1 (MQ=255)

            aTAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCt                                                                             <  1:7386544/93‑1 (MQ=255)

            aTAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCt                                                                             >  2:8491904/1‑93 (MQ=255)

            aTAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCt                                                                             <  2:1105395/93‑1 (MQ=255)

             tAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCTGTTCCTGCAATGCtt                                                                            >  1:4905346/1‑93 (MQ=255)

             tAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCtt                                                                            <  2:7496904/93‑1 (MQ=255)

             tAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCtt                                                                            >  1:1248324/1‑93 (MQ=255)

             catccccACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCtt                                                                            >  1:6123781/4‑93 (MQ=255)

              aGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTc                                                                           <  2:3617606/93‑1 (MQ=255)

               gCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCa                                                                          <  2:7754785/93‑1 (MQ=255)

               gCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCa                                                                          <  2:5076267/93‑1 (MQ=255)

               gCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCa                                                                          >  2:6439453/1‑93 (MQ=255)

                ccccACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCAt                                                                         >  2:6674106/1‑93 (MQ=255)

                 cccACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATa                                                                        <  1:3023174/93‑1 (MQ=255)

                  ccACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATag                                                                       >  2:4169193/1‑93 (MQ=255)

                  ccACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATag                                                                       >  2:7999099/1‑93 (MQ=255)

                   cacaAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATaga                                                                      >  1:5767491/1‑93 (MQ=255)

                   cacaAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATaga                                                                      <  1:8675177/93‑1 (MQ=255)

                   cacaAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATaga                                                                      >  2:4891335/1‑93 (MQ=255)

                    acaAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGATAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaa                                                                     >  1:3478558/1‑93 (MQ=255)

                    acaAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaa                                                                     <  1:5197459/93‑1 (MQ=255)

                       aaaGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATcgaaaga                                                                  >  1:2269968/1‑92 (MQ=255)

                       aaaGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaagg                                                                  >  2:6753631/1‑93 (MQ=255)

                         aGGTTTCATATTCACGACTGCTAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggag                                                                <  1:10014048/93‑1 (MQ=255)

                          ggTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggaga                                                               <  1:6708325/93‑1 (MQ=255)

                          ggTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggaga                                                               <  1:7908651/93‑1 (MQ=255)

                          ggTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggaga                                                               <  1:6288080/93‑1 (MQ=255)

                          ggTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggaga                                                               >  1:2746757/1‑93 (MQ=255)

                           gtgtCATATTCACGACTGCCAATACCGATTGGGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggagaa                                                              <  1:3333067/93‑1 (MQ=255)

                            ggtCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggagaaa                                                             <  1:7987745/92‑1 (MQ=255)

                             gtCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggagaaag                                                            >  2:1513271/1‑93 (MQ=255)

                               cATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGt                                                          <  1:7170479/93‑1 (MQ=255)

                               cATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGt                                                          <  1:5520787/93‑1 (MQ=255)

                               cATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGt                                                          <  1:204637/93‑1 (MQ=255)

                                atatTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGtt                                                         <  1:6579444/93‑1 (MQ=255)

                                 tatTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTg                                                        >  2:697932/1‑93 (MQ=255)

                                 tatTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTg                                                        >  2:8522992/1‑93 (MQ=255)

                                   ttCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTc                                                      >  2:7120060/1‑93 (MQ=255)

                                    tCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTcc                                                     <  2:985402/93‑1 (MQ=255)

                                    tCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTcc                                                     >  1:3888108/1‑93 (MQ=255)

                                       cGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGa                                                  <  1:5326465/93‑1 (MQ=255)

                                        gACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGaa                                                 <  1:3116300/93‑1 (MQ=255)

                                          cTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAata                                               <  2:8910672/93‑1 (MQ=255)

                                           tGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAatat                                              >  1:6016679/1‑93 (MQ=255)

                                            gCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATAtc                                             <  1:471173/93‑1 (MQ=255)

                                             ccAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATAtct                                            <  1:9483852/93‑1 (MQ=255)

                                              cAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATAtctc                                           <  2:1316457/93‑1 (MQ=255)

                                              cAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATAtctc                                           >  1:1174044/1‑93 (MQ=255)

                                               aaTACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTcc                                          <  2:9720166/93‑1 (MQ=255)

                                                aTACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCg                                         >  2:2720851/1‑93 (MQ=255)

                                                aTACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCg                                         <  2:2989036/93‑1 (MQ=255)

                                                aTACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCg                                         <  2:3109795/93‑1 (MQ=255)

                                                 tACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCgg                                        <  2:9009692/93‑1 (MQ=255)

                                                   ccGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGcc                                      >  2:4263655/1‑92 (MQ=255)

                                                   ccGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCa                                      <  1:375654/93‑1 (MQ=255)

                                                    cGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCAc                                     <  2:9557378/93‑1 (MQ=255)

                                                     gATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCAcc                                    <  1:4612581/93‑1 (MQ=255)

                                                        tGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCgtg                                 <  2:3961228/93‑1 (MQ=255)

                                                         gcgcCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCgtgg                                <  2:3287128/93‑1 (MQ=255)

                                                         gcgcCAAAGCGGACGGCAGAAAGATCCGGCTTCCGTTCCTGCCATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCgtgg                                <  1:1625336/93‑1 (MQ=255)

                                                          cgcCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCgtggg                               <  2:9272386/93‑1 (MQ=255)

                                                          cgcCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCgtggg                               <  2:2568984/93‑1 (MQ=255)

                                                           gcCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCgtgggt                              <  1:6870923/93‑1 (MQ=255)

                                                           gcCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCgtgggt                              <  2:2289110/93‑1 (MQ=255)

                                                           gcCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAAAGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCgtgggt                              <  2:5195445/93‑1 (MQ=255)

                                                             cAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCgtgggtgg                            <  2:5355277/93‑1 (MQ=255)

                                                              aaaGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGa                           >  2:1757894/1‑93 (MQ=255)

                                                              aaaGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGa                           <  1:5873651/93‑1 (MQ=255)

                                                              aaaGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGa                           >  2:3397975/1‑93 (MQ=255)

                                                              aaaGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGa                           <  2:2813299/93‑1 (MQ=255)

                                                               aaGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAg                          <  1:7591125/93‑1 (MQ=255)

                                                               aaGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAg                          >  2:8257578/1‑93 (MQ=255)

                                                                aGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGc                         >  2:1608873/1‑93 (MQ=255)

                                                                 gCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTTATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCt                        >  1:9724916/1‑93 (MQ=255)

                                                                 gCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCt                        >  1:6489520/1‑93 (MQ=255)

                                                                   ggACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGa                      >  2:570491/1‑93 (MQ=255)

                                                                    gACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGAt                     <  1:891202/93‑1 (MQ=255)

                                                                      cTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATaa                   <  1:4415173/93‑1 (MQ=255)

                                                                      cTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATaa                   >  1:1752667/1‑93 (MQ=255)

                                                                       tGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAAc                  <  2:7349666/93‑1 (MQ=255)

                                                                        gCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAAcc                 <  1:1387242/93‑1 (MQ=255)

                                                                              aGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAAccagccag           >  2:8208023/1‑93 (MQ=255)

                                                                               gATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGa          >  1:8782953/1‑93 (MQ=255)

                                                                                     gCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTg    <  1:3665939/93‑1 (MQ=255)

                                                                                     gCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTg    >  1:4090702/1‑93 (MQ=255)

                                                                                     gCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTg    <  2:6109581/93‑1 (MQ=255)

                                                                                     gCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTg    >  2:4541321/1‑93 (MQ=255)

                                                                                       ttCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAg  <  2:2021337/93‑1 (MQ=255)

                                                                                       ttCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAg  >  2:4721649/1‑93 (MQ=255)

TCGAGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCTGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAG  >  NC_000913/3923767‑3923946

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