BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	817,654	+G	coding (377/513 nt)	moaB →	molybdopterin biosynthesis protein B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	817,651	1	.	G	100.0%	609.5 / NA	194	A125G (GCG→GGG)	moaB	molybdopterin biosynthesis protein B
Reads supporting (aligned to +/- strand): ref base . (0/0); new base G (102/92); total (102/92)

Only 100 of 213 total aligned reads displayed.

GTTCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGC‑GGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGCGCCGCA  >  NC_000913/817559‑817738

   cGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGC‑ggg                                                                                      <  2:2047767/93‑1 (MQ=255)

    gACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGC‑gggg                                                                                     <  1:4832461/93‑2 (MQ=255)

     aCCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGT‑ggcgg                                                                                    >  1:4743968/1‑90 (MQ=255)

     aCCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGC‑ggggc                                                                                    >  1:2747147/1‑91 (MQ=255)

       cGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGC‑ggggcgt                                                                                  >  2:4421185/1‑89 (MQ=255)

       cGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGC‑ggggcgt                                                                                  <  1:3872207/93‑5 (MQ=255)

             gTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGGCGCcaa                                                                             >  1:3918226/1‑93 (MQ=255)

               tGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCcaaca                                                                           >  1:2474483/1‑93 (MQ=255)

                 aaGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACaaa                                                                         >  1:2028186/1‑93 (MQ=255)

                   ggTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAAc                                                                       >  1:1495535/1‑93 (MQ=255)

                   ggTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAAc                                                                       >  1:6286078/1‑93 (MQ=255)

                      tttGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCt                                                                    <  1:8815/93‑1 (MQ=255)

                        tGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGa                                                                  >  2:3873494/1‑93 (MQ=255)

                         ggTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGAt                                                                 <  1:5497189/93‑1 (MQ=255)

                          gTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGAtt                                                                <  2:1249913/93‑1 (MQ=255)

                           tGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGAttt                                                               <  1:1190017/93‑1 (MQ=255)

                           ggAAGTGTTCCGGATGTTGTCGTGTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGAttt                                                               <  2:4613591/92‑1 (MQ=255)

                            gAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGAtttt                                                              >  2:318543/1‑93 (MQ=255)

                            gAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGAtttt                                                              >  1:4563596/1‑93 (MQ=255)

                             gaGTGTGCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTc                                                             <  2:5684065/92‑1 (MQ=255)

                             aaGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTc                                                             >  2:2481152/1‑93 (MQ=255)

                             aaGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTc                                                             >  2:530595/1‑93 (MQ=255)

                             aaGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTc                                                             <  2:642362/93‑1 (MQ=255)

                              aGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCg                                                            <  1:3305808/93‑1 (MQ=255)

                              aGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCg                                                            >  2:4699890/1‑93 (MQ=255)

                              aGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCg                                                            >  1:1745784/1‑93 (MQ=255)

                                tgtTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGAATTTCGcc                                                          >  1:100756/1‑93 (MQ=255)

                                 gtTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCa                                                         <  1:404338/93‑1 (MQ=255)

                                 gtTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCa                                                         <  1:2227976/93‑1 (MQ=255)

                                  ttCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCAt                                                        <  2:1020371/93‑1 (MQ=255)

                                  ttCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCAt                                                        >  2:4908241/1‑93 (MQ=255)

                                  ttCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCAt                                                        <  2:2238276/93‑1 (MQ=255)

                                   tCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATg                                                       >  2:627360/1‑93 (MQ=255)

                                    ccGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGc                                                      <  2:1125548/93‑1 (MQ=255)

                                    ccGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGc                                                      >  1:6026390/1‑93 (MQ=255)

                                      ggATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCg                                                    <  1:4538872/93‑1 (MQ=255)

                                       gATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCgg                                                   >  1:1517975/1‑93 (MQ=255)

                                       gATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCgg                                                   <  1:3788830/93‑1 (MQ=255)

                                       gATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCgg                                                   >  1:2449555/1‑93 (MQ=255)

                                       gATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCgg                                                   >  2:4185160/1‑93 (MQ=255)

                                       gATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCACGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCgg                                                   >  1:4435132/1‑93 (MQ=255)

                                         tgttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGGTGATTTTCGCCATGCCGGGt                                                 <  1:4064935/93‑1 (MQ=255)

                                         tgttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGt                                                 <  2:5090481/93‑1 (MQ=255)

                                         tgttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGt                                                 >  1:5006337/1‑93 (MQ=255)

                                         tgttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGt                                                 >  1:3410386/1‑93 (MQ=255)

                                         tgttgtCGTTTGAAGAGATTGGCACTCCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGt                                                 >  1:655107/1‑93 (MQ=255)

                                          gttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCgggct                                                >  2:1605569/1‑91 (MQ=255)

                                           ttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTc                                               >  1:3161372/1‑93 (MQ=255)

                                           ttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTc                                               >  1:1420085/1‑93 (MQ=255)

                                           ttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTc                                               >  2:356962/1‑93 (MQ=255)

                                           ttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTc                                               >  2:5997188/1‑93 (MQ=255)

                                           ttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTc                                               >  1:3434237/1‑93 (MQ=255)

                                              tCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGAc                                            >  2:4027641/1‑93 (MQ=255)

                                              tCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGAc                                            <  2:532711/93‑1 (MQ=255)

                                              tCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGAc                                            <  2:5563464/93‑1 (MQ=255)

                                              tCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGAc                                            <  2:2507146/93‑1 (MQ=255)

                                                   tGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAg                                       <  2:4932201/93‑1 (MQ=255)

                                                   tGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAg                                       >  2:2190686/1‑93 (MQ=255)

                                                       gagaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTg                                   <  2:2723978/93‑1 (MQ=255)

                                                       gagaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTg                                   >  2:2696628/1‑93 (MQ=255)

                                                         gaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGcc                                 >  2:6102728/1‑93 (MQ=255)

                                                         gaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGcc                                 >  1:2838235/1‑93 (MQ=255)

                                                         gaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGcc                                 <  2:1591411/93‑1 (MQ=255)

                                                         gaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGcc                                 <  2:679531/93‑1 (MQ=255)

                                                          aTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGCTTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCg                                >  2:4262013/1‑93 (MQ=255)

                                                           ttGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGt                               <  2:1915793/93‑1 (MQ=255)

                                                           ttGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGt                               >  1:1359494/1‑93 (MQ=255)

                                                            tGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTa                              >  1:2959138/1‑93 (MQ=255)

                                                             ggCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTAc                             <  1:1908255/93‑1 (MQ=255)

                                                             ggCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTAc                             <  2:294493/93‑1 (MQ=255)

                                                              gCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTAcc                            <  2:3859711/93‑1 (MQ=255)

                                                              gCACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATATTCGCCATGCCGGGTTCGAACAAAGCGTGCCgagcc                            >  2:3131739/1‑89 (MQ=255)

                                                               cACTTCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCg                           >  1:1534443/1‑93 (MQ=255)

                                                                  ttCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCAt                        >  2:5504442/1‑93 (MQ=255)

                                                                  ttCCACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCAt                        <  2:1377741/93‑1 (MQ=255)

                                                                    ccACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATgg                      >  2:3727575/1‑93 (MQ=255)

                                                                    ccACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATgg                      >  2:3997814/1‑93 (MQ=255)

                                                                     cACGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATggg                     >  1:3306458/1‑93 (MQ=255)

                                                                      aCGTTGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGa                    <  2:5749836/93‑1 (MQ=255)

                                                                         ttGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGaaaa                 >  2:922789/1‑93 (MQ=255)

                                                                         ttGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGaaaa                 >  2:3711597/1‑93 (MQ=255)

                                                                         ttGCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGaaaa                 >  2:3334355/1‑93 (MQ=255)

                                                                           gCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAata               <  2:824708/93‑1 (MQ=255)

                                                                           gCAATCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAata               >  1:5146576/1‑93 (MQ=255)

                                                                             aaTCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatc             <  1:1082801/93‑1 (MQ=255)

                                                                              aTCTCGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatca            <  2:3832279/93‑1 (MQ=255)

                                                                               tctcGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatcat           >  1:6080185/1‑93 (MQ=255)

                                                                               tctcGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatcat           <  1:5793530/93‑1 (MQ=255)

                                                                               tctcGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatcat           <  1:5034636/93‑1 (MQ=255)

                                                                                ctcGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatcatc          <  2:2867488/93‑1 (MQ=255)

                                                                                ctcGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatcatc          <  2:4737491/93‑1 (MQ=255)

                                                                                 tcGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATcg         <  1:6226702/93‑1 (MQ=255)

                                                                                  cGTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATcgc        <  2:5207940/93‑1 (MQ=255)

                                                                                   gTGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATcgcg       >  1:2934397/1‑93 (MQ=255)

                                                                                    tGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCgcgc      >  2:3077651/1‑93 (MQ=255)

                                                                                    tGCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCgcgc      <  1:4124337/93‑1 (MQ=255)

                                                                                     gCGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGcgcc     <  2:2456643/93‑1 (MQ=255)

                                                                                      cGGTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGcgccg    <  1:3940231/93‑1 (MQ=255)

                                                                                       ggTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGcgccgc   >  1:5572051/1‑93 (MQ=255)

                                                                                       ggTAGCGGGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGcgccgc   <  1:2080170/93‑1 (MQ=255)

GTTCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGC‑GGGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGCGCCGCA  >  NC_000913/817559‑817738

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