BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,134,313	A→G	Q121R (CAG→CGG)	flgG →	FlgG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,134,313	0	A	G	83.3%	577.3 / 113.0	263	Q121R (CAG→CGG)	flgG	FlgG
Reads supporting (aligned to +/- strand): ref base A (22/22); new base G (126/93); total (148/115)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 264 total aligned reads displayed.

ACAGCAAAGATGTCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCAT  >  NC_000913/1134223‑1134403

acaGCAAAGATGTCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCggg                                                                                          >  1:1668639/1‑93 (MQ=255)

 caGCAAAGATGTCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGt                                                                                         <  1:195314/93‑1 (MQ=255)

 caGCAAAGATGTCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGt                                                                                         >  2:1577182/1‑93 (MQ=255)

   gCAAAGATGTCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTgg                                                                                       >  1:165924/1‑93 (MQ=255)

   gCAAAGATGTCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTgg                                                                                       <  1:2859152/93‑1 (MQ=255)

    cAAAGATGTCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGa                                                                                      >  1:1398658/1‑93 (MQ=255)

       aGATGTCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCa                                                                                   >  1:4145747/1‑93 (MQ=255)

        gATGTCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAg                                                                                  >  1:3095496/1‑93 (MQ=255)

         aTGTCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGa                                                                                 >  2:6470553/1‑93 (MQ=255)

         aTGTCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGa                                                                                 >  1:1149369/1‑93 (MQ=255)

          tGTCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGaa                                                                                <  2:1619514/93‑1 (MQ=255)

            tCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACg                                                                              <  1:7034937/93‑1 (MQ=255)

             cgcgATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACgg                                                                             <  1:3950852/93‑1 (MQ=255)

               cgATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGc                                                                           <  1:2675357/93‑1 (MQ=255)

               cgATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGc                                                                           <  2:2551486/93‑1 (MQ=255)

                gATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCa                                                                          <  1:4986062/93‑1 (MQ=255)

                 aTTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAg                                                                         >  2:4798183/1‑93 (MQ=255)

                  ttAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGc                                                                        <  1:5603656/93‑1 (MQ=255)

                   tAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCt                                                                       >  1:1404736/1‑93 (MQ=255)

                    aaaGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTg                                                                      >  1:3783509/1‑93 (MQ=255)

                     aaGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTgg                                                                     <  2:6178747/93‑1 (MQ=255)

                     aaGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTgg                                                                     <  1:6493135/93‑1 (MQ=255)

                      agggcagggcTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGt                                                                    >  2:3502677/1‑93 (MQ=255)

                        ggcagggcTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGa                                                                  >  2:6869632/1‑93 (MQ=255)

                         gcagggcTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGGCGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGAc                                                                 <  2:544897/93‑1 (MQ=255)

                           agggcTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACgg                                                               >  2:3280612/1‑93 (MQ=255)

                           agggcTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACgg                                                               >  2:889938/1‑93 (MQ=255)

                           agggcTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACgg                                                               >  1:32258/1‑93 (MQ=255)

                           agggcTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACgg                                                               >  1:6040276/1‑93 (MQ=255)

                           agggcTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACgg                                                               >  2:6746450/1‑93 (MQ=255)

                            gggcTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGc                                                              >  1:2046439/1‑93 (MQ=255)

                             ggcTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCt                                                             >  2:4615612/1‑93 (MQ=255)

                             ggcTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCt                                                             >  2:2251567/1‑93 (MQ=255)

                               cTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCtgg                                                           >  1:594453/1‑93 (MQ=255)

                               cTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCtgg                                                           <  1:6849330/93‑1 (MQ=255)

                                tttttCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCtggt                                                          >  1:5285021/1‑93 (MQ=255)

                                tttttCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCtggt                                                          >  2:4276604/1‑93 (MQ=255)

                                 ttttCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCtggtg                                                         >  2:5360300/1‑93 (MQ=255)

                                  tttCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCtggtgg                                                        <  2:6949544/93‑1 (MQ=255)

                                    tCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGtt                                                      <  1:5659023/93‑1 (MQ=255)

                                    tCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGtt                                                      <  2:7055191/93‑1 (MQ=255)

                                     ccAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGttt                                                     <  2:2902048/93‑1 (MQ=255)

                                      cAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGtttt                                                    >  1:2536924/1‑93 (MQ=255)

                                       aGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTc                                                   <  2:4012485/93‑1 (MQ=255)

                                        ggTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCa                                                  >  2:1041567/1‑93 (MQ=255)

                                        ggTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCa                                                  >  1:725876/1‑93 (MQ=255)

                                          tGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAgg                                                <  2:3432735/93‑1 (MQ=255)

                                          tGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAgg                                                <  2:2142270/93‑1 (MQ=255)

                                          tGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAgg                                                <  2:586898/93‑1 (MQ=255)

                                            aTGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTg                                              >  1:3634178/1‑93 (MQ=255)

                                              gTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGca                                            >  1:5057306/1‑93 (MQ=255)

                                               ttGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGcag                                           >  2:4303315/1‑93 (MQ=255)

                                                tGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGcagc                                          >  1:4600422/1‑93 (MQ=255)

                                                 gCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGcagcc                                         <  1:6249445/93‑1 (MQ=255)

                                                 gCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGcagcc                                         >  2:779494/1‑93 (MQ=255)

                                                  ccAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGcagcca                                        <  1:7141963/93‑1 (MQ=255)

                                                  ccAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGcagcca                                        <  2:2873850/93‑1 (MQ=255)

                                                  ccAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGcagcca                                        <  2:662043/93‑1 (MQ=255)

                                                   cAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGcagccag                                       >  1:5969363/1‑93 (MQ=255)

                                                    aGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGcagccagc                                      <  1:849665/93‑1 (MQ=255)

                                                       tGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGAt                                   <  1:6614011/93‑1 (MQ=255)

                                                         gTTCATCAGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCa                                 <  2:6837262/93‑1 (MQ=255)

                                                            catcaGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCa                              <  1:1485039/93‑1 (MQ=255)

                                                             atcaGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCAt                             <  2:5712782/93‑1 (MQ=255)

                                                             atcaGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCAt                             >  1:2207238/1‑93 (MQ=255)

                                                             atcaGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCAt                             <  1:4406742/93‑1 (MQ=255)

                                                             atcaGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCAt                             <  1:2975814/93‑1 (MQ=255)

                                                              tcaGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCAtt                            >  2:1880227/1‑93 (MQ=255)

                                                              tcaGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCAtt                            >  1:2446250/1‑93 (MQ=255)

                                                              tcaGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCAtt                            >  2:1378170/1‑93 (MQ=255)

                                                              tcaGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCAtt                            <  2:7422494/93‑1 (MQ=255)

                                                               caGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTc                           <  1:1225885/93‑1 (MQ=255)

                                                               caGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTc                           <  2:1197933/93‑1 (MQ=255)

                                                               caGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTc                           <  2:1413812/93‑1 (MQ=255)

                                                                aGCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTcc                          <  2:3004509/93‑1 (MQ=255)

                                                                 gCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCg                         >  2:1146971/1‑93 (MQ=255)

                                                                 gCCTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCg                         <  1:914425/93‑1 (MQ=255)

                                                                  ccTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCgg                        >  1:5099036/1‑93 (MQ=255)

                                                                   cTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGc                       >  1:5307962/1‑93 (MQ=255)

                                                                   cTATACCCGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGc                       >  1:783879/1‑93 (MQ=255)

                                                                    tataCCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCg                      >  1:3546778/1‑93 (MQ=255)

                                                                       aCCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAAt                   <  2:4856240/93‑1 (MQ=255)

                                                                        cccGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATg                  <  1:3388235/93‑1 (MQ=255)

                                                                         ccGTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGc                 >  1:5385462/1‑93 (MQ=255)

                                                                         ccGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGc                 <  2:5184728/93‑1 (MQ=255)

                                                                          cGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCg                >  2:3696586/1‑93 (MQ=255)

                                                                           gTGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGt               <  1:1705252/93‑1 (MQ=255)

                                                                            tGACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGtt              >  2:170568/1‑93 (MQ=255)

                                                                             gACGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTa             >  1:1845998/1‑93 (MQ=255)

                                                                              aCGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTGCGGCGAATGCGTTaa            >  2:3857947/1‑93 (MQ=255)

                                                                              aCGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTaa            >  2:7487373/1‑93 (MQ=255)

                                                                              aCGGCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTaa            <  1:6945951/93‑1 (MQ=255)

                                                                                ggCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGt          >  1:2277220/1‑93 (MQ=255)

                                                                                ggCTCTTTCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGt          >  2:7062040/1‑93 (MQ=255)

                                                                                     tttCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCAc     >  1:5747841/1‑93 (MQ=255)

                                                                                     tttCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCAc     >  2:267088/1‑93 (MQ=255)

                                                                                     tttCCGGGTCGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCAc     >  1:3140933/1‑93 (MQ=255)

                                                                                     tttCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCAc     >  1:5458727/1‑93 (MQ=255)

                                                                                     tttCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCAc     >  1:3696634/1‑93 (MQ=255)

                                                                                      ttCCGGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCAcc    <  2:1933935/93‑1 (MQ=255)

ACAGCAAAGATGTCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCAT  >  NC_000913/1134223‑1134403

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