BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,365,344	A→G	pseudogene (160/371 nt)	yhcE →	putative uncharacterized protein YhcE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,365,344	0	A	G	100.0%	1558.3 / NA	370	pseudogene (160/371 nt)	yhcE	putative uncharacterized protein YhcE
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (204/166); total (204/166)

Only 100 of 375 total aligned reads displayed.

GTAATATCTCCAGTAGCATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAACAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAA  >  NC_000913/3365246‑3365441

 tAATATCTCCAGTAGCATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCaa                                                                                                >  2:6485947/1‑101 (MQ=255)

 tAATATCTCCAGTAGCATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCAAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCaa                                                                                                >  1:6485946/1‑101 (MQ=255)

   atatCTCCAGTAGCATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATc                                                                                              >  1:6485954/1‑101 (MQ=255)

    tatCTCCAGTAGCATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCa                                                                                             >  2:6485955/1‑101 (MQ=255)

     atCTCCAGTAGCATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCAc                                                                                            <  2:6485116/101‑1 (MQ=255)

     atCTCCAGTAGCATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCAc                                                                                            <  1:6485375/101‑1 (MQ=255)

       ctcCAGTAGCATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGc                                                                                          >  2:6485963/1‑101 (MQ=255)

       ctcCAGTAGCATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGc                                                                                          <  1:6485156/101‑1 (MQ=255)

        tcCAGTAGCATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCt                                                                                         >  1:6485965/1‑101 (MQ=255)

        tcCAGTAGCATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCt                                                                                         <  1:6485628/101‑1 (MQ=255)

           aGTAGCATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGt                                                                                      <  2:6485285/101‑1 (MQ=255)

            gTAGCATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTa                                                                                     >  1:6485970/1‑101 (MQ=255)

             tAGCATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTaa                                                                                    <  1:6485401/101‑1 (MQ=255)

             tAGCATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAACAGACAGAGCAATCACGCTGGTaa                                                                                    <  2:6485143/101‑1 (MQ=255)

              aGCATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTaat                                                                                   >  2:6485972/1‑101 (MQ=255)

               gCATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTaatg                                                                                  <  2:6485350/101‑1 (MQ=255)

                cATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTaatgc                                                                                 >  1:6485976/1‑101 (MQ=255)

                 aTACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTaatgca                                                                                >  1:6485979/1‑101 (MQ=255)

                  tACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTaatgcaa                                                                               <  1:6484898/101‑1 (MQ=255)

                     aaaaTATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTaatgcaatgc                                                                            >  1:6485986/1‑101 (MQ=255)

                     aaaaTATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTaatgcaatgc                                                                            >  2:6485985/1‑101 (MQ=255)

                       aaTATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGa                                                                          >  2:6485990/1‑101 (MQ=255)

                       aaTATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGa                                                                          <  1:6485268/101‑1 (MQ=255)

                         tatCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATc                                                                        >  2:6485994/1‑101 (MQ=255)

                          atCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCa                                                                       >  2:6485999/1‑101 (MQ=255)

                           tCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAg                                                                      >  1:6486002/1‑101 (MQ=255)

                            cGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAgg                                                                     <  1:6485491/101‑1 (MQ=255)

                            cGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAgg                                                                     >  1:6486007/1‑101 (MQ=255)

                             gACTATGGCAAACGAAGTGGTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGa                                                                    >  2:6486011/1‑101 (MQ=255)

                                tATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGAc                                                                   <  1:6485014/99‑1 (MQ=255)

                                tATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGc                                                                 >  1:6486014/1‑101 (MQ=255)

                                 aTGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCt                                                                >  1:6486023/1‑101 (MQ=255)

                                 aTGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCt                                                                >  2:6486019/1‑101 (MQ=255)

                                 aTGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCt                                                                <  1:6485229/101‑1 (MQ=255)

                                    gCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCat                                                             >  1:6486027/1‑101 (MQ=255)

                                    gCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCat                                                             <  1:6485351/101‑1 (MQ=255)

                                      aaaCGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCatat                                                           >  2:6486030/1‑101 (MQ=255)

                                       aaCGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCatatc                                                          <  1:6485363/101‑2 (MQ=255)

                                       aaCGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATAtt                                                          >  2:6486035/1‑101 (MQ=255)

                                       aaCGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATAtt                                                          >  1:6486034/1‑101 (MQ=255)

                                         cGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCg                                                        >  1:6486040/1‑101 (MQ=255)

                                          gAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCgt                                                       >  2:6486044/1‑101 (MQ=255)

                                           aaGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCgtg                                                      >  1:6486046/1‑101 (MQ=255)

                                             gTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGtt                                                    <  2:6485750/101‑1 (MQ=255)

                                              tgctgcCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTc                                                   >  2:6486053/1‑101 (MQ=255)

                                              tgctgcCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTc                                                   <  2:6485386/101‑1 (MQ=255)

                                              tgctgcCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTc                                                   <  1:6485450/101‑1 (MQ=255)

                                               gctgcCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCa                                                  <  2:6485371/101‑1 (MQ=255)

                                                ctgcCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAg                                                 >  1:6486055/1‑101 (MQ=255)

                                                 tgcCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGt                                                <  2:6485507/101‑1 (MQ=255)

                                                 tgcCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGt                                                >  2:6486060/1‑101 (MQ=255)

                                                  gcCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGtt                                               >  1:6486064/1‑101 (MQ=255)

                                                  gcCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGtt                                               <  1:6485595/101‑1 (MQ=255)

                                                   ccATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGa                                                                    <  2:6485197/79‑1 (MQ=255)

                                                   ccATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTg                                              <  2:6485441/101‑1 (MQ=255)

                                                    cATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGa                                             <  2:6485305/101‑1 (MQ=255)

                                                    cATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGa                                             >  2:6486066/1‑101 (MQ=255)

                                                    cATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGa                                             >  1:6486070/1‑101 (MQ=255)

                                                      tGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAAt                                           <  1:6485339/101‑1 (MQ=255)

                                                       gcgcCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATa                                          <  1:6485121/101‑1 (MQ=255)

                                                       gcgcCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATa                                          <  2:6485533/101‑1 (MQ=255)

                                                        cgcCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATAc                                         >  2:6486074/1‑101 (MQ=255)

                                                         gcCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATAcc                                        >  1:6486077/1‑101 (MQ=255)

                                                           cAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGc                                      >  2:6486082/1‑101 (MQ=255)

                                                            aaGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCt                                     <  1:6485575/101‑1 (MQ=255)

                                                             aGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTa                                    <  1:6485392/101‑1 (MQ=255)

                                                             aGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTa                                    >  1:6486085/1‑101 (MQ=255)

                                                               ttGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAat                                  <  2:6485478/101‑1 (MQ=255)

                                                               ttGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAat                                  >  2:6486089/1‑101 (MQ=255)

                                                                tGACCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAata                                 >  2:6486092/1‑101 (MQ=255)

                                                                  aCCGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATAtt                               <  1:6485681/101‑1 (MQ=255)

                                                                   ccGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATAttt                              >  2:6486100/1‑101 (MQ=255)

                                                                   ccGGGGTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATAttt                              >  1:6486097/1‑101 (MQ=255)

                                                                     ggggTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTcc                            <  1:6485203/101‑1 (MQ=255)

                                                                      gggTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCca                           >  1:6486105/1‑101 (MQ=255)

                                                                       ggTAAAACAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCcaa                          >  2:6486106/1‑101 (MQ=255)

                                                                          aaaaCAACACAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCcaacaa                       >  2:6486110/1‑101 (MQ=255)

                                                                             acaacaCAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACgg                    >  1:6486115/1‑101 (MQ=255)

                                                                             acaacaCAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACgg                    <  1:6485426/101‑1 (MQ=255)

                                                                             acaacaCAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACgg                    <  2:6485663/101‑1 (MQ=255)

                                                                              caacaCAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACgga                   <  1:6485312/101‑1 (MQ=255)

                                                                               aacaCAATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACggat                  >  1:6486119/1‑101 (MQ=255)

                                                                                  acaATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACggatttg               >  2:6486123/1‑101 (MQ=255)

                                                                                   caATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACggatttgg              >  1:6486126/1‑101 (MQ=255)

                                                                                   caATTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACgg                    <  2:6485250/95‑1 (MQ=255)

                                                                                     aTTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACggatttggat            <  1:6485504/101‑1 (MQ=255)

                                                                                     aTTAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACggatttggat            <  1:6485763/101‑1 (MQ=255)

                                                                                      ttAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACggatttggatt           >  2:6486131/1‑101 (MQ=255)

                                                                                      ttAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACggatttggatt           <  1:6485213/101‑1 (MQ=255)

                                                                                       tAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACggatttggattt          <  2:6485613/101‑1 (MQ=255)

                                                                                       tAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACggatttggattt          <  1:6485534/101‑1 (MQ=255)

                                                                                       tAGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACggatttggattt          >  1:6486134/1‑101 (MQ=255)

                                                                                        aGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGgatttggatttg         >  1:6486139/1‑101 (MQ=255)

                                                                                        aGCAGACAGAGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGgatttggatttg         <  1:6485328/101‑1 (MQ=255)

                                                                                          cagacagaGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGc       >  2:6486141/1‑101 (MQ=255)

                                                                                          cagacagaGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGc       <  2:6485403/101‑1 (MQ=255)

                                                                                           agacagaGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGcc      >  1:6486145/1‑101 (MQ=255)

                                                                                            gacagaGCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGccc     <  1:6485474/101‑1 (MQ=255)

                                                                                            gacagaGCAATCACGCTGATAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGccc     <  1:6485703/101‑1 (MQ=255)

                                                                                               agagCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGaa  >  1:6486150/1‑101 (MQ=255)

GTAATATCTCCAGTAGCATACAAAATATCGACTATGGCAAACGAAGTGCTGCCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAACAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAA  >  NC_000913/3365246‑3365441

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`