BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,185,099	T→C	S40S (AGT→AGC)	b2107 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,185,099	0	T	C	99.7%	1089.2 / NA	331	S40S (AGT→AGC)	b2107	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (155/175); minor base G (1/0); total (156/175)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 344 total aligned reads displayed.

GTGCGCTTTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCA  >  NC_000913/2185008‑2185188

gTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAagcag                                                                                         <  2:5764247/93‑1 (MQ=255)

 tGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAagcagt                                                                                        <  2:1713711/93‑1 (MQ=255)

   gcgctttGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGtt                                                                                       >  1:7688771/5‑93 (MQ=255)

   gcgctttGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGtt                                                                                       <  1:3444706/89‑1 (MQ=255)

    cgctttGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGttt                                                                                      <  2:9417/90‑1 (MQ=255)

    cgctttGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGttt                                                                                      >  2:4182979/4‑93 (MQ=255)

      ctttGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCa                                                                                    <  2:3632138/92‑1 (MQ=255)

       tttGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCAt                                                                                   <  2:1312682/93‑1 (MQ=255)

         tGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTg                                                                                 >  1:1942540/1‑93 (MQ=255)

          gCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGc                                                                                >  2:6850481/1‑93 (MQ=255)

          gCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGc                                                                                <  2:479868/93‑1 (MQ=255)

            tGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTg                                                                              >  1:8591851/1‑93 (MQ=255)

            tGGTTACCAGTGCCGCCTGGGCCGCACCAGACACCGCGGGTTCGACCAATACCTCGGGAGTTTCTAAGTATGAGTTAAGAAGTTTCATTGCTg                                                                              >  1:6034389/1‑93 (MQ=255)

             ggTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGACTTTCTAAGTATGAGTTAAGCAGTTTCATTGCTgc                                                                             >  2:1003870/1‑92 (MQ=255)

              gTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGAc                                                                            >  1:5438469/1‑93 (MQ=255)

               ttACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACt                                                                           >  1:2253214/1‑93 (MQ=255)

               ttACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACt                                                                           >  2:2975339/1‑93 (MQ=255)

                tACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACtt                                                                          <  1:1664234/93‑1 (MQ=255)

                   cAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTaa                                                                       >  1:8868436/1‑93 (MQ=255)

                   cAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTaa                                                                       <  2:7151141/93‑1 (MQ=255)

                   cAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGAATTTaa                                                                       >  2:3942419/1‑93 (MQ=255)

                     gTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGc                                                                     <  2:1940603/93‑1 (MQ=255)

                      tGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCa                                                                    <  2:9234037/93‑1 (MQ=255)

                       gccgccTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCAt                                                                   <  2:8502282/93‑1 (MQ=255)

                       gccgccTGGGCCGCACCAGCCACCGCGGGTTCGACCAATAACTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCAt                                                                   >  1:5688440/1‑93 (MQ=255)

                         cgccTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCAttt                                                                 >  1:3701788/1‑93 (MQ=255)

                          gccTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTc                                                                <  2:3734763/93‑1 (MQ=255)

                          cgcTGGGCCGCCCCAGCCACCGCGGGTTCGACCAATGCCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTc                                                                <  1:8695927/91‑1 (MQ=255)

                             tGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCaaa                                                             <  1:3564267/93‑1 (MQ=255)

                             tGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCaaa                                                             >  1:6854568/1‑93 (MQ=255)

                              gggCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAAc                                                            <  1:9199726/93‑1 (MQ=255)

                              gggCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAAc                                                            >  1:4606179/1‑93 (MQ=255)

                              gggCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAAc                                                            >  2:6073781/1‑93 (MQ=255)

                               ggCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAAcc                                                           <  1:7953770/93‑1 (MQ=255)

                                gCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCa                                                          <  2:8604107/93‑1 (MQ=255)

                                gCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCa                                                          <  2:6375098/93‑1 (MQ=255)

                                gCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCa                                                          >  2:2779949/1‑93 (MQ=255)

                                  cGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAgg                                                        >  2:3042311/1‑93 (MQ=255)

                                   gCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAggg                                                       >  1:7276878/1‑93 (MQ=255)

                                   gCACCAGCCACCGCGAGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAggg                                                       <  1:7808054/93‑1 (MQ=255)

                                    cACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAgggg                                                      <  2:6690248/93‑1 (MQ=255)

                                     cccAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGa                                                     <  1:4080301/92‑1 (MQ=255)

                                      ccAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGac                                                    >  1:9296983/1‑93 (MQ=255)

                                      ccAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGac                                                    >  1:6973233/1‑93 (MQ=255)

                                      ccAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGac                                                    >  2:6527067/1‑93 (MQ=255)

                                       cAGCCACCGCGGGCTCGACCAATACCTCGGGCATTTCTCAGTATGAGTTCAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGaca                                                   >  2:5174331/1‑93 (MQ=255)

                                         gCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACAcc                                                 >  2:8771088/1‑93 (MQ=255)

                                          ccACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCg                                                >  1:1207581/1‑93 (MQ=255)

                                          ccACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCg                                                <  2:8204437/93‑1 (MQ=255)

                                          ccACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCg                                                <  1:5338472/93‑1 (MQ=255)

                                          ccACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCg                                                <  1:1409791/93‑1 (MQ=255)

                                           cACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGt                                               >  1:1347681/1‑93 (MQ=255)

                                            aCCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTa                                              >  1:1012287/1‑93 (MQ=255)

                                             ccGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTAc                                             <  1:9929108/93‑1 (MQ=255)

                                             ccGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTAc                                             >  2:1482015/1‑93 (MQ=255)

                                               gcgGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCa                                           <  2:3255222/93‑1 (MQ=255)

                                                cgGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAg                                          <  2:8376808/93‑1 (MQ=255)

                                                  ggTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGaa                                        >  1:1769646/1‑93 (MQ=255)

                                                  ggTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGaa                                        >  1:5033508/1‑93 (MQ=255)

                                                   gTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGaaa                                       >  1:7503209/1‑93 (MQ=255)

                                                      cGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgt                                    <  2:7478563/93‑1 (MQ=255)

                                                        aCCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtac                                  >  2:2259975/1‑93 (MQ=255)

                                                        aCCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtac                                  <  1:3941174/93‑1 (MQ=255)

                                                         ccAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtacc                                 <  2:8003902/93‑1 (MQ=255)

                                                          cAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtaccg                                >  2:7270586/1‑93 (MQ=255)

                                                          cAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtaccg                                <  1:826795/93‑1 (MQ=255)

                                                           aaTACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtaccgt                               >  1:4729100/1‑93 (MQ=255)

                                                           aaTACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtaccgt                               >  1:2781524/1‑93 (MQ=255)

                                                             tACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtaccgtac                             <  1:2454137/93‑1 (MQ=255)

                                                             tACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtaccgtac                             >  2:5917615/1‑93 (MQ=255)

                                                              aCCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtaccgtacc                            <  2:786493/93‑1 (MQ=255)

                                                              aCCTCGGGAATTTATAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtaccgtacc                            <  1:6093545/93‑1 (MQ=255)

                                                                cTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGa                          <  2:4632361/93‑1 (MQ=255)

                                                                cTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGa                          <  2:9810333/93‑1 (MQ=255)

                                                                   gggAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTTAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGa                       >  2:8176355/1‑93 (MQ=255)

                                                                       aTTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTac                   <  2:4958693/93‑1 (MQ=255)

                                                                       aTTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGCACCGTACCGATGAGTac                   >  1:8993045/1‑93 (MQ=255)

                                                                        tttCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCTTTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTaca                  >  2:9632810/1‑93 (MQ=255)

                                                                        tttCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTaca                  >  2:911367/1‑93 (MQ=255)

                                                                        tttCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTaca                  <  1:1750055/93‑1 (MQ=255)

                                                                          tCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTacaac                <  2:5574787/93‑1 (MQ=255)

                                                                          tCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTacaac                <  1:4274964/93‑1 (MQ=255)

                                                                          tCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTacaac                <  1:6338063/93‑1 (MQ=255)

                                                                          tCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTacaaa                >  1:6765893/1‑92 (MQ=255)

                                                                           cTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTacaaca               <  1:6628860/93‑1 (MQ=255)

                                                                            tAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACAt              >  1:5910064/1‑93 (MQ=255)

                                                                             aaGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACAtt             >  2:2648992/1‑93 (MQ=255)

                                                                              aGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTa            >  1:5655549/1‑93 (MQ=255)

                                                                              aGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTa            <  1:518906/93‑1 (MQ=255)

                                                                                 aTGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGc         <  2:2090432/93‑1 (MQ=255)

                                                                                   gAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAg       <  1:297166/93‑1 (MQ=255)

                                                                                   gAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAg       <  1:7095351/93‑1 (MQ=255)

                                                                                     gTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTg     >  2:693042/1‑93 (MQ=255)

                                                                                     gTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTg     <  1:5472751/93‑1 (MQ=255)

                                                                                     gTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAAGACAACATAAAGCAggg     >  1:2305888/1‑91 (MQ=255)

                                                                                      ttAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTgg    <  1:9352157/93‑1 (MQ=255)

                                                                                       tAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGc   <  1:2052863/93‑1 (MQ=255)

                                                                                       tAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGc   >  1:6218236/1‑93 (MQ=255)

                                                                                        aaGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCa  <  2:1636948/93‑1 (MQ=255)

                                                                                        aaGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCa  >  1:2674195/1‑93 (MQ=255)

GTGCGCTTTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCA  >  NC_000913/2185008‑2185188

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