BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	780,918	T→C	intergenic (+43/‑390)	lysQ → / → nadA	tRNA‑Lys/quinolinate synthetase A protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	780,918	0	T	C	100.0%	541.9 / NA	222	intergenic (+43/‑390)	lysQ/nadA	tRNA‑Lys/quinolinate synthetase A protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (114/108); total (114/108)

Only 100 of 225 total aligned reads displayed.

TGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGC  >  NC_000913/780828‑780999

tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGacac                                                                                 <  2:3527446/93‑1 (MQ=21)

tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGacac                                                                                 <  1:1178151/93‑1 (MQ=21)

 gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACt                                                                                <  1:1333373/93‑1 (MQ=18)

   cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCa                                                                              <  1:6121778/93‑1 (MQ=21)

    ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCaa                                                                             <  2:5174159/93‑1 (MQ=21)

     tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCaaa                                                                            <  1:29986/93‑1 (MQ=21)

     tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCaaa                                                                            >  1:1087403/1‑93 (MQ=18)

     tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCaaa                                                                            <  1:5037319/93‑1 (MQ=18)

     tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCaaa                                                                            >  2:4623316/1‑93 (MQ=18)

      ttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAg                                                                           >  2:2289445/1‑93 (MQ=18)

       tAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGa                                                                          <  1:3077621/93‑1 (MQ=18)

       tAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGa                                                                          <  2:2100523/93‑1 (MQ=21)

       tAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGa                                                                          <  1:3413533/93‑1 (MQ=18)

        aaTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGAt                                                                         >  2:6223013/1‑93 (MQ=18)

          tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTc                                                                       <  2:7827268/93‑1 (MQ=21)

          tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTc                                                                       >  1:498483/1‑93 (MQ=18)

            aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGa                                                                     >  2:2062999/1‑93 (MQ=18)

            aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGa                                                                     >  2:4584192/1‑93 (MQ=18)

            aaTTGGTCGCAGGCTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAAATTCGa                                                                     >  2:84444/1‑93 (MQ=255)

             aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGaa                                                                    >  2:3956220/1‑93 (MQ=18)

               tGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATGTGCGACACTCAAAGATTCGAAcc                                                                  >  2:4281807/1‑93 (MQ=16)

               tGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAAcc                                                                  <  2:6298475/93‑1 (MQ=21)

                  tCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTgc                                                               >  2:4225394/1‑93 (MQ=18)

                    gCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTgcag                                                             >  2:3786916/1‑93 (MQ=18)

                     cAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTgcagc                                                            >  1:5186055/1‑93 (MQ=18)

                        gTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAgg                                                         >  2:1884225/1‑93 (MQ=18)

                        gTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAgg                                                         >  2:6647021/1‑93 (MQ=18)

                        gTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAgg                                                         <  1:878436/93‑1 (MQ=21)

                        gTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAgg                                                         >  1:3982952/1‑93 (MQ=18)

                        gTTCGAATCCTACACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAgg                                                         <  1:7480096/93‑1 (MQ=16)

                         ttCTAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGt                                                        <  2:1817907/93‑1 (MQ=16)

                             aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTttga                                                    <  2:2964358/93‑1 (MQ=18)

                              aTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTttgag                                                   >  2:1489951/1‑93 (MQ=18)

                               tCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTttgagt                                                  >  2:291747/1‑93 (MQ=18)

                               tCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTttgagt                                                  >  2:6972406/1‑93 (MQ=18)

                                   gCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                              <  1:6571486/93‑1 (MQ=21)

                                    cACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgc                                             <  1:1514445/93‑1 (MQ=18)

                                      cGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcgc                                           <  1:4730616/93‑1 (MQ=18)

                                       gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                          >  1:2335853/1‑93 (MQ=18)

                                        aCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAg                                         <  2:5925826/93‑1 (MQ=18)

                                          ccaccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCg                                       >  1:4736715/1‑93 (MQ=18)

                                          ccaccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCg                                       >  1:7876354/1‑93 (MQ=18)

                                           caccaATGTAAAAAAGCGCCCTAAAGGCTCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGa                                      <  1:2097844/93‑1 (MQ=16)

                                           caccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGa                                      >  2:7453365/1‑93 (MQ=18)

                                           caccaATGTAAAAAAGCGCCCTAAAGGCGATTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGa                                      >  1:3683454/1‑93 (MQ=16)

                                            accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaa                                     >  1:4216396/1‑93 (MQ=18)

                                             ccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                    <  2:7650126/93‑1 (MQ=18)

                                              caATGAAAAAAAGCGCCCTAAGGGCGTTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                   <  2:1681439/93‑1 (MQ=11)

                                               aaTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaaca                                  >  1:5467233/1‑93 (MQ=21)

                                               aaTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaaca                                  <  1:2826851/93‑1 (MQ=18)

                                                aTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaa                                 <  1:4570124/93‑1 (MQ=18)

                                                aTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaa                                 <  1:1591171/93‑1 (MQ=18)

                                                 tGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaac                                >  1:784435/1‑93 (MQ=21)

                                                 tGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaac                                >  2:3339594/1‑93 (MQ=18)

                                                   tAAAACAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACgg                              <  1:5338879/93‑1 (MQ=16)

                                                   tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACgg                              >  2:6505402/1‑93 (MQ=18)

                                                   tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACgg                              <  2:4890022/93‑1 (MQ=18)

                                                   tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACgg                              >  1:2113204/1‑93 (MQ=18)

                                                    aaaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGa                             <  2:2857992/93‑1 (MQ=18)

                                                    aaaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGa                             <  1:2203174/93‑1 (MQ=18)

                                                     aaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAg                            >  2:6852939/1‑93 (MQ=21)

                                                       aaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGcc                          <  2:387498/93‑1 (MQ=18)

                                                       aaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGcc                          >  1:2446874/1‑93 (MQ=21)

                                                       aaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGcc                          <  1:1105534/93‑1 (MQ=18)

                                                        aaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCg                         >  1:2771694/1‑93 (MQ=21)

                                                         aGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGc                        <  2:4000898/93‑1 (MQ=18)

                                                         aGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGc                        <  2:1245688/93‑1 (MQ=18)

                                                           cgcCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTc                      <  1:7741214/93‑1 (MQ=18)

                                                           cgcCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTc                      <  1:6374809/93‑1 (MQ=18)

                                                            gcCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcg                     >  1:5816305/1‑93 (MQ=21)

                                                            gcCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcg                     >  2:5024932/1‑93 (MQ=21)

                                                            gcCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcg                     >  1:4259487/1‑93 (MQ=18)

                                                             cccTAAAGGTGTTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                    <  1:245523/93‑1 (MQ=255)

                                                             cccTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                    <  1:6772137/93‑1 (MQ=18)

                                                               cTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcgg                  >  1:5130962/1‑93 (MQ=21)

                                                                tAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggc                 >  1:1030132/1‑93 (MQ=21)

                                                                tAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggc                 >  1:3137687/1‑93 (MQ=21)

                                                                  aaGGCTCTTTTTTGCTATCTGCGACTCTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGa               <  2:4732538/93‑1 (MQ=255)

                                                                  aaGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGa               <  2:5999095/93‑1 (MQ=18)

                                                                  aaGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGa               <  1:3932119/93‑1 (MQ=18)

                                                                   aGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc              <  1:7189888/93‑1 (MQ=18)

                                                                    ggCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACg             <  2:3127288/93‑1 (MQ=18)

                                                                    ggCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACg             <  1:3178362/93‑1 (MQ=18)

                                                                     gcgcTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg            >  2:3194692/1‑93 (MQ=21)

                                                                     gcgcTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg            <  2:2222492/93‑1 (MQ=18)

                                                                      cgcTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGc           <  1:6954010/93‑1 (MQ=18)

                                                                       gcTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGcc          >  2:7849847/1‑93 (MQ=21)

                                                                       gcTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGcc          >  1:6582015/1‑93 (MQ=21)

                                                                        cTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGccc         >  2:5479099/1‑93 (MQ=21)

                                                                          tttttGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGa       >  1:5215060/1‑93 (MQ=21)

                                                                          tttttGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGa       >  1:7308123/1‑93 (MQ=21)

                                                                            tttGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAg     <  1:2018843/93‑1 (MQ=18)

                                                                            tttGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAg     >  1:1003196/1‑93 (MQ=18)

                                                                            tttGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAg     <  2:3698550/93‑1 (MQ=18)

                                                                            tttGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAg     <  2:3310342/93‑1 (MQ=18)

                                                                             ttGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAgg    <  2:5358878/93‑1 (MQ=18)

                                                                              tGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg   >  2:6077435/1‑93 (MQ=17)

                                                                              tGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg   >  2:5282992/1‑93 (MQ=17)

                                                                              tGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg   <  2:2466539/93‑1 (MQ=17)

                                                                               gCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgc  <  2:840820/93‑1 (MQ=14)

TGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGC  >  NC_000913/780828‑780999

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: `‑`