BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,890,598	A→G	V221A (GTC→GCC)	b1808 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,890,598	0	A	G	86.2%	557.2 / 74.0	231	V221A (GTC→GCC)	b1808	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (19/13); new base G (110/89); total (129/102)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 232 total aligned reads displayed.

GCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCG  >  NC_000913/1890508‑1890690

 ccAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGt                                                                                           <  2:2341524/93‑1 (MQ=255)

   aGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTcc                                                                                         <  2:4866339/93‑1 (MQ=255)

   aGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTcc                                                                                         <  1:2544542/93‑1 (MQ=255)

    ggTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCg                                                                                        >  2:4544131/1‑93 (MQ=255)

    ggTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCg                                                                                        >  1:7085264/1‑93 (MQ=255)

     gTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGc                                                                                       >  1:2996154/1‑93 (MQ=255)

       cGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCtt                                                                                     <  2:6133237/93‑1 (MQ=255)

        gAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTc                                                                                    >  2:110276/1‑93 (MQ=255)

         aGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTcc                                                                                   >  2:315642/1‑93 (MQ=255)

         aGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTcc                                                                                   <  2:3860440/93‑1 (MQ=255)

         aGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTcc                                                                                   <  1:2403689/93‑1 (MQ=255)

         aGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTcc                                                                                   >  2:2090772/1‑93 (MQ=255)

          gCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCg                                                                                  >  2:4249594/1‑93 (MQ=255)

          gCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCg                                                                                  >  2:1732825/1‑93 (MQ=255)

           cAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCgg                                                                                 <  2:7284026/93‑1 (MQ=255)

           cAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCgg                                                                                 >  1:1925376/1‑93 (MQ=255)

            aGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCggg                                                                                >  1:5825225/1‑93 (MQ=255)

            aGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCggg                                                                                >  2:5931569/1‑93 (MQ=255)

             gTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGa                                                                               <  2:1068085/93‑1 (MQ=255)

             gTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGa                                                                               >  1:1498795/1‑93 (MQ=255)

             gTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGa                                                                               >  2:2482329/1‑93 (MQ=255)

              ttGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGAt                                                                              >  2:5674357/1‑93 (MQ=255)

               tGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATc                                                                             <  1:5207822/93‑1 (MQ=255)

               tGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATc                                                                             >  2:635450/1‑93 (MQ=255)

                gTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCa                                                                            >  2:2919217/1‑93 (MQ=255)

                 tCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAg                                                                           <  2:1174539/93‑1 (MQ=255)

                 tCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAg                                                                           <  2:462021/93‑1 (MQ=255)

                   gACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGtt                                                                         >  1:1202731/1‑93 (MQ=255)

                   gACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGtt                                                                         >  1:2003194/1‑93 (MQ=255)

                    aCTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTc                                                                        >  1:469729/1‑93 (MQ=255)

                     cTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCg                                                                       <  1:5161736/93‑1 (MQ=255)

                      tGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGc                                                                      >  2:1981734/1‑93 (MQ=255)

                       ggAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGcc                                                                     <  1:7296429/93‑1 (MQ=255)

                       ggAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGcc                                                                     <  2:2158306/93‑1 (MQ=255)

                        gagaGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCa                                                                    >  2:6499525/1‑93 (MQ=255)

                           agTGACTGACCATAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAAt                                                                 <  2:6162392/93‑1 (MQ=255)

                           agTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAAt                                                                 <  2:5399509/93‑1 (MQ=255)

                            gTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATc                                                                >  2:3472426/1‑93 (MQ=255)

                            gTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATc                                                                <  2:575767/93‑1 (MQ=255)

                            gTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATc                                                                >  1:4879827/1‑93 (MQ=255)

                            gTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATc                                                                <  2:6798449/93‑1 (MQ=255)

                             tgactgacCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATcc                                                               <  1:5642142/93‑1 (MQ=255)

                             tgactgacCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATcc                                                               >  2:2574516/1‑93 (MQ=255)

                                ctgacCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCAct                                                            >  2:4100373/1‑93 (MQ=255)

                                 tgacCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCActc                                                           >  1:1378920/1‑93 (MQ=255)

                                  gacCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCActct                                                          >  2:1346728/1‑93 (MQ=255)

                                  gacCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCActct                                                          <  2:7177875/93‑1 (MQ=255)

                                  gacCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCActct                                                          >  2:4469178/1‑93 (MQ=255)

                                     cAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTTCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCtt                                                       <  1:5983031/93‑1 (MQ=255)

                                     cAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCtt                                                       <  1:3856310/93‑1 (MQ=255)

                                     cAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCtt                                                       >  2:6285095/1‑93 (MQ=255)

                                      aaaaTACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCttt                                                      >  1:5301455/1‑93 (MQ=255)

                                        aaTACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTaa                                                    <  2:5200905/93‑1 (MQ=255)

                                        aaTACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTaa                                                    <  1:4987893/93‑1 (MQ=255)

                                            ctggctggCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAaccac                                                <  2:6649766/93‑1 (MQ=255)

                                            ctggctggCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAaccac                                                >  2:3101511/1‑93 (MQ=255)

                                                ctggcAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCata                                            >  2:5498161/1‑93 (MQ=255)

                                                ctggcAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCata                                            <  1:7146087/93‑1 (MQ=255)

                                                 tggcAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCatat                                           <  1:1733989/93‑1 (MQ=255)

                                                 tggcAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCatat                                           >  1:3738955/1‑93 (MQ=255)

                                                 tggcAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCatat                                           >  2:6999380/1‑93 (MQ=255)

                                                  ggcAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATc                                          <  1:6679127/93‑1 (MQ=255)

                                                     aaTGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACGATATCCGc                                       >  2:3335389/1‑93 (MQ=255)

                                                     aaTGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGc                                       >  1:2762445/1‑93 (MQ=255)

                                                      aTGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGcc                                      >  1:1007183/1‑93 (MQ=255)

                                                      aTGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGcc                                      >  1:455800/1‑93 (MQ=255)

                                                      aTGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGcc                                      >  1:3346035/1‑93 (MQ=255)

                                                       tGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCa                                     >  2:6718357/1‑93 (MQ=255)

                                                       tGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCa                                     >  2:1253648/1‑93 (MQ=255)

                                                           cGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCCCTCTCTTTAACCACCATATCCGCCAGaaa                                 <  1:4441206/93‑1 (MQ=255)

                                                           cGGTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGaaa                                 >  1:6473908/1‑93 (MQ=255)

                                                            ggTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAag                                <  2:6899846/93‑1 (MQ=255)

                                                            ggTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAag                                <  1:2540052/93‑1 (MQ=255)

                                                             gTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAaga                               >  2:953351/1‑93 (MQ=255)

                                                             gTAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAaga                               <  2:1756838/93‑1 (MQ=255)

                                                             gTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAaga                               <  2:1832698/93‑1 (MQ=255)

                                                             gTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAaga                               >  2:6318726/1‑93 (MQ=255)

                                                              tAGCTGGTGGGCTTTGTCGAAGATCATGGCGTCCGCTTTCGGGGTCAGGTCACCAAATCCACTCTCTTTAACCACCATATCCGCCAGCAagag                              >  1:1856455/1‑93 (MQ=255)

                                                              tAGCTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAagag                              >  2:4943988/1‑93 (MQ=255)

                                                                 cTGGTGGGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAgatg                           <  1:3006533/93‑1 (MQ=255)

                                                                  tggtggGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAgatga                          <  1:419360/93‑1 (MQ=255)

                                                                   ggtggGCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAgatgat                         <  2:5333659/93‑1 (MQ=255)

                                                                      gggCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGt                      >  1:5115835/1‑93 (MQ=255)

                                                                       ggCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTAAACCACCATATCCGCCAGAAAGAGATGATGGtt                     >  2:6877725/1‑93 (MQ=255)

                                                                        gCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGttt                    <  1:104729/93‑1 (MQ=255)

                                                                        gCTTCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGttt                    <  1:5908712/93‑1 (MQ=255)

                                                                        gCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGttt                    >  1:6175641/1‑93 (MQ=255)

                                                                          ttCGTCGAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTac                  >  1:6021233/1‑93 (MQ=255)

                                                                             gtcgAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTaccac               >  1:3133628/1‑93 (MQ=255)

                                                                             gtcgAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTaccac               >  2:823259/1‑93 (MQ=255)

                                                                                gAAGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTaccaccac            <  1:1276267/93‑1 (MQ=255)

                                                                                 aaGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTaccaccacc           <  2:4987975/93‑1 (MQ=255)

                                                                                  aGATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTAccaccacca          >  2:1429421/1‑93 (MQ=255)

                                                                                   gATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACcaccaccac         <  2:3913068/93‑1 (MQ=255)

                                                                                   gATCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACcaccaccac         >  2:7470493/1‑93 (MQ=255)

                                                                                    aTCATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACcaca        >  1:3632198/1‑93 (MQ=255)

                                                                                      cATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATc      <  1:624992/93‑1 (MQ=255)

                                                                                      cATGGCGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATc      >  2:352233/1‑93 (MQ=255)

                                                                                      cATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATc      <  1:6727669/93‑1 (MQ=255)

                                                                                         gACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGc   >  2:7230983/1‑93 (MQ=255)

GCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCG  >  NC_000913/1890508‑1890690

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