BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	985,663	A→G	G181G (GGT→GGC)	ompF ←	outer membrane protein 1a

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	985,663	0	A	G	100.0%	813.8 / NA	256	G181G (GGT→GGC)	ompF	outer membrane protein 1a
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (127/129); total (127/129)

Only 100 of 264 total aligned reads displayed.

ACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAAC  >  NC_000913/985573‑985752

aCCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGcc                                                                                         <  2:4181237/93‑1 (MQ=255)

 ccAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGccc                                                                                        >  2:5691181/1‑93 (MQ=255)

  cAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCa                                                                                       <  1:4338638/93‑1 (MQ=255)

   aaCGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAg                                                                                      <  1:3157644/93‑1 (MQ=255)

   aaCGATACCAAAGCCTTCGTATTCGTAGCTCATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAg                                                                                      <  2:2509163/93‑1 (MQ=255)

    aCGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAgg                                                                                     <  1:4467985/93‑1 (MQ=255)

    aCGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAgg                                                                                     <  2:290557/93‑1 (MQ=255)

    aCGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAgg                                                                                     <  1:2817605/93‑1 (MQ=255)

     cGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGt                                                                                    <  2:548563/93‑1 (MQ=255)

      gATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTa                                                                                   <  2:3837404/93‑1 (MQ=255)

         aCCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTg                                                                                >  2:2060858/1‑93 (MQ=255)

          ccAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGa                                                                               >  1:3730199/1‑93 (MQ=255)

          ccAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGa                                                                               <  1:3348302/93‑1 (MQ=255)

           cAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGaa                                                                              <  2:946403/93‑1 (MQ=255)

           cAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGaa                                                                              >  1:1103274/1‑93 (MQ=255)

           cAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGaa                                                                              >  2:5033638/1‑93 (MQ=255)

             aaGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACTGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACa                                                                            <  2:5136353/93‑1 (MQ=255)

             aaGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACa                                                                            >  2:1132372/1‑93 (MQ=255)

             aaGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACa                                                                            <  1:3845821/93‑1 (MQ=255)

              aGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAg                                                                           <  1:6149825/93‑1 (MQ=255)

                 cTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGa                                                                        <  2:3812390/93‑1 (MQ=255)

                 cTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGa                                                                        >  1:4063911/1‑93 (MQ=255)

                  ttcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGaa                                                                       <  1:454428/93‑1 (MQ=255)

                  ttcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGaa                                                                       <  2:122495/93‑1 (MQ=255)

                    cgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGt                                                                     <  2:1669674/93‑1 (MQ=255)

                    cgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGt                                                                     >  2:3059587/1‑93 (MQ=255)

                     gtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGtt                                                                    <  1:1723732/93‑1 (MQ=255)

                      tattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTc                                                                   >  2:562252/1‑93 (MQ=255)

                      tattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTc                                                                   >  1:6227731/1‑93 (MQ=255)

                       attcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCa                                                                  >  1:3033278/1‑93 (MQ=255)

                        ttcgtaGCTGATAGACCCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAg                                                                 <  1:5426129/93‑1 (MQ=255)

                        ttcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAg                                                                 >  2:3349125/1‑93 (MQ=255)

                        ttcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAg                                                                 >  2:1522305/1‑93 (MQ=255)

                         tcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAgg                                                                >  1:356262/1‑93 (MQ=255)

                           gtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGcc                                                              >  2:5084409/1‑93 (MQ=255)

                            taGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCa                                                             >  1:875484/1‑93 (MQ=255)

                             aGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCAt                                                            >  1:5823229/1‑93 (MQ=255)

                              gCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATc                                                           <  1:3721318/93‑1 (MQ=255)

                              gCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATc                                                           <  2:2712620/93‑1 (MQ=255)

                                tGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCaa                                                         >  1:2790857/1‑93 (MQ=255)

                                 gATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAc                                                        >  1:2113607/1‑93 (MQ=255)

                                 gATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAc                                                        >  1:702224/1‑93 (MQ=255)

                                   tAGACCCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCa                                                      <  1:3971508/93‑1 (MQ=255)

                                   tAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCa                                                      <  2:2237518/93‑1 (MQ=255)

                                   tAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCa                                                      >  2:2157408/1‑93 (MQ=255)

                                   tAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCa                                                      <  2:5821279/93‑1 (MQ=255)

                                     gAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGa                                                    <  2:87787/93‑1 (MQ=255)

                                     gAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGa                                                    >  2:1413708/1‑93 (MQ=255)

                                     gAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCACAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGa                                                    >  1:3707916/1‑93 (MQ=255)

                                       aCCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAcc                                                  >  1:1208555/1‑93 (MQ=255)

                                       aCCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAcc                                                  <  2:4568700/93‑1 (MQ=255)

                                       aCCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGGCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAcc                                                  >  2:2805667/1‑93 (MQ=255)

                                        ccGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCa                                                 >  1:3454290/1‑93 (MQ=255)

                                        ccGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCa                                                 <  2:1031850/93‑1 (MQ=255)

                                         cGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCaa                                                <  2:403769/93‑1 (MQ=255)

                                          gCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCaaa                                               >  2:660697/1‑93 (MQ=255)

                                          gCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCaaa                                               >  2:4403524/1‑93 (MQ=255)

                                          gCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCaaa                                               >  2:2186474/1‑93 (MQ=255)

                                            cAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCGCGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaaga                                             <  1:1519429/93‑1 (MQ=255)

                                             aaCACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaagaa                                            >  1:1357402/1‑93 (MQ=255)

                                              acacCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaagaag                                           >  2:3599751/1‑93 (MQ=255)

                                              acacCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaagaag                                           >  1:1788415/1‑93 (MQ=255)

                                                acCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGtt                                         <  2:2562977/93‑1 (MQ=255)

                                                acCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGtt                                         <  2:1604067/93‑1 (MQ=255)

                                                 ccGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTg                                        <  1:3867640/93‑1 (MQ=255)

                                                  cGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTgg                                       <  1:5500436/93‑1 (MQ=255)

                                                   gTCGCCGTGAGAACGGCGTGCAGTCTCACCCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGa                                      <  2:2337521/93‑1 (MQ=25)

                                                     cGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGt                                    <  2:3193858/93‑1 (MQ=255)

                                                     cGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGt                                    <  1:768054/93‑1 (MQ=255)

                                                      gCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGtt                                   <  1:5345289/93‑1 (MQ=255)

                                                      gCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGtt                                   >  1:5917712/1‑93 (MQ=255)

                                                       ccGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTa                                  <  1:2989982/93‑1 (MQ=255)

                                                         gTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGGTACg                                >  2:3444106/1‑93 (MQ=255)

                                                          ttAGAACGGCGTGCAGTGTCACGCTCGTTTTTGGCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGa                               <  2:5016091/93‑1 (MQ=255)

                                                          ttAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACgt                               >  1:5127512/1‑92 (MQ=255)

                                                          ttAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGa                               <  1:1704595/93‑1 (MQ=255)

                                                           tAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGAt                              <  1:1137530/93‑1 (MQ=255)

                                                           tAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGAt                              <  2:6068589/93‑1 (MQ=255)

                                                            aGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATa                             <  1:2456076/93‑1 (MQ=255)

                                                               aCGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGt                          >  2:3669003/1‑93 (MQ=255)

                                                               aCGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGt                          <  1:481317/93‑1 (MQ=255)

                                                               aCGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGATGGAGATACGATAGGt                          >  2:4329535/1‑93 (MQ=255)

                                                                cggcggGCGGTGTCCCCCTCTTTTTTGCCCAGGAACGAAACAGCAAAGTCAAGGCCACAAACAAGACAAAAGAAGTGGGAGTTACGATAGGTa                         <  2:5831922/93‑1 (MQ=255)

                                                                cGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTa                         <  1:4580796/93‑1 (MQ=255)

                                                                cGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTa                         <  2:1899261/93‑1 (MQ=255)

                                                                  gCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGc                       <  2:5397530/93‑1 (MQ=255)

                                                                  gCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGc                       >  1:5022881/1‑93 (MQ=255)

                                                                     tGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAc                    >  2:3512750/1‑93 (MQ=255)

                                                                      gCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAcg                   >  2:4768496/1‑93 (MQ=255)

                                                                          tgtCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAcgccgc               <  2:1549378/93‑1 (MQ=255)

                                                                          tgtCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAcgccgc               <  1:4606275/93‑1 (MQ=255)

                                                                             cACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCaa            <  2:2986635/93‑1 (MQ=255)

                                                                              aCGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAac           >  1:4213325/1‑93 (MQ=255)

                                                                               cGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAaca          <  1:5732245/93‑1 (MQ=255)

                                                                                gCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAacac         >  1:493979/1‑93 (MQ=255)

                                                                                  tCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGa       >  2:1744730/1‑93 (MQ=255)

                                                                                   cGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGAc      >  1:2057029/1‑93 (MQ=255)

                                                                                     tttttGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCa    <  2:4122995/93‑1 (MQ=255)

                                                                                      ttttGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCaa   >  1:546300/1‑93 (MQ=255)

                                                                                      ttttGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCaa   <  2:1951067/93‑1 (MQ=255)

ACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAAC  >  NC_000913/985573‑985752

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